Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,522,272 A→C 22.4% Q80P (CAG→CCG)  ftsZ → Cell division protein FtsZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,522,2720AC22.4% 144.4 / 5.2 88Q80P (CAG→CCG) ftsZCell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base A (52/11);  new base C (3/19);  total (55/33)
Fisher's exact test for biased strand distribution p-value = 1.17e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCAAGGAA  >  ACB122/1522138‑1522402
                                                                                                                                        |                                                                                                                                  
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  3:178534/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  1:264124/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  6:124980/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  2:89871/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  6:326377/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGGCa                                                                                                                                    >  5:127767/1‑135 (MQ=255)
cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCACCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    >  5:137128/1‑135 (MQ=255)
  aaGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAgg                                                                                                                                  >  8:193388/1‑135 (MQ=255)
  aaGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTCCAATTGGGCCCGGGTGTGACCATGGGCCTGGGGGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                                                                                                    <  2:322668/133‑1 (MQ=255)
      cGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCcgc                                                                                                                              >  2:134492/1‑135 (MQ=255)
         ggAATTCATCTGCGCAAACACGGATGCCCAGGCGCTTAAAAACTTTGGCGCGCGCACAATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGGGCCGGTGCCAATCCGTAAGGCGGCCGTCAGGCCGCGCt                                                                                                                           <  3:446815/135‑1 (MQ=255)
          gAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCt                                                                                                                           >  3:178388/1‑134 (MQ=255)
          gAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAAATGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTAGGCGGGCACGCCgcgcgg                                                                                                                          >  5:347799/1‑133 (MQ=255)
           aaTTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCgcgcgg                                                                                                                          >  1:366708/1‑132 (MQ=255)
           aaTTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTg                                                                                                                          >  1:133749/1‑134 (MQ=255)
             ttCAGCTGGCCCAACCCCGATGCCCAGGCGCTTAAAAACATGGGCGCGCGAACCATCCTGAATTTGGGCAGGGGTGTTCCCACGGGCCTGGGGGCGGGGGCCAATCCGGACGTCGGCCGTCAGGCCGCGCTGGaa                                                                                                                       <  8:195936/135‑1 (MQ=255)
                        cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAgc                                                                                                             >  4:409424/1‑134 (MQ=255)
                        cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAgc                                                                                                             >  5:195247/1‑134 (MQ=255)
                        cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAgc                                                                                                             >  3:355080/1‑134 (MQ=255)
                             ccGATCCCCCGGCGCTCAAAAACTTTGGCGCGCGCACCACCCTGCATTGGGGCACGGGTGTCACCACGGGCCGGGGGGCCGGTGCCAATCCGTAAGGCGGCCGTCTGGCCGCGCTGGAAGACCGTGAGCGCATTg                                                                                                       <  7:307137/135‑1 (MQ=255)
                                agccccGGGCGCTTAAAAACTTGGGCGCGCCCACCATCGTCAATTGGGGCCCGGGTGTGACCATGGGCCGGGGGGCGGGTGCAAATCGGTAAGTGGGCCGACTGGCCGCGCTGGAAGACCGTGAGCGCATTGCCg                                                                                                    <  6:364295/132‑1 (MQ=255)
                                  gCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGaa                                                                                                  >  8:307137/1‑135 (MQ=255)
                                  gCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGaa                                                                                                  >  4:81699/1‑135 (MQ=255)
                                  gCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCACTCCGCGCTGGAAGCCCCTGCGCGCATTGCCGaa                                                                                                  >  7:210304/1‑135 (MQ=255)
                                  gCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCCGCCCGCCGCGCGGGCAGACCCTGCGCGCAATTCCGaa                                                                                                  >  1:444896/1‑135 (MQ=255)
                                   cccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCCGGCCGCGCTGGAAGCCCGGGAGCGCATTGCCGAAg                                                                                                 >  3:64622/1‑135 (MQ=255)
                                   cccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCCGGCCGCCCGGGCAGCCCGTGACCCCAGTGCCGAAg                                                                                                 >  6:249643/1‑135 (MQ=255)
                                   cccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCCTTGCCGCAg                                                                                                 >  6:51159/1‑135 (MQ=255)
                                          cctcaaaaaCATGGGCGCGCGCACCATCCTCCAATTGGCCCCGGGGGGGCCCACGGGCCGGGGTGCCGGGGCAAATCCGCACGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCa                                                                                          <  5:196277/131‑1 (MQ=255)
                                           cTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAg                                                                                         >  3:205820/1‑135 (MQ=255)
                                               aaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGGGTGACCAAGGGCCTGGGGGCCGGTGCCAATCCGGACGTCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAggg                                                                                       <  7:104791/133‑1 (MQ=255)
                                               aaaaCATTGGCGCGCGCACCATCCTGCAATGGGGCAGGGGTGTGACCACGGGCGGGGGGGCGGGTGCAAATCCGGAAGCGGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCa                                                                                     <  4:317932/135‑1 (MQ=255)
                                               aaaaCATTGGCGCGCGCACCAGCCTGCAATGGGCCACGGGGGTGACCACGGGCCTGGGGGCCGGGGCCAATCCGCACGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAggg                                                                                       <  4:234283/133‑1 (MQ=255)
                                                aaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCAATGCCGAAGTGCTGCAGGGccc                                                                                    >  3:455842/1‑133 (MQ=255)
                                                 aaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCAc                                                                                    >  7:15106/1‑134 (MQ=255)
                                                  acttgggCGCGCGCACCATCCTGCATTTGGGCACGGGGGGGACCACGGGCCGGGGGGCGGGGGCAAATCGGTACGGCGGCCGGCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCa                                                                                  <  4:430752/130‑1 (MQ=255)
                                                    aTTGGCGCGCGCACCATCGTGCAATGGGCCAGGGGTGTCACCAAGGGCCGGGGGGCGGGTGCAAATCCGGAAGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAAc                                                                                <  7:110226/135‑1 (MQ=255)
                                                    aTTGGCCGGCGCACCACCCTCAATTGGGGCCGGGGGGTGACAATGGGCGTGGGGGCGGGGGCCAAGCCGGACGGCGCCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAAc                                                                                <  3:294014/135‑1 (MQ=255)
                                                         cgcgcgcACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCGGGCAGACCGGGCGCGCAGTGCCGAAGTGCTGCAGCGCACCAACATggg                                                                           >  2:248128/1‑134 (MQ=255)
                                                            gcgcACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGtt                                                                        <  3:81699/135‑1 (MQ=255)
                                                                aCCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGCAGTGCTGCAGGGCACCAACAAGGTGTtcatc                                                                    >  1:339539/1‑135 (MQ=255)
                                                                      gtccAATGGGGCCCGGGGGGGACCACGGGCCGGGGGGCGGGGGCCATGCGGGAAGGCGCCCGTCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCaccacc                                                              <  6:184094/132‑1 (MQ=255)
                                                                      ctccAATTGGGCCGGGGGGTGACCACGGGCCTGGGTGCCGGGGCAAATCCGGACGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCaccacc                                                              <  5:326377/132‑1 (MQ=255)
                                                                       cgCAACTGGGCCCGGGGGGTGCCCAGGGGCCGGGGGGCCGGGCACAAGCGGGACGCCGGCCGCCCGCCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCaccacc                                                              <  5:86472/134‑1 (MQ=255)
                                                                                ggCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGt                                                     >  3:161388/1‑135 (MQ=255)
                                                                                ggCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGt                                                     >  7:123318/1‑135 (MQ=255)
                                                                                gccccGGGGGTCACCAAGGGCCGGGGGGCGGGGGCCAATCCGCACGGCGGCCGTCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGt                                                     <  1:288186/131‑1 (MQ=255)
                                                                                  gacaggggGTGACGACGGGCCTGGGTGCGGGTGCCAATCCGTAAGTCGGCGGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGTCACCAACATGGTGTTCATCACCACCGGCATGGGTg                                                    <  8:454739/131‑1 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGATCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCTCCGGCATGGGTgg                                                   >  3:314865/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCCCCAACATGGTGGTCATCACCACCGGCATGGGTgg                                                   >  7:201887/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTGATCACCACCGGCATGGGTgg                                                   >  5:341566/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTgg                                                   >  6:184085/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTgg                                                   >  2:59188/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTgg                                                   >  3:403768/1‑135 (MQ=255)
                                                                                  cACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACAAGGTGGTCATCACCACCGGCATGGGTgg                                                   >  4:261618/1‑135 (MQ=255)
                                                                                       ggTGTGACCATGGGCCGGGGTGCCGGTGCAAATCCGGAAGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGcggt                                               <  4:268768/135‑1 (MQ=255)
                                                                                       ggTGTCACCATGGGCCGGGGGGCGGGGGCAAATCCGGACGGCGGCCGCCCGCCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGcggt                                               <  2:427704/135‑1 (MQ=255)
                                                                                         tgtgACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGcggtac                                             >  1:252116/1‑135 (MQ=255)
                                                                                           tgACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCCGCCGTCAGGCCGCGCTGGAAGACCGGGAGCGCATTTGCGAAGGGCTGCAGGGCCCCCACAAGGTGTTGCTAACCACCCGGATGGGTGGcggtaccg                                           >  3:354840/1‑135 (MQ=255)
                                                                                              ccAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGcggtaccggta                                        >  4:140285/1‑135 (MQ=255)
                                                                                                 tgggCCGGGGTGCCGGTGCAATTCGGTAAGGGGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGgtaccggtaccg                                     <  3:409424/134‑1 (MQ=255)
                                                                                                 cgggccggggggccgggCCCACCCCGCACGGCGGCGGCCCGCCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGgtaccggtaccg                                     <  3:392963/134‑1 (MQ=255)
                                                                                                 aGGGCCGGGGTGCGGGGGCCAATCCGGAAGGCGGCCGCCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGgtaccggtaccg                                     <  1:345609/135‑1 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGGCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCCACATGGTGTTCATCACCACCGGCAAGGGTGGCGGTaccgggaccggt                                   >  6:384450/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCCACATGGTGTTCATCACCACCGGCATGGGTGGCGGTaccggtaccggt                                   >  7:442623/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGtt                                                                        <  2:358426/96‑1 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGtt                                                                        >  1:358426/1‑96 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTaccggtaccggt                                   >  8:36235/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTaccggtaccggt                                   >  4:429419/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTaccggtaccggt                                   >  4:172722/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGCCATCGGTGGCggtaccggtaccctgtc                                 >  7:454739/1‑130 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGAACGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTaccggtaccggt                                   >  8:457456/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCCCGCTGGAAGACCGGGAGCGCATTTCAGAAGTGCAGCAGGGCACCAACAAGGTGTTCATAACCACCGGCAGGGGTGGCGGTgccgggaccggt                                   >  1:121571/1‑133 (MQ=255)
                                                                                                     ccTGGGTGCCGGTGCCAATCCGGAAGTCCGCCGTCAGGCCGCGGTGGGAGCCCGGGCGCCCATTTCCGAACTGCTGCAGGGCCCCAACAAGGTGTTCATCACCACCCGCAGTGGTGGCGGTACccgtgcccgtgc                                 >  6:310782/1‑135 (MQ=255)
                                                                                                          gtgccggtgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTTATCGCATTTCCCAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCCCCGGCATGGGTGGCCGTACCGGTACCCGTGCAgcgc                            >  4:349715/1‑135 (MQ=255)
                                                                                                          gtgccggtgccAATCCGGAAGTCGGCCGTCAGGCCGCGCAGGGAGACCCTGAGCGCATTGCCGCACTGCTGCATGGCCCCCACATTGTGTTCATCACCACCCGCAAGGGTGGCGGTTCCGGATCCCGTCCAgcgc                            >  5:46502/1‑135 (MQ=255)
                                                                                                          gtgccggtgccAATCCGGAAGTCGGCCGTCAGGCCGCACTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTTCCGGTCCCGGTGCAgcgc                            >  3:351777/1‑135 (MQ=255)
                                                                                                          gtgccggtgccAATCCGGAAGTCGGCCGTCACGCCCCCCTGGAACCCCGTGAGCGCATTGCCGAAGAGCAGCAGGGCACCACCATAGTGTTCATCACCACCCGCATGGGTGGCCGTTCCCGTACCGGTGCAgcg                             >  2:416354/1‑134 (MQ=255)
                                                                                                           tgccggtgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGCCCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCCCCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGcc                           >  5:373136/1‑135 (MQ=255)
                                                                                                           ggccggggccAAGCGGGAAGGCGGCGGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGcc                           <  4:456639/134‑1 (MQ=255)
                                                                                                            gccggtgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCCCCAACATGGTGTTCATCCCCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCa                          >  3:388433/1‑135 (MQ=255)
                                                                                                                gtgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCCCCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCa                      >  6:215927/1‑135 (MQ=255)
                                                                                                                 tgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCAt                     >  2:323396/1‑135 (MQ=255)
                                                                                                                 tgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCAt                     >  5:423966/1‑135 (MQ=255)
                                                                                                                 tgccAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCAt                     >  7:44321/1‑135 (MQ=255)
                                                                                                                     aaTCGGTAAGGCGGCCGGCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGcc                 <  5:83683/135‑1 (MQ=255)
                                                                                                                     aaTCCGGACGTCGGCCGGCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGcc                 <  6:137128/135‑1 (MQ=255)
                                                                                                                           gacgTCGGCCGTCGGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTg           <  4:455842/132‑1 (MQ=255)
                                                                                                                           gAAGTCGGCCGGCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTg           <  2:366708/135‑1 (MQ=255)
                                                                                                                           gAAGGCGGCGGTCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTg           <  6:127767/135‑1 (MQ=255)
                                                                                                                            aggTCGGCCGTCCGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTgg          <  1:248324/133‑1 (MQ=255)
                                                                                                                            acgtgggCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTgg          <  3:213474/130‑1 (MQ=255)
                                                                                                                              gTCGGCCGTCAGGCCGCGCTGGAAGCCCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCCATCATTGCCGAAGTGGc         >  5:390338/1‑134 (MQ=255)
                                                                                                                                gggccggccggccGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCaa      <  7:193388/134‑1 (MQ=255)
                                                                                                                                cggccggcAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCaa      <  6:430468/135‑1 (MQ=255)
                                                                                                                                cGGCCGTACGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACAAAGATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCaa      <  8:210304/135‑1 (MQ=255)
                                                                                                                                    cGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCAAGGaa  >  3:452961/1‑135 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGGTACCGGTACCGGTGCAGCGCCAATCATTGCCGAAGTGGCCAAGGAA  >  ACB122/1522138‑1522402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: