Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 899,961 T→A 16.5% intergenic (‑38/‑88) ybjP ← / → ybjQ lipoprotein/UPF0145 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273899,9610TA16.5% 940.5 / 162.4 334intergenic (‑38/‑88)ybjP/ybjQlipoprotein/UPF0145 family protein
Reads supporting (aligned to +/- strand):  ref base T (131/148);  new base A (25/30);  total (156/178)
Fisher's exact test for biased strand distribution p-value = 8.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 341 total aligned reads displayed.

CAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTGTGGT  >  CP009273/899812‑900110
                                                                                                                                                     |                                                                                                                                                     
cAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCt                                                                                                                                                       <  1:2655054/150‑1 (MQ=255)
cAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCt                                                                                                                                                       <  1:2144481/150‑1 (MQ=255)
    gCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGa                                                                                                                                                   <  2:3925405/150‑1 (MQ=255)
    gCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGa                                                                                                                                                   >  1:1643453/1‑150 (MQ=255)
      aaGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGata                                                                                                                                                 <  2:4913318/150‑1 (MQ=255)
      aaGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGata                                                                                                                                                 >  2:1236240/1‑150 (MQ=255)
      aaGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGata                                                                                                                                                 >  1:4629787/1‑150 (MQ=255)
      aaGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGata                                                                                                                                                 <  1:4141638/150‑1 (MQ=255)
        ggACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGataat                                                                                                                                               <  2:3626128/150‑1 (MQ=255)
           ccACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATaataat                                                                                                                                            <  1:1530184/150‑1 (MQ=255)
           ccACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATaataat                                                                                                                                            <  1:206997/150‑1 (MQ=255)
               tGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGata                                                                                                                                        <  2:1872021/150‑1 (MQ=255)
               tGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGata                                                                                                                                        <  1:3993642/150‑1 (MQ=255)
                  gtgtGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATcc                                                                                                                                     <  1:4758117/150‑1 (MQ=255)
                    gtgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCtt                                                                                                                                   >  1:1647136/1‑150 (MQ=255)
                    gtgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCtt                                                                                                                                   >  1:1005315/1‑150 (MQ=255)
                      gCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTc                                                                                                                                 <  2:4622927/150‑1 (MQ=255)
                      gCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTc                                                                                                                                 <  2:4622972/150‑1 (MQ=255)
                      gCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTc                                                                                                                                 <  1:1282468/150‑1 (MQ=255)
                      gCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTc                                                                                                                                 <  2:3461090/150‑1 (MQ=255)
                          ttATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCaata                                                                                                                             >  2:760181/1‑150 (MQ=255)
                            aTCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCaataat                                                                                                                           <  1:2924721/150‑1 (MQ=255)
                              cTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAg                                                                                                                         >  1:2397713/1‑150 (MQ=255)
                                ttATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCg                                                                                                                       <  2:3415177/150‑1 (MQ=255)
                                     aGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCa                                                                                                                  <  1:3140759/150‑1 (MQ=255)
                                      gCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAg                                                                                                                 >  2:2257237/1‑150 (MQ=255)
                                          gAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTg                                                                                                             <  1:1568117/150‑1 (MQ=255)
                                              gACTGTGGTTCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTATCGCATAATTATTTCCTTATAAGCTATCGCTCTTAAAGCGTTCTACGATAATCATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         <  1:3129116/150‑1 (MQ=255)
                                              gACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         >  2:92996/1‑150 (MQ=255)
                                              gACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         >  1:221153/1‑150 (MQ=255)
                                              gACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         >  2:3748429/1‑150 (MQ=255)
                                               aCTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAAt                                                                                                        <  2:1326173/150‑1 (MQ=255)
                                                cTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATg                                                                                                       <  1:4391644/150‑1 (MQ=255)
                                                   tggtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCtt                                                                                                    >  1:788771/1‑150 (MQ=255)
                                                   tggtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCtt                                                                                                    >  2:4763831/1‑150 (MQ=255)
                                                    ggtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCttt                                                                                                   >  1:3664121/1‑150 (MQ=255)
                                                         aGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTgaga                                                                                              <  2:2191958/150‑1 (MQ=255)
                                                         aGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTgaga                                                                                              <  2:4600741/150‑1 (MQ=255)
                                                         aGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTgaga                                                                                              >  2:2130672/1‑150 (MQ=255)
                                                             aCTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGa                                                                                          <  1:419346/150‑1 (MQ=255)
                                                               tGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAg                                                                                        <  1:659363/150‑1 (MQ=255)
                                                               tGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAg                                                                                        <  1:1439638/150‑1 (MQ=255)
                                                                 aggagCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGc                                                                                      >  1:644475/1‑150 (MQ=255)
                                                                  ggagCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCt                                                                                     >  1:1959494/1‑150 (MQ=255)
                                                                    agCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTa                                                                                   <  1:4482166/150‑1 (MQ=255)
                                                                    agCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTa                                                                                   <  2:412357/150‑1 (MQ=255)
                                                                     gCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAg                                                                                  <  2:1855372/150‑1 (MQ=255)
                                                                       aGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCa                                                                                >  2:1405595/1‑150 (MQ=255)
                                                                         tGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCaaa                                                                              >  2:1021592/1‑150 (MQ=255)
                                                                         tGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCaaa                                                                              >  1:3448781/1‑150 (MQ=255)
                                                                          gcgcAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAAc                                                                             <  1:2591720/150‑1 (MQ=255)
                                                                          gcgcAGGGGATAAGCATTGTCAATGTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAAc                                                                             <  1:3945381/150‑1 (MQ=255)
                                                                            gcAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACaa                                                                           >  1:5141677/1‑150 (MQ=255)
                                                                            gcAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACaa                                                                           <  2:3159661/150‑1 (MQ=255)
                                                                              aGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAgg                                                                         >  1:583139/1‑150 (MQ=255)
                                                                                 ggATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGgaga                                                                      >  2:2628529/1‑150 (MQ=255)
                                                                                 ggATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAAAAAGgaga                                                                      >  2:562581/1‑150 (MQ=255)
                                                                                  gATAAGCATTGTCTATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgat                                                                     <  1:1074774/150‑1 (MQ=255)
                                                                                  gATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgat                                                                     >  2:2972241/1‑150 (MQ=255)
                                                                                   aTAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatc                                                                    <  2:4988784/150‑1 (MQ=255)
                                                                                   aTAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatc                                                                    >  2:4243306/1‑150 (MQ=255)
                                                                                       gCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcgatc                                                                <  1:3383978/150‑1 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               <  2:2510445/150‑1 (MQ=255)
                                                                                         aTTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                              <  1:4336409/150‑1 (MQ=255)
                                                                                         aTTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                              >  2:2300635/1‑150 (MQ=255)
                                                                                         aTTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTATGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                              <  1:2040379/150‑1 (MQ=255)
                                                                                                aTTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAAttt                                                       <  1:4839187/150‑1 (MQ=255)
                                                                                                  ttGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                                                     <  2:1911130/150‑1 (MQ=255)
                                                                                                      tGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                                                     <  2:3357659/146‑1 (MQ=255)
                                                                                                      tGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                                                     >  1:3357659/1‑146 (MQ=255)
                                                                                                       gTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACaa                                                <  1:4428999/150‑1 (MQ=255)
                                                                                                       gTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACaa                                                <  1:3802405/150‑1 (MQ=255)
                                                                                                         aGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                              >  2:2843347/1‑150 (MQ=255)
                                                                                                           cgcATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAcccc                                            <  1:1146069/150‑1 (MQ=255)
                                                                                                            gcATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCa                                           <  2:2495074/150‑1 (MQ=255)
                                                                                                              aTAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAc                                         <  2:2054872/150‑1 (MQ=255)
                                                                                                              aTAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAc                                         <  1:181710/150‑1 (MQ=255)
                                                                                                                 aTGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActct                                      <  2:2892422/150‑1 (MQ=255)
                                                                                                                  tGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTg                                     >  1:2476060/1‑150 (MQ=255)
                                                                                                                    aTTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGa                                   >  1:1344567/1‑150 (MQ=255)
                                                                                                                     tttCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                                  <  2:4403547/150‑1 (MQ=255)
                                                                                                                     tttCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                                  <  2:2695489/150‑1 (MQ=255)
                                                                                                                        ccTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGc                               >  1:5020468/1‑150 (MQ=255)
                                                                                                                         cTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGcc                              <  2:1163520/150‑1 (MQ=255)
                                                                                                                          ttATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa                             >  2:4011284/1‑150 (MQ=255)
                                                                                                                             taAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGAc                          >  2:507776/1‑150 (MQ=255)
                                                                                                                                   aTCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGt                    >  2:4153081/1‑150 (MQ=255)
                                                                                                                                    tCGCTCTGAAAGCGTTCTACGATAATTATGATATCCTTTCAATATTCGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGtt                   <  1:3869357/150‑1 (MQ=255)
                                                                                                                                      gCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGa                 >  1:2805313/1‑150 (MQ=255)
                                                                                                                                        tctGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAat               >  2:2437627/1‑150 (MQ=255)
                                                                                                                                          tGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTGTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAatat             >  2:1664466/1‑150 (MQ=255)
                                                                                                                                            aaaGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTg           <  1:3748563/150‑1 (MQ=255)
                                                                                                                                                cGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgt       >  2:3121682/1‑150 (MQ=255)
                                                                                                                                                cGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgt       >  2:857057/1‑150 (MQ=255)
                                                                                                                                                 gTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgtg      >  1:3549136/1‑150 (MQ=255)
                                                                                                                                                  ttCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgtgt     >  2:4313429/1‑150 (MQ=255)
                                                                                                                                                  ttCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgtgt     <  2:2154510/150‑1 (MQ=255)
                                                                                                                                                    cTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTgtgg   >  1:1824228/1‑150 (MQ=255)
                                                                                                                                                     tACGATATTAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGATGGCTTAGCAAACAAGGAGATCTATCATGCATTTTTCTTCAACACCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGtgtgtt  <  2:710360/150‑3 (MQ=255)
                                                                                                                                                     tACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTgtggt  >  2:1215596/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTGTGGT  >  CP009273/899812‑900110

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