Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,607:1 +G 7.0% coding (443/822 nt) yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6071.G 7.0% 455.5 / 28.8 114coding (443/822 nt)yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (50/56);  new base G (3/5);  total (53/61)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 115 total aligned reads displayed.

CCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATC  >  NC_000913/2449459‑2449754
                                                                                                                                                     |                                                                                                                                                   
cccATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACt                                                                                                                                                    >  2:710172/1‑150 (MQ=255)
 ccATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTa                                                                                                                                                   >  2:1550585/1‑150 (MQ=255)
   aTCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATc                                                                                                                                                 >  1:1360742/1‑150 (MQ=255)
    tCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGTGTGAGCAGTTTTGACTTACTTACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCa                                                                                                                                                <  1:467838/150‑1 (MQ=255)
    tCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCa                                                                                                                                                <  1:610135/150‑1 (MQ=255)
     cTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAg                                                                                                                                               >  2:1158179/1‑150 (MQ=255)
      tttCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAgg                                                                                                                                              >  2:930700/1‑150 (MQ=255)
      tttCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAgg                                                                                                                                              >  1:1686894/1‑150 (MQ=255)
        tcatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGAcc                                                                                                                                            <  1:996467/150‑1 (MQ=255)
        tcatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGGCGTGACTATCAGGcc                                                                                                                                            <  2:1243839/150‑1 (MQ=255)
         catcaAATACGTTGAAAGCGTTTTCTCAGCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACTTAATTCTCTTTCCATGAGAGGTCAGGAGGGTGATCGGAGACGGATGGACATATCGTGTCGTGGGTAGCAGGccc                                                                                                                                           >  2:1598415/1‑150 (MQ=255)
          atcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCATGGTGAGCAGTTTTGCCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTAGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcccc                                                                                                                                          >  2:255224/1‑150 (MQ=255)
           tcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCt                                                                                                                                         >  1:540851/1‑150 (MQ=255)
            caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                                                        >  1:1300806/1‑150 (MQ=255)
            caAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                                                        <  2:1385275/150‑1 (MQ=255)
                tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATAc                                                                                                                                    <  1:133829/150‑1 (MQ=255)
                tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATAc                                                                                                                                    <  1:1536635/150‑1 (MQ=255)
                   gattaaaGCGTTTTCTCAACCACACCGCCCGTATACTGCATGGGGAGCAGTTTTGACTTACTGACATAATTGTCATTCCATGAAACGTCATAATGTTGAACATAAACGGATGGACAACTCGTGTCGTGACGTATCAGGGGCCTGATagtc                                                                                                                                  <  1:1072645/146‑4 (MQ=255)
                   gTTCAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATagtc                                                                                                                                  <  1:887442/150‑4 (MQ=255)
                    ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGGGCCTGATAGTCAc                                                                                                                               <  1:1207757/150‑1 (MQ=255)
                     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                                                               <  1:1376728/150‑1 (MQ=255)
                     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGACGTGACTATAAGGCCCCTGATACGTAAc                                                                                                                               >  2:334563/1‑150 (MQ=255)
                     tGAAAGCGTTTTCTCAACCACACCGCCCGGCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                                                               <  1:1370043/150‑1 (MQ=255)
                      gAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACg                                                                                                                              >  1:1170038/1‑150 (MQ=255)
                        aaGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGcc                                                                                                                            >  1:1581997/1‑150 (MQ=255)
                         aGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                                                           >  2:1402592/1‑150 (MQ=255)
                           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacacg                                                                                                                         <  1:1170852/150‑6 (MQ=255)
                           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                                                         >  1:1095889/1‑150 (MQ=255)
                             ttttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTg                                                                                                                       <  1:188447/150‑1 (MQ=255)
                                  tcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                                                  >  1:880274/1‑150 (MQ=255)
                                    aaCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                                                <  1:1642446/150‑1 (MQ=255)
                                     aCCACACCGCCCGGCGACTGCAGGGTGATCAGTTTTGAATTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                                               <  2:683340/150‑1 (MQ=255)
                                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                                                          >  1:705577/1‑150 (MQ=255)
                                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACATCACGCCCCCTGGAGAGGCTTCaga                                                                                                          <  1:38697/150‑1 (MQ=255)
                                           ccGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                         >  1:175012/1‑150 (MQ=255)
                                           ccGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                         <  2:1624841/150‑1 (MQ=255)
                                            cGCCCGGCGAATTCAGGCTGAGCAGTTTTGACTTAATGACATAATTCGCATTCCATGAAACGTCATAAGGTTGAACCGACACGGATGGACAAATGGGGTCGTGACTATCAGGCCCATGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                        <  2:289366/150‑1 (MQ=255)
                                               ccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGCACCCTGCAACGTCACGCCCCCGGGAGAGGCTTCAGAGtcct                                                                                                     >  1:461199/1‑148 (MQ=255)
                                                 gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                   <  1:1648071/150‑1 (MQ=255)
                                                   cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                 >  1:1383415/1‑150 (MQ=255)
                                                    gACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCata                                                                                                <  1:448792/150‑1 (MQ=255)
                                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATTAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                              <  1:931323/150‑1 (MQ=255)
                                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                              <  1:975799/150‑1 (MQ=255)
                                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                              <  2:185021/150‑1 (MQ=255)
                                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                              >  1:359173/1‑150 (MQ=255)
                                                          aGGGTGAGCAGTTTTGACTTACTGCCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGACCCCTGGAGAGGATTCAGAGTCATTCATAATaacc                                                                                          >  1:268842/1‑150 (MQ=255)
                                                          aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                                                          <  2:880274/150‑1 (MQ=255)
                                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATTGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                       <  1:320628/150‑1 (MQ=255)
                                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                       <  2:673052/150‑1 (MQ=255)
                                                              tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                                                      >  1:153237/1‑150 (MQ=255)
                                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACATATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                   <  2:546108/150‑1 (MQ=255)
                                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                   >  2:465322/1‑150 (MQ=255)
                                                                        tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGa                                                                            >  2:1023422/1‑150 (MQ=255)
                                                                         gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                                                           >  2:1636805/1‑150 (MQ=255)
                                                                          aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCACCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaa                                                                          <  1:220887/150‑1 (MQ=255)
                                                                            ttACTGACCTAGTTGTCATGCCATGAACAGTGATCAGGTTGATCCGAACCGGCGGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCTTTGATATTAACCTACCATGGACAAGAAAAg                                                                        <  2:131530/150‑1 (MQ=255)
                                                                            ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                                                        <  1:1314856/150‑1 (MQ=255)
                                                                              aCTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                                                                      <  1:1412246/150‑1 (MQ=255)
                                                                                 gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                   >  1:460674/1‑150 (MQ=255)
                                                                                 gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                   >  1:1661135/1‑150 (MQ=255)
                                                                                     tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                                               <  2:603266/150‑1 (MQ=255)
                                                                                      aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                              >  1:749594/1‑150 (MQ=255)
                                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                        >  1:550415/1‑150 (MQ=255)
                                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCAATCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                        >  1:131213/1‑150 (MQ=255)
                                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                        <  2:1242356/150‑1 (MQ=255)
                                                                                             aTTCCATGCAACGTCATAAGGTTGAACAGGTGCGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCGCGGTGGAGAGGCTTCTGAGTCATTGATTATAACCAGCCTGGGCCTAGATATGTTAGGTTCAGTCGCTgg                                                       >  2:279128/1‑149 (MQ=255)
                                                                                             aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                                                       >  1:366820/1‑150 (MQ=255)
                                                                                              ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                      <  1:930700/150‑1 (MQ=255)
                                                                                               tCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                                                     <  1:373393/150‑1 (MQ=255)
                                                                                                ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGTCCACTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCata                                                    >  2:174697/1‑148 (MQ=255)
                                                                                                ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                                                    >  1:859563/1‑150 (MQ=255)
                                                                                                 cctGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                   <  2:1435488/148‑1 (MQ=255)
                                                                                                 cctGAAACGTCATAACGTGCCACCTACACGGACGTACAACTCGGGGCGTGACTATCAGGCCCCGGATACGTCACGCCCCCGGTAGAGGCTGCAGAGTCATTCATAATAACCAACCATGGACAAGAAATGTAAGATTCAATCTCTGCagcg                                                   <  2:120733/148‑3 (MQ=255)
                                                                                                 cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                   <  1:477923/150‑1 (MQ=255)
                                                                                                  aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga                                                  >  2:928517/1‑150 (MQ=255)
                                                                                                   tGAAACGTCATAAGGTTTAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                 <  1:1072671/150‑1 (MQ=255)
                                                                                                   tGAAACGTCATAAGGATGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCATCACAGTCATTCATAATAACCAACCAAGGACAAGAAAAATAATATACAATCTCTGCAGacac                                                 >  1:351443/1‑150 (MQ=255)
                                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                           <  2:830136/150‑1 (MQ=255)
                                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                           <  2:546270/150‑1 (MQ=255)
                                                                                                          tCATAAGGTTGAAGAGAATCTGATTGACAAATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCACCCTGTAGAGGCTTAAGAGTCATTCTTAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                          <  1:255224/150‑1 (MQ=255)
                                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                          >  2:707630/1‑150 (MQ=255)
                                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                          <  1:1636805/150‑1 (MQ=255)
                                                                                                            aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTc                                        <  2:175012/150‑1 (MQ=255)
                                                                                                             aaaGGTTGAACATAAACGGAATGTAAAATCGTGTAGTTTCTAATATGGCCCTGATACCTCACGGCCCCTGGAGATGCTTCACATTGATTTAAAAAAACCAAACATTGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                       <  2:882532/149‑1 (MQ=255)
                                                                                                                  ttGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt                                  >  1:203768/1‑150 (MQ=255)
                                                                                                                   tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                 >  1:1176550/1‑150 (MQ=255)
                                                                                                                    gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTGCATAATAACCCACCATGGACAAGAAAAGTAAGATTAAATCTCTGCAGAGACAATTAACTCAAATGtct                                >  2:320663/1‑148 (MQ=255)
                                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                             <  2:898811/150‑1 (MQ=255)
                                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAGACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAACTCAAATGTataat                             >  2:1455826/1‑150 (MQ=255)
                                                                                                                             cGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGc                       <  1:664651/150‑1 (MQ=255)
                                                                                                                              ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCa                      <  1:1346813/150‑1 (MQ=255)
                                                                                                                                  ggACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTc                  <  2:1170038/150‑1 (MQ=255)
                                                                                                                                   gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                 >  1:379151/1‑150 (MQ=255)
                                                                                                                                   gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                 >  2:1026294/1‑150 (MQ=255)
                                                                                                                                     cAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAc               >  2:467843/1‑150 (MQ=255)
                                                                                                                                        aTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACAAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa            <  2:1581997/150‑1 (MQ=255)
                                                                                                                                               cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa     >  1:1373630/1‑150 (MQ=255)
                                                                                                                                               cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa     >  1:1627113/1‑150 (MQ=255)
                                                                                                                                                gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa    <  2:1095889/150‑1 (MQ=255)
                                                                                                                                                  gACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATc  <  2:203768/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATC  >  NC_000913/2449459‑2449754

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: