Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,594 G→T 7.9% T168T (ACC→ACA secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,5940GT 7.9% 457.6 / 22.9 152T168T (ACC→ACAsecYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (67/73);  new base T (7/5);  total (74/78)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 155 total aligned reads displayed.

AGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGGTGT  >  NC_000913/3443445‑3443742
                                                                                                                                                     |                                                                                                                                                    
aGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATg                                                                                                                                                      >  2:277553/1‑150 (MQ=255)
       gCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCGCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtg                                                                                                                                               >  1:1686297/1‑150 (MQ=255)
       gCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtg                                                                                                                                               <  1:898915/150‑1 (MQ=255)
       gCAGGTCGCCTTGACGCGCTTGCGCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtg                                                                                                                                               <  2:618993/150‑1 (MQ=255)
          ggTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtgatg                                                                                                                                            <  2:1364940/150‑4 (MQ=255)
          ggTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtgatg                                                                                                                                            <  2:212665/150‑4 (MQ=255)
           gTCGCCTTGACGCGCTTGCTCGATAGTATGTGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGAc                                                                                                                                           <  1:106090/150‑1 (MQ=255)
            tCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGAcc                                                                                                                                          >  1:1277542/1‑150 (MQ=255)
              gCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGCCCac                                                                                                                                        >  1:852636/1‑149 (MQ=255)
              gCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAg                                                                                                                                        >  1:1101583/1‑150 (MQ=255)
               ccTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGa                                                                                                                                       >  2:1210835/1‑150 (MQ=255)
                     cgcgCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGAca                                                                                                                                          <  2:1090079/141‑2 (MQ=255)
                     cgcgCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTaca                                                                                                                                 <  2:884992/150‑1 (MQ=255)
                      gcgcTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaa                                                                                                                                <  1:1049666/150‑1 (MQ=255)
                       cgcTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaac                                                                                                                               >  1:686225/1‑150 (MQ=255)
                            gCTCGATAGTATGGTCAATGTCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCgg                                                                                                                          <  1:314371/150‑1 (MQ=255)
                            gCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCgg                                                                                                                          <  2:57511/150‑1 (MQ=255)
                            gCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCgg                                                                                                                          >  2:1385899/1‑150 (MQ=255)
                             cTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGt                                                                                                                         >  1:122901/1‑150 (MQ=255)
                              tCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGAACAGACTTCCACCAGCGGTg                                                                                                                        >  2:89143/1‑150 (MQ=255)
                                  tAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGt                                                                                                                    <  2:1570699/150‑1 (MQ=255)
                                    gTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAg                                                                                                                  >  1:684138/1‑150 (MQ=255)
                                        gggCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATg                                                                                                              <  1:260966/150‑1 (MQ=255)
                                        gggCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATg                                                                                                              <  2:431419/150‑1 (MQ=255)
                                                gctggcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTac                                                                                                                                  <  2:240236/122‑1 (MQ=255)
                                                gctggcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTac                                                                                                                                  >  1:240236/1‑122 (MQ=255)
                                                    gcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGt                                                                                                  >  1:10630/1‑150 (MQ=255)
                                                     cGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGtt                                                                                                 >  1:1390062/1‑150 (MQ=255)
                                                     cGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGtt                                                                                                 >  1:418396/1‑150 (MQ=255)
                                                     cGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGtt                                                                                                 <  2:353043/150‑1 (MQ=255)
                                                      gggAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTa                                                                                                <  2:1311258/150‑1 (MQ=255)
                                                          ggcccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCa                                                                                            <  1:142663/148‑1 (MQ=255)
                                                          gTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCATAGCCCGGGTTAATCa                                                                                            <  1:401301/150‑1 (MQ=255)
                                                             ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCa                                                                                         >  2:1115608/1‑150 (MQ=255)
                                                             ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCa                                                                                         <  1:1017609/150‑1 (MQ=255)
                                                                   cAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCtt                                                                                   >  1:410803/1‑150 (MQ=255)
                                                                      tACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCa                                                                                <  2:799927/150‑1 (MQ=255)
                                                                       aCCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCAt                                                                               >  1:173263/1‑150 (MQ=255)
                                                                        ccGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTTATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATa                                                                              <  1:13272/150‑1 (MQ=255)
                                                                        ccGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACAACATCAGGAACATGGTTCCTGTGACCAGACTTACAACATCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATa                                                                              <  1:1115608/150‑1 (MQ=255)
                                                                            cGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAACGCCCGGGTTAATCACCAGGCCTTGCATACCAg                                                                          >  2:946380/1‑150 (MQ=255)
                                                                            cGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAg                                                                          >  2:1624805/1‑150 (MQ=255)
                                                                            cGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAg                                                                          <  1:537861/150‑1 (MQ=255)
                                                                                  tAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCatat                                                                    >  1:737058/1‑150 (MQ=255)
                                                                                  tAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCatat                                                                    >  2:1717783/1‑150 (MQ=255)
                                                                                    aTGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTc                                                                  <  1:1399445/150‑1 (MQ=255)
                                                                                     tGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCg                                                                 >  1:1523941/1‑150 (MQ=255)
                                                                                     tGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCg                                                                 >  2:267422/1‑150 (MQ=255)
                                                                                     tGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCg                                                                 >  2:1560419/1‑150 (MQ=255)
                                                                                      gATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCCCATCAGGAACATGGTTCTTGTGACCAGACTTCCAACAGCGGTGAAGTAGAATGCAAGGCCCGGGTTAAGCACCAGGCCTGGCATACCAGGCATATTCgg                                                                >  2:164063/1‑150 (MQ=255)
                                                                                       aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGc                                                               <  2:1277542/150‑1 (MQ=255)
                                                                                       aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGc                                                               <  1:955511/150‑1 (MQ=255)
                                                                                       aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGc                                                               <  2:758020/150‑1 (MQ=255)
                                                                                       aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGc                                                               <  1:1013549/150‑1 (MQ=255)
                                                                                          gAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGa                                                            >  1:957810/1‑150 (MQ=255)
                                                                                          gAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGa                                                            >  2:1328818/1‑150 (MQ=255)
                                                                                           aaaTACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCCCCAGGCCTTGCATACCAGGCATATTCGGCAGAc                                                           >  1:108342/1‑150 (MQ=255)
                                                                                           aaaTACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGAc                                                           <  2:684138/150‑1 (MQ=255)
                                                                                            aaTACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGAcc                                                          <  1:606825/150‑1 (MQ=255)
                                                                                              tACCGTTGCCGATACCTCGTTCAGTAATCTGGTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTCAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCgg                                                        <  2:937379/150‑1 (MQ=255)
                                                                                                  gTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGGGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGc                                                    <  2:1105218/150‑1 (MQ=255)
                                                                                                  gTTGCCGAAACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGc                                                    <  1:1103568/150‑1 (MQ=255)
                                                                                                     gCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAAt                                                 >  2:1534819/1‑150 (MQ=255)
                                                                                                      ccGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATa                                                <  2:586830/150‑1 (MQ=255)
                                                                                                      ccGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATa                                                >  1:1607074/1‑150 (MQ=255)
                                                                                                        gATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTg                                                                                  <  2:441066/114‑1 (MQ=255)
                                                                                                        gATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTg                                                                                  >  1:441066/1‑114 (MQ=255)
                                                                                                         aTACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATCCCg                                             >  2:762158/1‑150 (MQ=255)
                                                                                                            ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                          >  1:13534/1‑150 (MQ=255)
                                                                                                            ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                          <  1:1210835/150‑1 (MQ=255)
                                                                                                            ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                          >  2:1037756/1‑150 (MQ=255)
                                                                                                            ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                          >  2:819135/1‑150 (MQ=255)
                                                                                                            ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAAAGCAAAGCCCGGGTTAATCACCACGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                          <  1:97664/150‑1 (MQ=255)
                                                                                                              tCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGGGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGa                                        <  1:1218421/150‑1 (MQ=255)
                                                                                                               cGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGAc                                       >  1:1186162/1‑150 (MQ=255)
                                                                                                                gTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAACAGGCCTTGCATACCAGGCATATTAGTCAGACCGGTAGAAATACGGATCGACt                                      >  1:290513/1‑150 (MQ=255)
                                                                                                                  tCAGTAATCTGTTCGCGCAACCACATCAGGAACATGGTTCCTGTGACCAGACTAACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGTCCTTGCATACCAGGCATTTTCGGCAGACCGGTAGCAATACCGATCGACTgg                                    <  2:1580635/150‑1 (MQ=255)
                                                                                                                  tCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTgg                                    >  2:591486/1‑150 (MQ=255)
                                                                                                                   cAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGACATATTCGGCAGACCGGTAGCAATACCGATCGACTGGa                                   <  1:238389/150‑1 (MQ=255)
                                                                                                                   cAGTAATCTGTTCGCCAAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGACATATTCGGCAGACCGGTAGCAATACCGATCGACTGGa                                   >  2:238389/1‑150 (MQ=255)
                                                                                                                    aGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGaa                                  >  1:480774/1‑150 (MQ=255)
                                                                                                                           tGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGTTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGcc                           <  2:1526157/150‑1 (MQ=255)
                                                                                                                            gTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGAt                                           >  2:966241/1‑133 (MQ=255)
                                                                                                                            gTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCa                          <  1:1037756/150‑1 (MQ=255)
                                                                                                                            gTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCa                          <  1:525211/150‑1 (MQ=255)
                                                                                                                                gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCAc                      >  1:790135/1‑150 (MQ=255)
                                                                                                                                gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGAACGACTGGAATATTGCCAGCAc                      <  2:627544/150‑1 (MQ=255)
                                                                                                                                 cccAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCAcc                     <  2:545471/150‑1 (MQ=255)
                                                                                                                                 cccAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCAcc                     >  1:1252256/1‑150 (MQ=255)
                                                                                                                                   cAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCag                   <  1:1534819/150‑1 (MQ=255)
                                                                                                                                   cAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCag                   <  1:1560419/150‑1 (MQ=255)
                                                                                                                                          cacAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACAAGAGTACCg            >  1:1505503/3‑150 (MQ=255)
                                                                                                                                          aTCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCg            >  1:913695/1‑150 (MQ=255)
                                                                                                                                          aTCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCg            >  1:644339/1‑150 (MQ=255)
                                                                                                                                           acAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGt           >  1:1029829/2‑150 (MQ=255)
                                                                                                                                           acAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGAACCGt           >  1:1700992/2‑150 (MQ=255)
                                                                                                                                           acAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATAGCCAGCACCAGAGAACCGt           >  1:255639/2‑150 (MQ=255)
                                                                                                                                             aGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAg         <  1:86446/150‑1 (MQ=255)
                                                                                                                                               gAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCg       >  1:758025/1‑150 (MQ=255)
                                                                                                                                                    tGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGgtgt  >  1:1662140/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGGTGT  >  NC_000913/3443445‑3443742

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