Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 639,654 T→C 5.2% intergenic (+146/‑99) ahpC → / → ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913639,6540TC 5.2% 460.8 / 11.1 153intergenic (+146/‑99)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (70/74);  new base C (4/4);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-02

Only 100 of 159 total aligned reads displayed.

TAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTG  >  NC_000913/639506‑639803
                                                                                                                                                    |                                                                                                                                                     
tAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGtgctg                                                                                                                                                      <  2:1557770/150‑1 (MQ=255)
 aaaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGtgctgc                                                                                                                                                     <  2:1041373/150‑1 (MQ=255)
     gtCCTTCGTCTTTCACGCGCTAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCTCTTACTTGTGTAAGCTCCCTGGGATTCACCGCCTAGCCGCCTTGCTCTGACGCGAAATACTGCGGCAATGCACCTAATTCTTCGGGTGCTTCGGCa                                                                                                                                                 <  2:297781/149‑1 (MQ=255)
      tCCTTCTTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGGCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                                                                                                                <  2:563881/150‑1 (MQ=255)
       ccTTCGTCTTTCACGCCATAGCGGCGTTTGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCAGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAcc                                                                                                                                               >  2:1766148/1‑150 (MQ=255)
       ccTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAcc                                                                                                                                               >  1:959202/1‑150 (MQ=255)
        cTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCAccc                                                                                                                                              <  1:1774408/150‑1 (MQ=255)
           cGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAt                                                                                                                                           >  1:49467/1‑150 (MQ=255)
            ggcTTTCACGCCATAGCGGCGGTGGCGGCGCCCGCTCACCCCGGGCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCGCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaag                                                                                                                                          <  2:330081/148‑3 (MQ=255)
             tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAttt                                                                                                                                         <  1:389234/150‑1 (MQ=255)
              cTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCGCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaagaa                                                                                                                                        <  1:1314899/150‑5 (MQ=255)
              cTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtttt                                                                                                                                        >  1:1297263/1‑150 (MQ=255)
                ttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTATTTCTTCGGGTGCTGCGGCACCCGATTttct                                                                                                                                      <  1:1004775/150‑1 (MQ=255)
                ttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttct                                                                                                                                      <  2:914592/150‑1 (MQ=255)
                ttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttct                                                                                                                                      <  1:1730671/150‑1 (MQ=255)
                ttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCGGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCCCCCGAGTttcg                                                                                                                                      >  1:1219879/1‑149 (MQ=255)
                   ccGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                   <  1:773959/149‑1 (MQ=255)
                   aCGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                   >  2:41377/1‑150 (MQ=255)
                   aCGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                   <  1:767185/150‑1 (MQ=255)
                     cccATAGCGGCGTTGGCGTAGCCCGCTCCCCCCGGTCACTGACATGTGTAAGCTCCCGGGGATTCACAGGCTCGCACCCTAGCTCTGACGCGAAATACTTCGGAAATTCACATAATTCTTCGGGTGCTGCGGCAACCGATGttctgcacc                                                                                                                                 >  2:306165/2‑145 (MQ=255)
                      ccATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTccccc                                                                                                                                >  2:1842096/1‑149 (MQ=255)
                      ccATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCg                                                                                                                                <  1:841656/150‑1 (MQ=255)
                          aGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCAcc                                                                                                                            <  1:1631431/150‑1 (MQ=255)
                           ccggcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCCCCa                                                                                                                           >  1:90377/2‑150 (MQ=255)
                             ggcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCatg                                                                                                                         >  2:1313262/1‑150 (MQ=255)
                             ggcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCatg                                                                                                                         >  1:1603142/1‑150 (MQ=255)
                              gcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAGTACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCatga                                                                                                                        <  1:1026168/150‑1 (MQ=255)
                                gTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGCTTCACAGGCTAGCCGCCTTGCTCTGGCGCGAAATACTTCGGAATTTCACCTAATTCTTCGGGTGCTGCGGCACCCGAGGTTCTTCCCCGCACCatgatg                                                                                                                      >  1:1231303/1‑150 (MQ=255)
                                   ggCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCaa                                                                                                                   >  1:995411/1‑150 (MQ=255)
                                     cGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGACTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttct                                                                                                                                      >  1:1855257/1‑129 (MQ=255)
                                      ggcGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCt                                                                                                                <  2:187251/148‑1 (MQ=255)
                                        cGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAGTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                              >  2:1167881/1‑150 (MQ=255)
                                         gCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCa                                                                                                             <  2:188068/150‑1 (MQ=255)
                                          cccGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATTATTCAAGCTGCAt                                                                                                            >  2:1501174/1‑150 (MQ=255)
                                          cccGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGAGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCAt                                                                                                            <  2:1448701/150‑1 (MQ=255)
                                           ccGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGTCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATc                                                                                                           <  1:1670486/150‑1 (MQ=255)
                                             gCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCa                                                                                                         <  2:238647/150‑1 (MQ=255)
                                              cTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg                                                                                                        <  1:1498877/150‑1 (MQ=255)
                                               tCACCCCCGTCACTTACTTGTGTCAGCTCCCTGGGATTCACAGGCTAGCCGCCTTGCTCGGACGCGAAATACTTCGGACATTCACCTACTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAGCTTGCATcatgg                                                                                                       <  2:1236236/150‑5 (MQ=255)
                                                 aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTa                                                                                                     >  1:838564/1‑150 (MQ=255)
                                                  ccccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAg                                                                                                    <  1:261109/150‑1 (MQ=255)
                                                   cccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGc                                                                                                   <  2:90377/150‑1 (MQ=255)
                                                   ccaggTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAAATCACCTAATTCTTCGGGTGCTGCGGAACCCGATTTTCTTCCCCGCACCATGATGCAAGCTTCATCAAGGTAGc                                                                                                   >  2:1382143/4‑150 (MQ=255)
                                                        tCACTTACTTGTGTAAGCTCCCGGTGATTCACAGGCTAGCCTCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCag                                                                                              <  2:1298687/150‑1 (MQ=255)
                                                         cACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCGAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTGCCCGCATCATGATGCAAGCTGCATCCAGGTAGCCGCagc                                                                                             >  2:741773/1‑149 (MQ=255)
                                                          acttacttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCagag                                                                                            >  1:1439887/1‑150 (MQ=255)
                                                           cttacttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAgg                                                                                           >  1:1289859/1‑150 (MQ=255)
                                                             tacttGTGTAAGCTCCCGGGGATTCACCGGCTAGCCGCCGTGCTCTGACGCGAAATACTTCGGAAATTAACCTAATTCTTCGGGGGCTGAGGCACCCGAGTTTCGTCCCCGCACCAGGCAGCAACCACCAACCCGGTAGCCGCCGCGGcc                                                                                         >  2:51638/1‑150 (MQ=255)
                                                               cttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGc                                                                                       >  2:954647/1‑150 (MQ=255)
                                                                    gtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCTCCGCACCATGATGCATGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCa                                                                                  >  1:1708185/1‑150 (MQ=255)
                                                                        gCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgat                                                                              <  1:1337857/150‑1 (MQ=255)
                                                                         cTCCCGGGGATTCACTGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatg                                                                             <  2:1433352/150‑1 (MQ=255)
                                                                         cTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGCGGCCGCTTGCatgatg                                                                             >  1:590842/1‑150 (MQ=255)
                                                                          gcccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatga                                                                            <  2:1602656/149‑1 (MQ=255)
                                                                           cccGGGGATTCACAGGCTAGCCGCCTTGCTCTGCCGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCCCCCGATTTTCTTCCCCGCCCCATGATGCAAGCTGCATCCAGGTGGCCGCAGAGGCCGCTTGCatgatgat                                                                           >  1:151209/1‑150 (MQ=255)
                                                                           cccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCCAGCTGCTTCCAGGTAGCCGCCGCGGCCGCTGGCatgatgct                                                                           >  1:732547/1‑148 (MQ=255)
                                                                             cGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                         >  2:688315/1‑150 (MQ=255)
                                                                             cGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                         >  2:1514966/1‑150 (MQ=255)
                                                                             cGGGGATTCACAGGCTAGCCGCCTTGCGCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                         <  2:750565/150‑1 (MQ=255)
                                                                              ggggATTCACAGGCTCGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCAGCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGtt                                                                        <  2:988552/150‑1 (MQ=255)
                                                                                 gATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTaa                                                                     <  2:1791148/150‑1 (MQ=255)
                                                                                 gATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTaa                                                                     <  2:480451/150‑1 (MQ=255)
                                                                                   ttCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAg                                                                   >  1:1536114/1‑150 (MQ=255)
                                                                                   ttCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAg                                                                   <  1:1131352/150‑1 (MQ=255)
                                                                                    tCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAATTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGAATCCAGGTAGCCGCAGATGCCGCTGGCATGATGATGTTTaaatc                                                                  >  2:1813456/1‑148 (MQ=255)
                                                                                    tCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTACCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGGTTAAAGc                                                                  >  1:1854179/1‑150 (MQ=255)
                                                                                         ggCTAGCCGCCTTGCTCTGACGCGAAATACTTCTGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAgg                                                             <  2:210239/150‑1 (MQ=255)
                                                                                         ggCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGGTGCAAGCTGCATCCAGGTAGCCGCAGCGGCCGCTTGCATGATGATGTTTAAAGCCCAgg                                                             >  1:893048/1‑150 (MQ=255)
                                                                                          gCTAGCCGCCTTGCGCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGga                                                            <  1:339047/150‑1 (MQ=255)
                                                                                           cTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGgag                                                           >  1:1511380/1‑150 (MQ=255)
                                                                                             aGCCGCCTAGCCCTGACCCGAAATACATCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATATTCTACCCCGCACCATGATGCCCGCGGCATCCAGGTAGCCGCCGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGAt                                                         <  2:1334762/150‑1 (MQ=255)
                                                                                               ccgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaa                                                       <  2:667062/150‑1 (MQ=255)
                                                                                                   cTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACAt                                                   >  1:736036/1‑150 (MQ=255)
                                                                                                        tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCg                                              <  2:668341/150‑1 (MQ=255)
                                                                                                        tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCAACCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTACAGCCCAGGAGATAAACATGCTCg                                              >  1:624723/1‑150 (MQ=255)
                                                                                                         ctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGGATGATGATGTTTAAAGCACAGGTGATAAACATGCTCgc                                             >  2:1737408/1‑149 (MQ=255)
                                                                                                          ttaCGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                            <  2:1401645/148‑1 (MQ=255)
                                                                                                          tGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                            >  2:1057447/1‑150 (MQ=255)
                                                                                                          tGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                            >  2:578469/1‑150 (MQ=255)
                                                                                                           gACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGaca                                           >  1:990057/1‑150 (MQ=255)
                                                                                                            aCGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGacac                                          <  2:959202/150‑1 (MQ=255)
                                                                                                            aCGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCGTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGacac                                          >  1:1120266/1‑150 (MQ=255)
                                                                                                             cgcgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGAcaca                                         >  1:1474131/1‑150 (MQ=255)
                                                                                                               cgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACaaa                                       <  1:659931/150‑1 (MQ=255)
                                                                                                                      cTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAACGCCCAGGAGATAAACATGCTCGACACAAATATGaaa                                >  2:28001/1‑150 (MQ=255)
                                                                                                                       ttCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATGCAGGTAGCCGCAGAGGCCGCGTGCAGGAGGATGTTTAACGACCAGGAGATAAACATGAGCGCCACAAATATGtaaa                               >  2:334604/1‑150 (MQ=255)
                                                                                                                        tCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaac                              <  1:201511/150‑1 (MQ=255)
                                                                                                                           gAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactca                           <  1:743179/150‑1 (MQ=255)
                                                                                                                           gAAATTCACCTAATGCTTCGCGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTTCATCCAGGTAGCCGCAGAGGCAGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactca                           <  1:285684/150‑1 (MQ=255)
                                                                                                                            aaaTTCACCTAATTCTTCTCGTGCTGCGGCACCCGATTTTCTTACCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTTCATGATTATGTTTTAATCCCAGGAGATAAACAAGCTCGACACAAATATGAAaactcaa                          >  2:1223613/1‑150 (MQ=255)
                                                                                                                            aaaTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                          >  2:1862606/1‑150 (MQ=255)
                                                                                                                              aTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaacg                        >  1:1340424/1‑149 (MQ=255)
                                                                                                                                       aTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTa               <  2:1289559/150‑1 (MQ=255)
                                                                                                                                         aaatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTAcc             >  2:612403/4‑150 (MQ=255)
                                                                                                                                           ttcGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCtt           >  1:702406/1‑150 (MQ=255)
                                                                                                                                              gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACTCAAATATGAAAACTCAACTCAAGGCTTACCTTgag        >  2:1613399/1‑150 (MQ=255)
                                                                                                                                               ggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCATTTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAATGTCAACTCAAGGCTTACCTTgaga       <  1:1570316/150‑1 (MQ=255)
                                                                                                                                               ggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTgaga       >  2:385833/1‑150 (MQ=255)
                                                                                                                                                    tgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTg  >  2:1282171/1‑150 (MQ=255)
                                                                                                                                                    tgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTg  >  2:590167/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTG  >  NC_000913/639506‑639803

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: