Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,505 G→A 5.4% intergenic (+101/‑81) 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,5050GA 5.4% 335.8 / 5.8 112intergenic (+101/‑81)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (54/51);  new base A (2/4);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 113 total aligned reads displayed.

AACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTTAAA  >  NC_000913/1530357‑1530654
                                                                                                                                                    |                                                                                                                                                     
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)
     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)
                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:1640135/1‑150 (MQ=255)
                        gACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTa                                                                                                                              >  1:1190565/1‑150 (MQ=255)
                         aCTTGATTGGCAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTAc                                                                                                                             <  2:501856/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                                                                                                                  >  1:1717674/1‑150 (MQ=255)
                                    aaaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                                                                                  <  2:129447/150‑1 (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:1704559/150‑1 (MQ=255)
                                           cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                                                                           <  2:495694/150‑1 (MQ=255)
                                           cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                                                                           <  2:738792/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)
                                                                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)
                                                                                   tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGa                                                                   >  1:1517646/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)
                                                                                                 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)
                                                                                                       aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                               <  2:108401/150‑1 (MQ=255)
                                                                                                            cGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTatat                                          >  1:779809/1‑150 (MQ=255)
                                                                                                              gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAttt                                        >  2:390634/1‑150 (MQ=255)
                                                                                                                   gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                                   >  2:1518749/1‑150 (MQ=255)
                                                                                                                    ggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGa                                  >  2:402556/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:490671/1‑150 (MQ=255)
                                                                                                                         gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                             >  1:1330248/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)
                                                                                                                                        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)
                                                                                                                                                 cACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGtt     >  1:1436309/1‑150 (MQ=255)
                                                                                                                                                   cGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTTaa   >  2:1139434/1‑150 (MQ=255)
                                                                                                                                                    gATTATCATCTTTGCTTTAACGCCTACAGTGTATCAAAATTTAACCTTTGGAAAAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTTaaa  <  2:289107/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACAATTTGAAAAATGAAGTTTTTGACTGGATTGGTAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTTAAA  >  NC_000913/1530357‑1530654

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: