Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,499 A→C 5.5% intergenic (+95/‑87) ydcD → / → yncI putative immunity protein for RhsE/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,4990AC 5.5% 334.2 / 7.3 110intergenic (+95/‑87)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base A (52/52);  new base C (2/4);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 110 total aligned reads displayed.

ATAACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAA  >  NC_000913/1530355‑1530646
                                                                                                                                                |                                                                                                                                                   
aTAACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCATAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAc                                                                                                                                                <  2:2900/150‑1 (MQ=255)
  aaCAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                                                                              >  1:1733452/1‑150 (MQ=255)
   aCAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGATAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAt                                                                                                                                             <  2:174681/150‑1 (MQ=255)
    cAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                                                                            >  1:108401/1‑150 (MQ=255)
    aaaTTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAAAACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAtt                                                                                                                                            <  2:860139/149‑1 (MQ=255)
       tttGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                                                                                         <  1:517249/150‑1 (MQ=255)
       tttGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCACAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAGATTGTGATGTAAATCACGATTTTc                                                                                                                                         >  2:11727/1‑150 (MQ=255)
       tttGAAAAATGAAGTTTTTGACTGGATTGGAAATTGGTCAAAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTc                                                                                                                                         <  2:341938/150‑1 (MQ=255)
         tcaaaaaTGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                                                                       <  1:712491/148‑1 (MQ=255)
           aaaaaTGAAGTTTTTGACTTTATTGGTAAATGGTCAAAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACACAAAGACATGAGAGCAAGGCACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCt                                                                                                                                     <  1:870921/150‑1 (MQ=255)
           aaaaaTGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCt                                                                                                                                     >  2:1273167/1‑150 (MQ=255)
               aTGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAATAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGc                                                                                                                                 <  2:1729363/150‑1 (MQ=255)
                  aaGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATTAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGATAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                                                                              >  2:159945/1‑150 (MQ=255)
                   aGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                                                                             <  2:740905/150‑1 (MQ=255)
                     tttttGACTGGATTGGTAAATGGTCATAATTAACAACACCAGAAAGAGAAAAAAATAAAATACATACACAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                                                                                           <  1:880975/150‑1 (MQ=255)
                      ttttGACTGCATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATTTAAAATACATACAGAAAGACATGAGAGAAAGGTACGGGGCAGGGCAAGATCATCAAAATTGTGATCTAAATCACGATTTTCATCTTTGCTTTAACg                                                                                                                          >  2:887720/1‑150 (MQ=255)
                       tttGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGc                                                                                                                         <  2:545752/150‑1 (MQ=255)
                        ttGACTGGATTGGTAAATGGTCACAATGAACAACACCAGAAAGAGAAAAATATAAAATACATTCAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAGATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                                                                                        <  1:11727/150‑1 (MQ=255)
                         tGACTGGATTCGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAACACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCt                                                                                                                       <  2:965086/150‑1 (MQ=255)
                          gACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTa                                                                                                                      >  1:1190565/1‑150 (MQ=255)
                          gACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTa                                                                                                                      >  1:1640135/1‑150 (MQ=255)
                           aCTTGATTGGCAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAc                                                                                                                     <  2:501856/150‑1 (MQ=255)
                           aCTGGATTGGTAAAAGGTCATAATCAACAACACCAGAAAGAGAAAAATATAAAAAACATACAGAAAGACATGAGAGCAAGGGCCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAc                                                                                                                     <  1:140025/150‑1 (MQ=255)
                             tGGATTGGTAAATGGTCATAATGAACAACACCATAAAGATAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                                                                   <  2:1177591/150‑1 (MQ=255)
                              ggATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGTAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                                                                  >  2:314187/1‑150 (MQ=255)
                              ggATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACAAACAGAAAGAAATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                                                                  <  2:293331/150‑1 (MQ=255)
                               gATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATTAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                                                                 <  1:1215471/150‑1 (MQ=255)
                               gATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                                                                 <  2:957142/150‑1 (MQ=255)
                                aTTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAAAATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTg                                                                                                                >  1:1313285/1‑150 (MQ=255)
                                 ttGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATATTCATCTATGCTTTAACGCCTACAGgaga                                                                                                               >  1:1694060/1‑150 (MQ=255)
                                   ggTAAATGGTCATAATGAACAACACCAGAAAGAGAAAATTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACTATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                                                                             >  1:6451/1‑150 (MQ=255)
                                      aaaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                                                                          <  2:129447/150‑1 (MQ=255)
                                      aaaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                                                                          >  1:1717674/1‑150 (MQ=255)
                                       aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaaa                                                                                                         <  1:391847/150‑1 (MQ=255)
                                           gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGCCGGGGCTGGGCAAGATCATGAAAATTGTGATGTAAATCAAGATTTTCATCTTTGCGTTAACGTCGACAGGTGATCAAAATTTa                                                                                                     >  2:169580/1‑150 (MQ=255)
                                           gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                                                                                     >  2:1014294/1‑150 (MQ=255)
                                             cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                                                                   <  2:495694/150‑1 (MQ=255)
                                             cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                                                                   <  2:1704559/150‑1 (MQ=255)
                                             aaaaaTGAACATCAACAGAAAGATAAAAATAAAAAAAACATACAGAAAGAAATGAGAGCAAGGGACGGGGCTGGGCAAGATCTTGTAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCATATTTTAAc                                                                                                   <  2:814486/147‑1 (MQ=255)
                                                aaTGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCCTGAAAATTGTGATGTAAATCACGATTTACAACTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                                                                >  1:279874/1‑150 (MQ=255)
                                                aaTGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTATCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                                                                >  1:314166/1‑150 (MQ=255)
                                                 aTGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGCGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAGCGTTGCTTTAACGCCTACAGGTGATCAAAATTTAACAttt                                                                                               >  1:1039730/1‑150 (MQ=255)
                                                          acCAGAAAGAGAAAAATATAAAATACATACATAAAGACATGAGAGCAAGTGACTGGGCAGGGAAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                                                                                      <  1:264401/150‑1 (MQ=255)
                                                           ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTa                                                                                     >  1:180845/1‑150 (MQ=255)
                                                            cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaa                                                                                    >  2:589649/1‑150 (MQ=255)
                                                              gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                                                                                  >  2:1224988/1‑150 (MQ=255)
                                                               aaaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAAGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                                                                 <  1:1152131/150‑1 (MQ=255)
                                                                aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                                                                                <  1:1424952/150‑1 (MQ=255)
                                                                aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                                                                                <  1:1448267/150‑1 (MQ=255)
                                                                aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                                                                                >  1:1371971/1‑150 (MQ=255)
                                                                 agagAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                                                                               <  2:1485708/150‑1 (MQ=255)
                                                                 agagAAAAATATAAAATACAAACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGa                                                                               <  2:75744/150‑1 (MQ=255)
                                                                  tagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGAt                                                                              <  1:1491185/149‑1 (MQ=255)
                                                                    gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                                                                            >  2:1252336/1‑150 (MQ=255)
                                                                      aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGAGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                                                                          >  2:1470402/1‑150 (MQ=255)
                                                                       aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaa                                                                         <  1:380231/150‑1 (MQ=255)
                                                                         atatAAAATACATACAGAAAGACATGAGAGCAAGTGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                                                                       <  1:134676/150‑1 (MQ=255)
                                                                         atatAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                                                                       <  2:1523854/150‑1 (MQ=255)
                                                                           ataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGa                                                                     <  1:1601639/150‑1 (MQ=255)
                                                                            taAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGAc                                                                    >  1:552042/1‑150 (MQ=255)
                                                                            aaaaaTACATACAGAAAGACATTAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                                                                    <  2:713823/149‑1 (MQ=255)
                                                                                    atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGTAATAACTAAAAAGATAAAAAAGCACGCCAGGTg                                                            >  2:1137541/1‑150 (MQ=255)
                                                                                    atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg                                                            >  2:1256382/1‑150 (MQ=255)
                                                                                      acaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                                                          >  2:389495/1‑150 (MQ=255)
                                                                                            aGACATTAGAGAAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGTAATAACTAAAAAGATAATAAAGGACGCCAGGTGAGTATTCa                                                    <  2:1119490/150‑1 (MQ=255)
                                                                                            aGACATGAGAGCAAGGGACGGGGCTGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTATAACGCCTACAGGTGAACAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGACAGTATTCa                                                    >  2:753240/1‑150 (MQ=255)
                                                                                                  gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGttt                                              <  2:1385314/150‑1 (MQ=255)
                                                                                                   agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTg                                             <  1:1156689/150‑1 (MQ=255)
                                                                                                    gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTTCTTTAACGCCTACAGGTTATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGc                                            <  2:1263191/150‑1 (MQ=255)
                                                                                                     agCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCt                                           <  2:695870/150‑1 (MQ=255)
                                                                                                      gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAATGTTTGCtt                                          >  1:289107/1‑150 (MQ=255)
                                                                                                       cAAGGGTCGGGGCAGGGCAAGATCATGAAAATTTTTATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGTTGATCAAAATTTAACCTTTGTAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTg                                         <  1:475027/150‑1 (MQ=255)
                                                                                                       cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTg                                         <  2:1733452/150‑1 (MQ=255)
                                                                                                        aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGa                                        >  1:794650/1‑150 (MQ=255)
                                                                                                         aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTTAt                                       >  1:134448/1‑150 (MQ=255)
                                                                                                         aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                       >  1:139749/1‑150 (MQ=255)
                                                                                                              cGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTatat                                  >  1:779809/1‑150 (MQ=255)
                                                                                                                gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAttt                                >  2:390634/1‑150 (MQ=255)
                                                                                                                   cagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAg                             <  2:555026/150‑1 (MQ=255)
                                                                                                                     gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                           >  2:1518749/1‑150 (MQ=255)
                                                                                                                       gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAc                         >  1:1370996/1‑150 (MQ=255)
                                                                                                                       gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCGTTGCTTTAACGCCGACAGGTGATCAAAATTTAACCTTTGGAATAACAAAAAAGATAAAAACGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAc                         >  1:380682/1‑150 (MQ=255)
                                                                                                                         aaGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGCAATAACTAAAACGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                       <  2:683736/150‑1 (MQ=255)
                                                                                                                          caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcccc                      >  2:926053/6‑150 (MQ=39)
                                                                                                                          caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCAAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcccc                      >  2:824757/6‑150 (MQ=39)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                     >  1:1330248/1‑150 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                     >  1:1607746/1‑150 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                     <  2:318777/150‑1 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                     >  1:490671/1‑150 (MQ=255)
                                                                                                                              cATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGat                  >  2:578806/1‑150 (MQ=255)
                                                                                                                              aaTGAAAATTGTTATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAATATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGat                  <  1:609901/149‑1 (MQ=255)
                                                                                                                                tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGatat                <  2:1717674/150‑1 (MQ=255)
                                                                                                                                  aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGCCCCCTGATATAc              >  1:603830/1‑150 (MQ=255)
                                                                                                                                   aaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGTCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACg             <  2:199997/150‑1 (MQ=255)
                                                                                                                                     aTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGac           <  2:486618/150‑1 (MQ=255)
                                                                                                                                        gttaTGTAAATAACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAATAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaag        <  1:852927/147‑2 (MQ=255)
                                                                                                                                          gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGAcaacaa      >  1:1494819/1‑150 (MQ=255)
                                                                                                                                            tGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCATGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAgg    <  1:24963/150‑1 (MQ=255)
                                                                                                                                             gTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGa   >  1:1575110/1‑150 (MQ=255)
                                                                                                                                              tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGaa  >  1:1163313/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATAACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAA  >  NC_000913/1530355‑1530646

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: