Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,854,381 T→C 5.7% F15S (TTC→TCC)  fhlA → DNA‑binding transcriptional activator FhlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,854,3810TC 5.7% 321.9 / 7.1 105F15S (TTC→TCC) fhlADNA‑binding transcriptional activator FhlA
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base C (3/3);  total (53/52)
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.71e-01

Only 100 of 108 total aligned reads displayed.

CCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCTTATTACGCGTCGC  >  NC_000913/2854232‑2854530
                                                                                                                                                     |                                                                                                                                                     
 cacaCGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTc                                                                                                                                                      >  2:21530/1‑150 (MQ=255)
 cacaCGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTc                                                                                                                                                      >  2:1152507/1‑150 (MQ=255)
 cacaCGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTc                                                                                                                                                      <  1:1581009/150‑1 (MQ=255)
  acacGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCg                                                                                                                                                     <  1:636366/150‑1 (MQ=255)
     cGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACa                                                                                                                                                  >  1:55217/1‑150 (MQ=255)
      gCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACAt                                                                                                                                                 <  2:1326830/150‑1 (MQ=255)
         gAACGGCTTCCGCGTATATGCTAATAAATTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTTCTTTCCTTGAAGAACAAAATGTCATATACACCGATGAGTGATGTCGAACAACCCTTGTTGTCCGAGATCAc                                                                                                                                              <  2:312134/150‑1 (MQ=37)
         gAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCAc                                                                                                                                              <  1:1428733/150‑1 (MQ=255)
             cGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTcat                                                                                                                                          <  2:1229972/150‑3 (MQ=38)
             cGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCgg                                                                                                                                          <  1:1265717/150‑1 (MQ=255)
              gCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGa                                                                                                                                         >  1:588432/1‑150 (MQ=255)
               cTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGac                                                                                                                                        <  1:825950/150‑1 (MQ=255)
                 tCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGacac                                                                                                                                      >  1:41094/1‑150 (MQ=255)
                   cgcgTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTAt                                                                                                                                   >  2:730678/1‑150 (MQ=255)
                     cgTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTAtt                                                                                                                                  >  1:469012/1‑150 (MQ=255)
                      gTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTg                                                                                                                                 <  2:1408307/150‑1 (MQ=255)
                         tatGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTgcag                                                                                                                              <  2:769883/150‑1 (MQ=255)
                             cTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCagcagc                                                                                                                          >  1:352805/1‑150 (MQ=255)
                               aataaAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGccc                                                                                                                        >  2:97939/1‑150 (MQ=255)
                                 taaAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGa                                                                                                                      >  1:990747/1‑150 (MQ=255)
                                    aaTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCt                                                                                                                   >  1:286795/1‑150 (MQ=255)
                                      ttCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTgg                                                                                                                 <  1:1604720/150‑1 (MQ=255)
                                      ttCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTgg                                                                                                                 <  2:1412241/150‑1 (MQ=255)
                                      ttCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTgg                                                                                                                 <  2:1065135/150‑1 (MQ=255)
                                         tAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGccg                                                                                                              >  2:1611143/1‑149 (MQ=255)
                                         tAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCt                                                                                                              >  2:916436/1‑150 (MQ=255)
                                          aaaTCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTc                                                                                                             >  1:576872/1‑150 (MQ=255)
                                              ctcCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtg                                                                                                         <  2:394138/150‑1 (MQ=255)
                                                 cTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGc                                                                                                           <  2:1134930/145‑1 (MQ=255)
                                                 cTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGc                                                                                                           >  1:1134930/1‑145 (MQ=255)
                                                  tataGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                                                                                                     >  1:840608/1‑150 (MQ=255)
                                                  tataGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                                                                                                     <  1:339543/150‑1 (MQ=255)
                                                  tataGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                                                                                                     <  2:142035/150‑1 (MQ=255)
                                                  tataGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                                                                                                     <  1:161347/150‑1 (MQ=255)
                                                  tataGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                                                                                                     <  2:1277914/150‑1 (MQ=255)
                                                     agttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAgg                                                                                                  >  2:1016727/1‑150 (MQ=255)
                                                        tagtCAATGACCTTTTGCACCGCTTTGAGGTGCTTTCATGGAAGAACAAAAAGTCATATACAACGATGAATGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGctc                                                                                               >  2:915031/1‑150 (MQ=255)
                                                         agtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGctct                                                                                              >  1:436332/1‑150 (MQ=255)
                                                          gtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCtt                                                                                             <  1:489623/150‑1 (MQ=255)
                                                           tCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCttt                                                                                            >  2:179186/1‑150 (MQ=255)
                                                            cAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTc                                                                                           >  2:1018697/1‑150 (MQ=255)
                                                             aaTGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCg                                                                                          >  2:1144593/1‑150 (MQ=255)
                                                              aTGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGc                                                                                         <  1:1396420/150‑1 (MQ=255)
                                                               tGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCa                                                                                        >  1:1382851/1‑150 (MQ=255)
                                                               tGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCa                                                                                        <  1:1053474/150‑1 (MQ=255)
                                                                 caggaaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAAc                                                                                      >  1:1135071/8‑150 (MQ=255)
                                                                 caggaaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAAc                                                                                      >  1:1614613/8‑150 (MQ=255)
                                                                  ccTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACt                                                                                     >  1:111116/1‑150 (MQ=255)
                                                                   ggaaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTg                                                                                    >  1:432281/6‑150 (MQ=255)
                                                                   cTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTg                                                                                    <  1:1499604/150‑1 (MQ=255)
                                                                         cACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAg                                                                              <  2:840608/150‑1 (MQ=255)
                                                                         cACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAg                                                                              >  1:863525/1‑150 (MQ=255)
                                                                          aCCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGc                                                                             <  1:1018697/150‑1 (MQ=255)
                                                                          aCCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGc                                                                             <  1:730678/150‑1 (MQ=255)
                                                                              cTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTc                                                                         <  1:1152507/150‑1 (MQ=255)
                                                                                     gTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTc                                                                  >  1:835864/1‑150 (MQ=255)
                                                                                        cTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGcc                                                               <  2:352805/150‑1 (MQ=255)
                                                                                          ttCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGa                                                             <  1:1593120/150‑1 (MQ=255)
                                                                                             cTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACaa                                                          >  2:1116057/1‑150 (MQ=255)
                                                                                             cTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACaa                                                          >  1:1284802/1‑150 (MQ=255)
                                                                                              tGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAAc                                                         >  2:658795/1‑150 (MQ=255)
                                                                                               ggAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGGCAAcg                                                        >  2:1070391/1‑150 (MQ=255)
                                                                                               ggAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAAcg                                                        >  2:1392027/1‑150 (MQ=255)
                                                                                                  aaaaaaaaaAGTCATATACAACAATGAGTCATATAGTACAACAAGGTTTTATCCAAATCAATAAAAAAATATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCgg                                                     >  2:914850/6‑150 (MQ=255)
                                                                                                   gaaCAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGc                                                    >  1:944455/1‑150 (MQ=255)
                                                                                                    aaCAAAATGTCATATACACCGATGAGTGATCTCGGACAACATGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCt                                                   <  2:835864/150‑1 (MQ=255)
                                                                                                     aCAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTa                                                  >  2:687763/1‑150 (MQ=255)
                                                                                                        aaaTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttg                                               <  2:1434192/150‑1 (MQ=255)
                                                                                                        aaaTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttg                                               <  2:286795/150‑1 (MQ=255)
                                                                                                        aaaTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttg                                               <  2:249758/150‑1 (MQ=255)
                                                                                                        aaaTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttg                                               <  2:576872/150‑1 (MQ=255)
                                                                                                          aTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttgtg                                             >  1:433482/1‑150 (MQ=255)
                                                                                                          aTGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTAttgtg                                             <  1:687763/150‑1 (MQ=255)
                                                                                                                 ataCACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGc                                      <  1:411543/150‑1 (MQ=255)
                                                                                                                 ataCACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGc                                      <  1:1110884/150‑1 (MQ=255)
                                                                                                                 ataCACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGACCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGc                                      <  1:494645/150‑1 (MQ=255)
                                                                                                                  taCACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCa                                     >  2:1419024/1‑150 (MQ=255)
                                                                                                                   acacCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCaa                                    >  1:821903/1‑150 (MQ=255)
                                                                                                                      tgtccgAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAgcg                                 >  2:125762/6‑150 (MQ=37)
                                                                                                                      ccGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAgcg                                 <  2:1135071/150‑1 (MQ=255)
                                                                                                                       cGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAgcgc                                >  1:473876/1‑150 (MQ=255)
                                                                                                                          tGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGa                             >  2:22625/1‑150 (MQ=255)
                                                                                                                          cgAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAgcgc                                <  2:1623139/146‑1 (MQ=38)
                                                                                                                          cgAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAgcgc                                >  1:1623139/2‑147 (MQ=38)
                                                                                                                             gTGATGTCGAACAACCCTTGTTGTCCGAGATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTc                          >  1:105656/1‑150 (MQ=38)
                                                                                                                               gATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCaa                        <  1:1016727/150‑1 (MQ=255)
                                                                                                                               gATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCaa                        <  1:689065/150‑1 (MQ=255)
                                                                                                                                 tctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACg                      <  1:1392027/150‑1 (MQ=255)
                                                                                                                                 tctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACg                      <  2:407641/150‑1 (MQ=255)
                                                                                                                                 tctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACg                      <  2:55217/150‑1 (MQ=255)
                                                                                                                                 tctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACg                      >  1:425246/1‑150 (MQ=255)
                                                                                                                                  ctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGt                     >  1:1603602/1‑150 (MQ=255)
                                                                                                                                  ctcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGt                     >  2:1185538/1‑150 (MQ=255)
                                                                                                                                   tcGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCATCCCCATCTGGCCTCGCTGTGTGAGGCTCGTTCGCAACTGGTAACGCGTTCTGCGCTCCCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACgag                    >  2:123274/1‑148 (MQ=255)
                                                                                                                                     ggACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAAGTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCg                  <  1:1231066/150‑1 (MQ=255)
                                                                                                                                      gACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGt                 >  1:948318/1‑150 (MQ=255)
                                                                                                                                          acaaGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCtta             >  2:1266047/1‑150 (MQ=255)
                                                                                                                                          acaaGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCtta             >  1:1174291/1‑150 (MQ=255)
                                                                                                                                            aaGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCttatt           <  2:863525/150‑1 (MQ=255)
                                                                                                                                                     tCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCTTATTAcgcgtcgc  >  2:688936/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCGTTCTGCGCTCGCCGACAACGCGGCTATTGTGTTGTGGCAAGCGCAGACTCAACGTGCGTCTTATTACGCGTCGC  >  NC_000913/2854232‑2854530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: