Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,584 A→T 12.5% P787P (CCA→CCT ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,5840AT12.5% 119.2 / 3.8 69P787P (CCA→CCTftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base A (37/20);  new base T (3/6);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGA  >  NC_000913/935435‑935732
                                                                                                                                                     |                                                                                                                                                    
ccTGGACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCca                                                                                                                                                      <  1:1518297/150‑1 (MQ=255)
  tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                    <  2:991754/150‑1 (MQ=255)
  tGTACAGCAGCCGCAACAACCGGATGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                    <  1:36339/150‑1 (MQ=255)
     aCAGCAGCCGCAACAACCGCTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca                                                                                                                                                 <  2:1071199/150‑1 (MQ=255)
      caccacccGCAACACCCGGATGCACCGCAGCAGCAAGAGCACCAGCCGCACCAACCACGTCCGCCGCAGCCGCAGTCTCAGCAGCCACAACAGCCGCTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag                                                                                                                                                <  1:348754/144‑1 (MQ=21)
       agcagcCGCATGATCCGTTTGCAGGGCAGGAGCATTATCAGCAGCCTCATCAACATCTTCCGCCGCAGCCGCAGTATCAGCAGCCACACCAGCCGGTTGCTCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagcagc                                                                                                                                               <  2:132121/150‑1 (MQ=18)
       agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagtagc                                                                                                                                               <  1:17568/150‑1 (MQ=255)
        gcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACACTATCAGCAGCCGCAACAACCGGATGCGCCCCCgcagca                                                                                                                                              >  1:28658/1‑150 (MQ=255)
          atacGCATCATCCGTTTGCTCTGCTGCAGGATTATCATCAGCATCAACAAATATATCCGCAGCAGGATGATTTTCAGCAGCCACAACAGCCGGTTGCGCCGCATCCACAATTTCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAat                                                                                                                                            <  2:71448/147‑1 (MQ=12)
           gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAGCCAGTTCCGCCGCAGCCGCAGTATCAGCCGCCACCACAGCCCGTTGCGCCGCAGCCACTTGATCAGCGGAAGCCGCCAGCGGGTGCGCCAcagcagggggg                                                                                                                                           >  1:1293009/1‑145 (MQ=17)
           gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCAAGATCAGCAGCAGCCACACCCGGTTGCGCCACACCAGCaaga                                                                                                                                           >  1:1009478/1‑148 (MQ=21)
             cGCAACAACCGGTTGCAGCGCAGCAGCAATATCAGCAGCCGCCACAACCATATCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCTGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATc                                                                                                                                         <  2:1794030/150‑1 (MQ=21)
             cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATc                                                                                                                                         <  2:12771/150‑1 (MQ=255)
                 acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAAGAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATcagca                                                                                                                                     <  1:495552/150‑1 (MQ=255)
                 acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACCGCCGGTTGCGCCGCAGCCCCACTATCCGCAGCCGCAACACCCGGTTGCGCCACAGCCGCAAAATccgcc                                                                                                                                     >  2:2001324/1‑149 (MQ=18)
                            cACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCAACAACAGCCGGTTGCGCCTCATCCCCAATATCAGCAGCAGCAACAGCCGGTTGAGCGAAATCAGCAATAACAGCTGAAGcaacagc                                                                                                                          >  2:2068221/1‑148 (MQ=18)
                            cACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCAACAACAGCCGGTTGCGCCTCATCCACACTATCCGCAGCAGCATAAACCGGATACGCCGCAGCCGGAAGTTCAGCAGCTGCaaccac                                                                                                                          >  2:2154489/1‑150 (MQ=18)
                                gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACCGCCGGTTGCGCCGCATCACCAATATCCGCAGCATCAAAAACCGGTGGCGCCACAGCAGCAAAATCAGAAGCAGccacccccgac                                                                                                                      >  2:1584846/1‑143 (MQ=18)
                                  agcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTg                                                                                                                    >  1:1262645/1‑150 (MQ=255)
                                       caATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGCTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGcgccg                                                                                                               <  1:1870549/150‑1 (MQ=255)
                                         atatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCAATATCAGCAGCAGCAACCACCGGTTGCGCCACAGCAGCAAAAACAGCAGCAGCAACAACCGCATACGCCgca                                                                                                             >  1:1780976/1‑150 (MQ=25)
                                                    ctcaacaaCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATc                                                                                                  <  1:1653699/148‑1 (MQ=255)
                                                     ggctcaacCTCTTCCGCCGGAGCAGCAGTATCATCACCAACCTCAGCTGTATGCGTTGCAGGAACATTATCAGCAGCAGCAACATCCTTTTGCGCCACAGCACCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATca                                                                                                 <  2:222647/146‑1 (MQ=12)
                                                      catcaacCACATCCGCCGCAGCCGCAGTATCATCAGCCACCTCATCAGTTTGCGCCGCAGCCACAATCTCAGCAGCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATcag                                                                                                <  1:1744791/147‑1 (MQ=21)
                                                      caacaaCCTCATCCGCCGCAGCCGCAGGATCAGCAGCCACCACAGCCGTTTGCGACGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATcag                                                                                                <  1:1736093/150‑1 (MQ=21)
                                                         caacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCAATATCAGCAGCATCAACACCCGGTTCAGCCTAAGCGCGACAATCGGCCGAGGCATCAACCGGGGGCGCCGCAGCAGAAGGAGcagaag                                                                                             >  2:1202675/1‑150 (MQ=16)
                                                           acCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCGACAGCTGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGcc                                                                                           >  1:999118/1‑150 (MQ=255)
                                                             cAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACACCCGGTTGCGCCACATCACCAATATCAGCAGCCGCAACAACCGGATGCGCCGCAGCCGCCGGATCAGCAGccgc                                                                                         >  2:1289663/1‑148 (MQ=18)
                                                               gTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAACCCCAGTATCCGCAGCCGCAACCCCAGGTTGCGCCTAAGCACCACAAACAGCCGCAGCAACAACCGCTTACGCCGCAGCCTCAGAagcagcagcagcca                                                                                       >  2:793340/1‑145 (MQ=16)
                                                                 tCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTTGCAGCAGCCacaacc                                                                                     >  1:708729/1‑149 (MQ=255)
                                                                      cgcagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAg                                                                                >  1:706554/1‑150 (MQ=255)
                                                                            cgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCCGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                                                          >  2:706637/1‑150 (MQ=255)
                                                                            cgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCAGCAACAACAGGTTGAGCCACAGCCGGAATATCAGCAGCCGCAACAACCGGATGCGCCGCAGCCGAAGAATCATCAGCCGCCACAGCCAGTTgcgc                                                                          >  2:626991/1‑150 (MQ=18)
                                                                             gcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACTGGCTGTTgtggc                                                                         <  2:1948874/150‑5 (MQ=255)
                                                                                   tCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                                                   <  1:1962788/150‑1 (MQ=255)
                                                                                     agcagcCACAACCGCCGGGTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTCCGCCACAGCCCCAATAACAAAAGGCGCCTCAACCGGTTGCGCCGCAGCAGGATGATCAGCATCCGCACCAGCCTGTGACGCCTCAACAGc                                                                 >  2:888290/1‑150 (MQ=12)
                                                                                     agcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCTCAACAACAGGTTGCGCCTAAGCAGCAAAATCAGCAGACGCAACAACCGTATGCGCCGCAGCCGAAGAATCAGAAGCACCCACAGCCAGTTACGCAACAACAGc                                                                 >  2:1172411/1‑150 (MQ=255)
                                                                                      tgagcCAGAACACCCGGATGCGCTGCTGGCTCAATATCAGCAGCCGCAACAACCTTTTGCGCCTCAGCAGCAATATAAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                                <  2:1022703/148‑1 (MQ=18)
                                                                                      gcagcCACAACAGCCTGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATAACAGCAGCCGCATCAACCGGTTGCGCCGAAGCATCAGAAAAAGCAGCCACAACAGCCAGATGCGCACCAACCGCa                                                                >  2:1451860/1‑150 (MQ=21)
                                                                                      gcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                                >  2:1230771/1‑150 (MQ=255)
                                                                                          ccACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACCGCCAGTTGCGCCACAACCGCCGGAt                                                            >  1:2210206/1‑150 (MQ=255)
                                                                                           cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCCGCAATATCAGCAGCCGCAACAACCGTTTGCGCCGCAGCAGCAGTATCCGCAGCCACAACAGCCAGTTGCGCCCCACCCGCAggctc                                                           >  2:1070265/1‑147 (MQ=18)
                                                                                                 aGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCCGCCACATCAGCCAGTTGCGCCCCAACCGCAGGATACCCTGc                                                     >  2:66405/1‑150 (MQ=255)
                                                                                                 aGCCGCTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                                                     <  1:681865/150‑1 (MQ=255)
                                                                                                  gCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                                                    >  2:603559/1‑150 (MQ=255)
                                                                                                     tgTTGCGCCGCAGCCACAATATCAGCAGACGAAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                                                 <  2:28658/149‑1 (MQ=255)
                                                                                                     ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                                                 >  1:700605/1‑150 (MQ=255)
                                                                                                           gccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCt                                           >  2:666674/1‑150 (MQ=255)
                                                                                                              gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt                                        >  2:1650114/1‑150 (MQ=255)
                                                                                                                aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa                                      <  1:1983583/150‑1 (MQ=255)
                                                                                                                  ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                    <  1:1243920/150‑1 (MQ=255)
                                                                                                                  ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                    >  2:387417/1‑150 (MQ=255)
                                                                                                                      aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCCGGATCCCCTGCTTCATCCGCTGTTGATGCGTa                                >  2:1807063/1‑150 (MQ=255)
                                                                                                                          accgcagcCGCAACAACCTGTTGCTCATCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATgg                            <  2:1480060/147‑1 (MQ=25)
                                                                                                                           cagcagCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACCACAGCCAGTTGCGCCACCACCGCAGGATACACTGCTTCATCCGCGGTTGAGGAGAAATGGc                           >  1:261028/1‑150 (MQ=21)
                                                                                                                            agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                          >  2:1648206/1‑150 (MQ=255)
                                                                                                                            agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                          >  1:1264720/1‑150 (MQ=255)
                                                                                                                              cagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCCGCAGCAGTATCAGCCGCCACCACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCGGCATCCGCTGTTGATGCTTAATGGCGAc                        >  1:1085883/1‑150 (MQ=25)
                                                                                                                              cagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCCTCCGCTGTTGATGCGTAATGGCGAc                        >  1:1092941/1‑150 (MQ=255)
                                                                                                                                gcCGCAACAACCGGTTGCGCCACAGAAGCAATATCCGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTATTGTTGCGTAATGGCGACAg                      >  2:1355029/1‑150 (MQ=255)
                                                                                                                                 actcaacaaCCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGc                     <  1:1575510/147‑1 (MQ=255)
                                                                                                                                  cGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGcc                    >  1:941633/1‑150 (MQ=255)
                                                                                                                                   gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccg                   >  1:1813294/1‑150 (MQ=255)
                                                                                                                                    caacaaCCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgt                  <  1:2068221/150‑1 (MQ=255)
                                                                                                                                       aaacCTGTAGCGCCCGCGCAGCTATATCCCCATCCGCATCCACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg               <  1:1550639/149‑1 (MQ=37)
                                                                                                                                          ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCCGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTg            >  2:987855/1‑150 (MQ=255)
                                                                                                                                          ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTg            <  2:1262645/150‑1 (MQ=255)
                                                                                                                                            ggTTGCGCCTCATCAGAAATATCAGCAGCCGCATCAACCTGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa          <  2:1840938/150‑1 (MQ=37)
                                                                                                                                            ggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACTGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa          >  2:585402/1‑150 (MQ=255)
                                                                                                                                             gTTGCTCCTGAGCATCAATTTCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCAt         <  2:576158/150‑1 (MQ=37)
                                                                                                                                                 cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc     >  1:725674/1‑150 (MQ=255)
                                                                                                                                                   ccACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATCCCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcg   >  2:207486/1‑150 (MQ=255)
                                                                                                                                                    cacaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcgc  >  1:1237143/1‑149 (MQ=255)
                                                                                                                                                    cacaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGATGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcga  >  1:480405/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGA  >  NC_000913/935435‑935732

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: