Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,291,581 A→C 5.6% L277R (CTC→CGC) ‡ ccmH ← holocytochrome c synthetase ‑ thiol:disulfide oxidoreductase CcmH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,5810AC 5.6% 322.0 / 8.1 104L277R (CTC→CGC) ‡ccmHholocytochrome c synthetase ‑ thiol:disulfide oxidoreductase CcmH
Reads supporting (aligned to +/- strand):  ref base A (48/50);  new base C (3/3);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 106 total aligned reads displayed.

CGTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAAC  >  NC_000913/2291432‑2291730
                                                                                                                                                     |                                                                                                                                                     
cGTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGa                                                                                                                                                       >  2:1020830/1‑150 (MQ=255)
 gTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAg                                                                                                                                                      >  2:166061/1‑150 (MQ=255)
    ggCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCg                                                                                                                                                   <  2:1091140/150‑1 (MQ=255)
      cAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGt                                                                                                                                                 <  2:595384/150‑1 (MQ=255)
        ggTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg                                                                                                                                               <  1:170486/150‑1 (MQ=255)
          tAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTc                                                                                                                                             <  2:485027/150‑1 (MQ=255)
           aaGAGTTTCACCATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCTCGGCGGtgaact                                                                                                                                            <  2:1616761/150‑6 (MQ=255)
           aaGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCTCGGCGGtgaact                                                                                                                                            <  1:1475668/150‑6 (MQ=255)
           aaGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCg                                                                                                                                            <  1:1788204/150‑1 (MQ=255)
            agagTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCTCGGCGGtgaactg                                                                                                                                           <  2:1022866/150‑7 (MQ=255)
               gTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                                                                                                        >  2:1781758/1‑150 (MQ=255)
               gTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                                                                                                        <  1:140703/150‑1 (MQ=255)
                 ttCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGa                                                                                                                                      >  2:1137153/1‑150 (MQ=255)
                   cAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATc                                                                                                                                    <  2:1859662/150‑1 (MQ=255)
                    aaCATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCa                                                                                                                                   <  1:1137153/150‑1 (MQ=255)
                     aCATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAg                                                                                                                                  <  1:1971971/150‑1 (MQ=255)
                        tcatcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATg                                                                                                                               <  1:1851499/150‑1 (MQ=255)
                            ctcCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACg                                                                                                                           <  2:1007408/150‑1 (MQ=255)
                            ctcCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACg                                                                                                                           <  1:1854500/150‑1 (MQ=255)
                             tcCCACGCTGCAACGGCTTAGCCAAATCTCTGCTGCTCAAAGGCATTAAACGCATACATGCATATCACACGCATAATGCTATGGTCCGTTCTCACCAGCTTACGTAGCAGTTCACCGCCGAGCCGGTTGTCGTTGGGATAAGATGAACGa                                                                                                                          >  2:600797/1‑150 (MQ=255)
                              cccACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAg                                                                                                                         >  1:1597713/1‑150 (MQ=255)
                              cccACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAg                                                                                                                         >  1:1892870/1‑150 (MQ=255)
                                cccgcgGCACCGGCTTCGCCACATCGCTGCTGCGCAACGGCATTAAACGCATACCGGCTTAGCCCACGGATATTGCTATGGTCCGTTCTCCCCAGCGGACGTAGCAGTTCACCGCCGAGGCGGGGGTCGGTGGGCGCCGATGACCGAGta                                                                                                                       >  1:1509245/3‑149 (MQ=255)
                                       cAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCtt                                                                                                                >  1:1068512/1‑150 (MQ=255)
                                          cGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAg                                                                                                             <  1:1781758/150‑1 (MQ=255)
                                          cGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAg                                                                                                             <  1:1595808/150‑1 (MQ=255)
                                            gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                                                                                                           <  1:1751134/150‑1 (MQ=255)
                                            gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                                                                                                           >  1:1076795/1‑150 (MQ=255)
                                            gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                                                                                                           <  2:85994/150‑1 (MQ=255)
                                            gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                                                                                                           <  1:487855/150‑1 (MQ=255)
                                            gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                                                                                                           <  1:1327788/150‑1 (MQ=255)
                                             cTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGt                                                                                                          >  1:530847/1‑150 (MQ=255)
                                                  ccAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATccc                                                                                                     <  1:340879/150‑1 (MQ=255)
                                                      aTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTg                                                                                                 >  1:1255469/1‑150 (MQ=255)
                                                      aTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTg                                                                                                 <  2:258479/150‑1 (MQ=255)
                                                       tCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgc                                                                                                <  1:1958165/150‑1 (MQ=255)
                                                           tgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagca                                                                                            <  2:204152/150‑1 (MQ=255)
                                                            gctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCAt                                                                                           >  1:1741441/1‑150 (MQ=255)
                                                             ctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATc                                                                                          <  2:1768564/150‑1 (MQ=255)
                                                                      ggCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                                                                                                        <  2:576623/95‑1 (MQ=255)
                                                                      ggCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                                                                                                        >  1:576623/1‑95 (MQ=255)
                                                                        cATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCg                                                                               <  1:820455/150‑1 (MQ=255)
                                                                         aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCgg                                                                              <  1:312432/150‑1 (MQ=255)
                                                                         aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCgg                                                                              >  2:1330369/1‑150 (MQ=255)
                                                                         aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCgg                                                                              >  1:123587/1‑150 (MQ=255)
                                                                          ttAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGa                                                                             <  1:1227440/150‑1 (MQ=255)
                                                                          ttAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGa                                                                             <  1:1220043/150‑1 (MQ=255)
                                                                             aaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCg                                                                          >  1:1338389/1‑150 (MQ=255)
                                                                             aaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCg                                                                          <  2:467435/150‑1 (MQ=255)
                                                                             aaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCg                                                                          <  2:1290414/150‑1 (MQ=255)
                                                                                 cATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCATCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGc                                                                      <  2:1068512/150‑1 (MQ=255)
                                                                                 cATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGc                                                                      >  1:1744374/1‑150 (MQ=255)
                                                                                   tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTg                                                                                      <  1:760840/132‑1 (MQ=255)
                                                                                   tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTg                                                                                      >  2:760840/1‑132 (MQ=255)
                                                                                   tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGa                                                                    >  2:284073/1‑150 (MQ=255)
                                                                                    aCATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGAt                                                                   >  2:588261/1‑150 (MQ=255)
                                                                                    aCATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGAt                                                                   <  2:1892870/150‑1 (MQ=255)
                                                                                     cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATa                                                                  >  1:468241/1‑150 (MQ=255)
                                                                                       tGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACg                                                                >  1:1836570/1‑150 (MQ=255)
                                                                                         cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGa                                                                    <  2:889625/144‑1 (MQ=255)
                                                                                         cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGa                                                                    >  1:889625/1‑144 (MQ=255)
                                                                                            aGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTg                                                           >  2:639726/1‑150 (MQ=255)
                                                                                               acacGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCg                                                        >  1:709862/1‑150 (MQ=255)
                                                                                                  cGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCg                                                                               <  1:635676/124‑1 (MQ=255)
                                                                                                  cGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCg                                                                               >  2:635676/1‑124 (MQ=255)
                                                                                                  cGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTAt                                                     >  2:1283451/1‑150 (MQ=255)
                                                                                                        ttGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGGTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCg                                               <  2:1282956/150‑1 (MQ=255)
                                                                                                             aTGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGc                                          <  1:166061/150‑1 (MQ=255)
                                                                                                             aTGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGc                                          <  2:1689877/150‑1 (MQ=255)
                                                                                                               ggTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCg                                         >  1:1860379/1‑149 (MQ=255)
                                                                                                               ggTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCg                                         <  2:1860379/149‑1 (MQ=255)
                                                                                                                   cGTTCTCACCAGCTGACGTAGCAGTTAACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACt                                    >  2:140699/1‑150 (MQ=255)
                                                                                                                    gTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTg                                   >  2:1688104/1‑150 (MQ=255)
                                                                                                                        tcACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGt                               >  2:1386953/1‑150 (MQ=255)
                                                                                                                         cACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGtt                              >  2:350541/1‑150 (MQ=255)
                                                                                                                         cACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGtt                              <  2:1338389/150‑1 (MQ=255)
                                                                                                                          aCCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCgctg                             >  2:671263/1‑147 (MQ=255)
                                                                                                                            cAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTAcc                           <  1:1020830/150‑1 (MQ=255)
                                                                                                                                tGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAgc                       >  2:1620025/1‑150 (MQ=255)
                                                                                                                                 gACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAgcg                      <  1:762737/150‑1 (MQ=255)
                                                                                                                                  aCGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAgcgc                     <  2:1628602/150‑1 (MQ=255)
                                                                                                                                   cGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGc                                          >  2:18264/1‑128 (MQ=255)
                                                                                                                                   cGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGc                                          <  1:18264/128‑1 (MQ=255)
                                                                                                                                    gTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCa                   <  1:1330369/150‑1 (MQ=255)
                                                                                                                                     tAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCAt                  >  1:658525/1‑150 (MQ=255)
                                                                                                                                     tAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCAt                  >  1:665122/1‑150 (MQ=255)
                                                                                                                                       gCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGc                <  1:284073/150‑1 (MQ=255)
                                                                                                                                       aacgacaaCCGCCTCGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGc                >  1:1878606/8‑150 (MQ=255)
                                                                                                                                        cAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGcc               >  1:643538/1‑150 (MQ=255)
                                                                                                                                          gacaaCCGCCTCGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCaa             >  1:898758/5‑150 (MQ=255)
                                                                                                                                          gTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCaa             >  1:39171/1‑150 (MQ=255)
                                                                                                                                           ttCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAc            >  1:633119/1‑150 (MQ=255)
                                                                                                                                           ttCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAc            <  2:1597713/150‑1 (MQ=255)
                                                                                                                                            tCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAcg           >  2:835478/1‑150 (MQ=255)
                                                                                                                                            caaCCGCCTCGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAcg           >  2:1464883/3‑150 (MQ=255)
                                                                                                                                             cACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAcgc          >  2:721598/1‑150 (MQ=255)
                                                                                                                                             cACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAAcgc          >  2:815732/1‑150 (MQ=255)
                                                                                                                                                cgccgAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGc       >  2:1408730/1‑150 (MQ=255)
                                                                                                                                                 gccgAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGcc      >  2:813650/1‑150 (MQ=255)
                                                                                                                                                     aGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAAc  <  1:600797/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAAC  >  NC_000913/2291432‑2291730

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: