Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,654,036 C→T A14T (GCA→ACA)  unknown ← RNA polymerase sigma factor RpoD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,654,0360CT100.0% 230.9 / NA 79A14T (GCA→ACA) unknownRNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/45);  total (34/45)

CGGTAGTCCGCCCTAATTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGGCTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGCGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGTTGCAGGTAGAAATAACG  >  1/1653790‑1654282
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cGGTAGTCCGCCCTAATTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATg                                                                                                                                                                                                                                                    <  1:144074/251‑1 (MQ=255)
   tAGTCCGCCCTAATTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGagt                                                                                                                                                                                                                                                 <  1:805978/251‑1 (MQ=255)
      tCCGCCCTAATTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGagtag                                                                                                                                                                                                                                               <  1:492200/250‑1 (MQ=255)
              aaTTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGAGCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTg                                                                                                                                                                                                                                      <  1:425110/251‑1 (MQ=255)
               aTTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGttt                                                                                                                                                                                                                                            <  1:507755/244‑1 (MQ=255)
                  cTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGAc                                                                                                                                                                                                                                  <  1:614737/251‑1 (MQ=255)
                          aCGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGAtt                                                                                                                                                                                                                          >  1:691585/1‑251 (MQ=255)
                              tGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGa                                                                                                                                                                                                                      <  1:598996/251‑1 (MQ=255)
                               ggATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTg                                                                                                                                                                                                                                      >  1:187555/1‑234 (MQ=255)
                               ggATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGAt                                                                                                                                                                                                                     >  1:577623/1‑251 (MQ=255)
                               ggATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGCGGATTTTGAt                                                                                                                                                                                                                     >  1:29117/1‑251 (MQ=255)
                               ggATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCAGCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTgaga                                                                                                                                                                                                                                   <  1:432768/237‑1 (MQ=255)
                                     taCTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCatat                                                                                                                                                                                                               >  1:185872/1‑251 (MQ=255)
                                     taCTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCatat                                                                                                                                                                                                               >  1:156083/1‑251 (MQ=255)
                                      aCTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGag                                                                                                                                                                                                                                                  >  1:645526/1‑215 (MQ=255)
                                       cTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATAttt                                                                                                                                                                                                             <  1:82929/251‑1 (MQ=255)
                                           ggTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:171184/250‑1 (MQ=255)
                                           ggTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:244763/250‑1 (MQ=255)
                                           ggTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:625346/250‑1 (MQ=255)
                                           ggTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:28392/250‑1 (MQ=255)
                                           ggTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGACTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:401542/250‑1 (MQ=255)
                                             ttGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATc                                                                                                                                                                                                       <  1:864724/251‑1 (MQ=255)
                                                            cATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGAtt                                                                                                                                                                                                                          >  1:161394/1‑217 (MQ=255)
                                                            cATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTa                                                                                                                                                                                        <  1:37943/251‑1 (MQ=255)
                                                               cTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTAttt                                                                                                                                                                                     >  1:175276/1‑251 (MQ=255)
                                                                                     gAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGa                                                                                                                                                               >  1:290697/1‑251 (MQ=255)
                                                                                            aGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGAt                                                                                                                                                                            <  1:7582/231‑1 (MQ=255)
                                                                                              atcatcGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTat                                                                                                                                                         <  1:854773/248‑1 (MQ=255)
                                                                                                       aTCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:391164/190‑1 (MQ=255)
                                                                                                               cGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAAc                                                                                                                                     <  1:788849/251‑1 (MQ=255)
                                                                                                                 gTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACt                                                                                                                                    <  1:48747/250‑1 (MQ=255)
                                                                                                                   ctcCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAat                                                                                                                                                                                                          <  1:873935/178‑1 (MQ=255)
                                                                                                                   ctcCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTAttt                                                                                                                                                                                     >  1:617211/1‑199 (MQ=255)
                                                                                                                   ctcCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTat                                                                                                                                                         >  1:406793/1‑227 (MQ=255)
                                                                                                                   ctcCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACAt                                                                                                                                             >  1:8851/1‑239 (MQ=255)
                                                                                                                           ccTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTa                                                                                                                                                                                                            <  1:770869/168‑1 (MQ=255)
                                                                                                                              gAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACt                                                                                                                                    >  1:181961/1‑237 (MQ=255)
                                                                                                                                         tCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGaaa                                                                                                                                                                                >  1:678084/1‑182 (MQ=255)
                                                                                                                                           aTTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACa                                                                                                         <  1:864639/251‑1 (MQ=255)
                                                                                                                                            ttAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAAt                                                                                                                                                                       <  1:22222/188‑1 (MQ=255)
                                                                                                                                            ttAATCATCTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACaa                                                                                                        >  1:454973/1‑251 (MQ=255)
                                                                                                                                                 catcTGAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAAtgtg                                                                                                   >  1:126521/1‑251 (MQ=255)
                                                                                                                                                      gAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGttt                                                                                                                                                                                                                                            <  1:506887/109‑1 (MQ=255)
                                                                                                                                                      gAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGAc                                                                                                                                                                                                                                  <  1:593372/119‑1 (MQ=255)
                                                                                                                                                      gAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATc                                                                                                                                                                                                       <  1:717126/146‑1 (MQ=255)
                                                                                                                                                      gAATAATATCTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTa                                                                                                                                   <  1:509762/214‑1 (MQ=255)
                                                                                                                                                               cTTCAATTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATc                                                                                                                                                                                                                    <  1:3180/124‑1 (MQ=255)
                                                                                                                                                                    aTTTGATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAAc                                                                                                                                     <  1:602073/198‑1 (MQ=255)
                                                                                                                                                                        gATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGa                                                                                                  <  1:643925/229‑1 (MQ=255)
                                                                                                                                                                        gATCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGa                                                                            >  1:232225/1‑251 (MQ=255)
                                                                                                                                                                        aaTCTGCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGa                                                                            >  1:232202/2‑251 (MQ=255)
                                                                                                                                                                             gCATCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTAc                                                                                                              <  1:239809/212‑1 (MQ=255)
                                                                                                                                                                               aTCGACTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGc                                                                     <  1:63421/251‑1 (MQ=255)
                                                                                                                                                                                  gacTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGAt                                                                                                                                                                            >  1:546633/3‑145 (MQ=255)
                                                                                                                                                                                        ttcttcGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCt                                                                   <  1:544393/244‑1 (MQ=255)
                                                                                                                                                                                             cGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGaaa                                                                                                                                                                                >  1:46364/1‑130 (MQ=255)
                                                                                                                                                                                                ggAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCttt                                                     <  1:888255/250‑1 (MQ=255)
                                                                                                                                                                                                    gCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGAGTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGc                                                  <  1:423492/249‑1 (MQ=255)
                                                                                                                                                                                                     cTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGa                                                <  1:607782/250‑1 (MQ=255)
                                                                                                                                                                                                       gTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGa                                             >  1:710816/1‑251 (MQ=255)
                                                                                                                                                                                                         cATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTAtctc                                                                                                                               >  1:406652/1‑167 (MQ=255)
                                                                                                                                                                                                          aTTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGttt                                                                                                                                                                                                                                            <  1:268897/57‑1 (MQ=255)
                                                                                                                                                                                                                  cGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTAtt                                  >  1:13701/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCg                                <  1:432492/250‑1 (MQ=255)
                                                                                                                                                                                                                         ggTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTAt                           >  1:117679/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTAt                           >  1:422180/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTAt                           >  1:101525/1‑251 (MQ=255)
                                                                                                                                                                                                                          gTTAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTa                                                                                                                                   >  1:514119/1‑146 (MQ=255)
                                                                                                                                                                                                                            tAAATAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTAGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACa                                                                                                                                              <  1:416813/133‑1 (MQ=255)
                                                                                                                                                                                                                             aaaTAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTAc                                                                                                              >  1:853211/1‑164 (MQ=255)
                                                                                                                                                                                                                             aaaTAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACt                                                                                                             <  1:739225/165‑1 (MQ=255)
                                                                                                                                                                                                                             aaaTAACCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCt                                                                   <  1:845000/207‑1 (MQ=255)
                                                                                                                                                                                                                                  aCCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATc                                                                                                                                                                                                       >  1:667267/1‑70 (MQ=255)
                                                                                                                                                                                                                                  aCCTTGCTCTTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGtt                  >  1:890161/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ctTTCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGc                                                  >  1:402060/1‑211 (MQ=255)
                                                                                                                                                                                                                                             tCCCTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGTTGCAGGTAGaaa       >  1:479782/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cccTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTa                                                                                                                                   >  1:713859/1‑126 (MQ=255)
                                                                                                                                                                                                                                              cccTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGt                               <  1:338401/226‑1 (MQ=255)
                                                                                                                                                                                                                                              cccTTGTGTGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGTTGCAGGTAg          <  1:433386/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  tgtgtgATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGTTGCAGGTAGAAATAACg  <  1:891056/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CGGTAGTCCGCCCTAATTCTGCTTCCACGCTGGATAATACTTGGGTTGCTTCTTCAACGACATCTTCATCGGTGATATTTTCATTGAGCATCAGATCATCGGCATCCGGTGCGGTCTCCAACACCTGAATCCCCATATCATTAATCATCTGAATAATATCTTCAATTTGATCTGCATCGGCTAATTCTTCGGGGAGGCTGTCATTGACCTCGGCGTAGGTTAAATAACCTTGCTCTTTCCCTTGTGCGATGAGTAGTTTCAATTGAGACTGTGGATTTTGATCCATATTTAATATCCGCCTTGCTTGCTTATTTGGAAAAGATAAAATGCTGGGGATTTTATCATCGTTAACATGGTTTAACTATCTCTTTTCTTGTTTTTTTACTAACAAATGTGAAAGTTCTTTTTTCTCGTTCAGATCCAAGCCTTCCGCGCGATCTTTAGCGATGAGCCAGTTTATTCGTTTATCAACGAGTTGCAGGTAGAAATAACG  >  1/1653790‑1654282

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑