Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,837,041 A→G 5.9% Y126Y (TAT→TAC EO946_RS14635 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349432,837,0410AG 5.9% 426.0 / 11.0 135Y126Y (TAT→TACEO946_RS14635diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (62/64);  new base G (4/4);  total (66/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 135 total aligned reads displayed.

ATCTTGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAGCCGTCCT  >  NZ_CP034943/2836903‑2837174
                                                                                                                                          |                                                                                                                                     
aTCTTGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGAt                                                                                                                                      >  1:1112637/1‑140 (MQ=255)
aTCTTGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGAt                                                                                                                                      >  2:2115021/1‑140 (MQ=255)
aTCTTGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTAATTTTGTATATCTTTTTTTGAt                                                                                                                                      >  1:498619/1‑140 (MQ=255)
   ttGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGa                                                                                                                                   <  1:2309813/140‑1 (MQ=255)
     ggTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATc                                                                                                                                 >  1:856873/1‑140 (MQ=255)
     ggTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATc                                                                                                                                 >  1:2296428/1‑140 (MQ=255)
     ggTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATc                                                                                                                                 >  1:20937/1‑140 (MQ=255)
      gTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCg                                                                                                                                >  2:1507533/1‑140 (MQ=255)
      gTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCg                                                                                                                                >  2:1139941/1‑140 (MQ=255)
        cGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtt                                                                                                                              <  1:362044/140‑1 (MQ=255)
          tGCGAAATAAACAATGAACACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  1:994602/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTTTCTGATCAATCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  2:994757/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGACTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  1:455487/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  1:2033397/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  1:2274687/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  2:937840/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGtttt                                                                                                                            >  1:675992/1‑140 (MQ=255)
          tGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGATTTCTAAAAGCTGCTTCAGCACCTCTGTATCATCAAGCAGTGAATGGACCAGCTCTGTCATATTGTATATCTTTTTTTGATATAACGattt                                                                                                                            >  2:1249201/1‑140 (MQ=255)
                    ataaagAAAACTCAAAAATATAAAATGTGATTGTCAAGCTTTCTAAAAGATGCTTCAGCACATCGGTCAGATCAAGCAGTGAATGGAACAGCTCTGTCATGTTCTATATCTTTTTTTGGTAGATCGTTTTGTTATAGAGt                                                                                                                  >  1:915189/6‑140 (MQ=255)
                    ataaagAAAACTCAAAAAGAGAAAATGTGCTTGTCAAGCTTTCTAAAAGCTGCTTCAGCACATCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGGTAGATCGTTTTGTTATAGAGt                                                                                                                  >  1:311424/6‑140 (MQ=255)
                    ataaagAAAACTCAAAAAGAGAAAATGTGCTTGTCAAGCTTTCTAAAAGCTGCTTCAGCACAGCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCAGTTTGTATATCTTTTTTTGGTCGATCGTTTTGTTATAGAGt                                                                                                                  >  1:1768081/6‑140 (MQ=255)
                    aCAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGt                                                                                                                  <  2:1704079/140‑1 (MQ=255)
                        tGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCttt                                                                                                              >  1:1770973/1‑140 (MQ=255)
                        tGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCttt                                                                                                              >  1:1301004/1‑140 (MQ=255)
                         caaaaCTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTc                                                                                                             <  1:782252/139‑1 (MQ=255)
                                      gagaAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGAc                                                                                                <  2:1368090/140‑1 (MQ=255)
                                         aaaaCGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTgg                                                                                             >  1:1077271/1‑140 (MQ=255)
                                         aaaaCGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATAGAGAGTCTTTCCTGTTTCATTGACTgg                                                                                             >  1:1828992/1‑140 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATag                                                                                                                     >  2:371787/1‑109 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATag                                                                                                                     <  1:371787/109‑1 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGAc                                                                                                <  1:1056796/130‑1 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGAc                                                                                                >  2:1056796/1‑130 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTaaaa                                                                                      >  2:1774178/1‑140 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTaaaa                                                                                      >  2:1908066/1‑140 (MQ=255)
                                                gCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTaaaa                                                                                      >  2:552098/1‑140 (MQ=255)
                                                          cTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAAc                                                                            <  2:74625/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACATCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:994602/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:2033397/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:1074773/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:2274687/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:1618067/140‑1 (MQ=255)
                                                            ttCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCa                                                                          <  2:1413351/140‑1 (MQ=255)
                                                               tAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGt                                                                       <  1:1104316/140‑1 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  >  1:1787181/1‑140 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  <  2:1301004/140‑1 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  <  2:1770973/140‑1 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  <  1:550380/140‑1 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  >  1:2278576/1‑140 (MQ=255)
                                                                    gctgctTCAGCACCTCTGTCGGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGCTAGATCGTTTTGATATCGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGaa                                                                  <  1:759034/140‑1 (MQ=255)
                                                                      tgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAg                                                                <  1:1139941/140‑1 (MQ=255)
                                                                      tgctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAg                                                                <  1:1507533/140‑1 (MQ=255)
                                                                       gctTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGc                                                               >  2:254110/1‑140 (MQ=255)
                                                                               aCCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAgg                                                       <  2:455487/140‑1 (MQ=255)
                                                                                 ctctGTCTTATCAAGCTGTGAATGGAACAGCTCTGTCATTTTTTATATCTTTTTTGGATATATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:1378155/140‑1 (MQ=255)
                                                                                 ctctGTCTGATCAAGCAGTGAATGGACCAGCTCTGTCATTGTGTCTATCTTTTTTTGATAGATCGTTTTGATATAGCGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:245412/140‑1 (MQ=255)
                                                                                 ctctGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:1787181/140‑1 (MQ=255)
                                                                                 ctctGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:2175524/140‑1 (MQ=255)
                                                                                 ctctGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:1112637/140‑1 (MQ=255)
                                                                                 ctctGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGaa                                                     <  2:275870/140‑1 (MQ=255)
                                                                                   ctgtctgATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAAtc                                                   >  1:1193967/1‑140 (MQ=255)
                                                                                   ctgtctgATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAAtc                                                   >  1:504765/1‑140 (MQ=255)
                                                                                      tctgATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAAtctct                                                >  1:684967/1‑140 (MQ=255)
                                                                                       ctgATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAAtctctc                                               <  1:1061194/140‑1 (MQ=255)
                                                                                          aTCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAc                                            >  2:952481/1‑140 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           <  1:587945/140‑1 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           >  1:410798/1‑140 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           <  1:732228/140‑1 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           >  1:2410485/1‑140 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           >  2:732228/1‑140 (MQ=255)
                                                                                           tCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           >  2:1700283/1‑140 (MQ=255)
                                                                                           acAACCAGTGAATGGAACAGCTCTGACATTTTGTAGATCCTTTTTCGATAGATCGTATTGATATCGCGTCTTTCCGGTTTCATTGACTGGTTTACAAATGAAGCAACCACGTGCGCAAGCCCTTCAGGAATCTCTCCAcc                                           <  1:1700283/139‑1 (MQ=255)
                                                                                            cAAGGAGTGAATTGAACAGCTGTGTCATTTTGTATATATTTTTTTGGTAGATCGTTTTGTTATAGAGTCTTTCCTGTTTCATCGACTGGTTTAAAAATGAAGCAAGCACGTCCGGAAGCCCTTCAGGAATCTCTTCAtcc                                          <  2:915189/140‑1 (MQ=255)
                                                                                            cAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGGTAGATCGTTTTGTTATAGAGTCTTTCCTGTTTCATCGACTGGTTTAAAAATGAAGCAACCACGTCCGGAAGCCCTTCAGGAATCTCTTCAtcc                                          <  2:855554/140‑1 (MQ=255)
                                                                                             aaGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTTGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcccg                                         <  2:335789/140‑2 (MQ=255)
                                                                                             aaGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCt                                         <  2:1317657/140‑1 (MQ=255)
                                                                                              aGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGATTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTg                                        >  1:2374421/1‑140 (MQ=255)
                                                                                              aGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGa                                                      >  2:2410171/1‑126 (MQ=255)
                                                                                                     aaTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGaa                                                                                 <  1:1205937/92‑1 (MQ=255)
                                                                                                     aaTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGaa                                                                                 >  2:1205937/1‑92 (MQ=255)
                                                                                                     aaTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           >  1:605847/1‑130 (MQ=255)
                                                                                                     aaTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCAcc                                           <  2:605847/130‑1 (MQ=255)
                                                                                                      aTGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATg                                <  1:737214/140‑1 (MQ=255)
                                                                                                        ggAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTTACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCGCTTCTTTAATCTATCCACCGTGTCTTCTTGTa                              <  2:1517381/140‑1 (MQ=255)
                                                                                                        ggAACAGCTCTGTAATTTTTTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCATGTTTCATTGACTGGTTTAAAAATGAAGCAACCACTTCAGAAAGCCCTTCTTTAATCTATCCACCCTGTCTTCTTGTa                              <  2:1002782/140‑1 (MQ=255)
                                                                                                          aaCAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAg                            >  2:739596/1‑140 (MQ=255)
                                                                                                               ctctGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGaa                       >  1:805197/1‑140 (MQ=255)
                                                                                                                   gTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcg                   <  2:1595902/140‑1 (MQ=255)
                                                                                                                   gTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcg                   <  2:2067416/140‑1 (MQ=255)
                                                                                                                    tCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcgc                  >  2:1289814/1‑140 (MQ=255)
                                                                                                                    tCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcgc                  >  2:775137/1‑140 (MQ=255)
                                                                                                                    tCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcgc                  >  1:1243591/1‑140 (MQ=255)
                                                                                                                    tCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcgc                  >  2:484729/1‑140 (MQ=255)
                                                                                                                    tCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAgcgc                  >  2:1897898/1‑140 (MQ=255)
                                                                                                                         ttGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGcc             >  1:417746/1‑140 (MQ=255)
                                                                                                                             atatCTTTTTTTTATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGATTAAAAATGAAGCAACCACGTCCGAACGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAg         <  2:364388/140‑1 (MQ=255)
                                                                                                                             atatCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAg         <  1:687784/140‑1 (MQ=255)
                                                                                                                             atatCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAg         <  2:20937/140‑1 (MQ=255)
                                                                                                                             atatCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAg         <  2:2296428/140‑1 (MQ=255)
                                                                                                                                  tttttttGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAGCCGTc    <  1:937840/140‑1 (MQ=255)
                                                                                                                                    tttttGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCCCGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAGCCGTCCt  <  2:2277174/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
ATCTTGGTCGTGCGAAATAAACAATGAAAACTCAATAGGAGAAAACGTGCGGGTCAAGCTTTCTAAAAGCTGCTTCAGCACCTCTGTCTGATCAAGCAGTGAATGGAACAGCTCTGTCATTTTGTATATCTTTTTTTGATAGATCGTTTTGATATAGAGTCTTTCCTGTTTCATTGACTGGTTTAAAAATGAAGCAACCACGTCCGAAAGCCCTTCAGGAATCTCTCCACCCTGTCTGCATGTAAGATGAAAGCGCCTGCCGGAGCCGTCCT  >  NZ_CP034943/2836903‑2837174

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