Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,496 G→T 5.5% intergenic (+92/‑90) 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,4960GT 5.5% 318.3 / 5.0 108intergenic (+92/‑90)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (49/52);  new base T (2/4);  total (52/56)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 109 total aligned reads displayed.

ATAACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGA  >  NC_000913/1530355‑1530645
                                                                                                                                             |                                                                                                                                                     
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)
          gAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                                                                                     >  1:555026/1‑150 (MQ=255)
           aaaaaTGAAGTTTTTGACTTTATTGGTAAATGGTCAAAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACACAAAGACATGAGAGCAAGGCACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCt                                                                                                                                    <  1:870921/150‑1 (MQ=255)
           aaaaaTGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCt                                                                                                                                    <  2:1731937/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)
                           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:738792/150‑1 (MQ=255)
                                             cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                                                                  <  2:495694/150‑1 (MQ=255)
                                             aaaaaTGAACATCAACAGAAAGATAAAAATAAAAAAAACATACAGAAAGAAATGAGAGCAAGGGACGGGGCTGGGCAAGATCTTGTAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCATATTTTAAc                                                                                                  <  2:814486/147‑1 (MQ=255)
                                               tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                                                                >  1:1331040/1‑150 (MQ=255)
                                               tAATGAACAACACAAGAAATATAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                                                                <  2:1331040/150‑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:1448267/150‑1 (MQ=255)
                                                                aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAg                                                                               <  1:1424952/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)
                                                                                    atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGTAATAACTAAAAAGATAAAAAAGCACGCCAGGTg                                                           >  2:1137541/1‑150 (MQ=255)
                                                                                     tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGa                                                          >  1:1517646/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)
                                                                                                   agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTGTAACGCCTACAGGTGATCTACATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGGATTCAAAGTTTg                                            >  1:265122/1‑150 (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)
                                                                                                         aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTTAt                                      >  1:134448/1‑150 (MQ=255)
                                                                                                         aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                      >  1:139749/1‑150 (MQ=255)
                                                                                                         aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                      <  2:108401/150‑1 (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)
                                                                                                                      ggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGa                         >  2:402556/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:490671/1‑150 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                    >  1:1607746/1‑150 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                    >  1:1330248/1‑150 (MQ=255)
                                                                                                                           gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                    <  2:318777/150‑1 (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)
                                                                                                                                             |                                                                                                                                                     
ATAACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGA  >  NC_000913/1530355‑1530645

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: