Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,887,147 G→A 5.0% A282V (GCG→GTG)  rnd ← RNase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,887,1470GA 5.0% 311.6 / 3.1 97A282V (GCG→GTG) rndRNase D
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base A (3/2);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGATGCCAGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGA  >  NC_000913/1887000‑1887296
                                                                                                                                                   |                                                                                                                                                     
cgcgATGCCAGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCt                                                                                                                                                     <  1:327076/150‑1 (MQ=255)
cgcgATGCCAGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCt                                                                                                                                                     >  1:657292/1‑150 (MQ=255)
        cAGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTccac                                                                                                                                             >  1:747015/1‑150 (MQ=255)
         aGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTccacc                                                                                                                                            <  2:249160/150‑1 (MQ=255)
          gCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTccacca                                                                                                                                           >  1:1420694/1‑150 (MQ=255)
          gCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTccacca                                                                                                                                           <  1:410838/150‑1 (MQ=255)
             aTTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAgcg                                                                                                                                        <  1:1576357/150‑1 (MQ=255)
               tCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCt                                                                                                                                      <  2:520698/150‑1 (MQ=255)
                cGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTa                                                                                                                                     >  1:1069822/1‑150 (MQ=255)
                      tGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCg                                                                                                                               >  1:775328/1‑150 (MQ=255)
                       gATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGt                                                                                                                              >  1:1672332/1‑150 (MQ=255)
                       gATCTTATGCGTTTAGCTCACGTCAGTAATCAGCGACTTAATCTCTATAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCAACAGCGCTAGCAGCGt                                                                                                                              >  2:1317890/1‑150 (MQ=255)
                        aTCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGtt                                                                                                                             <  1:1016362/150‑1 (MQ=255)
                          cTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGtttt                                                                                                                           >  2:609568/1‑150 (MQ=255)
                            tATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTAc                                                                                                                         >  1:1296272/1‑150 (MQ=255)
                              tGCTTTTCGCTCTCGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTACCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCg                                                                                                                       <  2:399985/150‑1 (MQ=255)
                                cGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTT‑AACGCGTTACGCTAACCCGGCATGTCCATCAGGTTAACCATCCGCAGCCGTAAGGCAACATACGGCAAAGTCAGCGCATTTTCCCACAGCGCAAGCAGCCTTTTACCCTGa                                                                                                                    >  1:16508/1‑150 (MQ=255)
                                 gTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGa                                                                                                                    <  1:1382158/150‑1 (MQ=255)
                                    tCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAg                                                                                                                 >  2:43321/1‑150 (MQ=255)
                                    tCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAg                                                                                                                 <  1:980878/150‑1 (MQ=255)
                                     cGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGc                                                                                                                >  1:1279887/1‑150 (MQ=255)
                                      gCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                                                                                               <  1:609568/150‑1 (MQ=255)
                                         cACGTCAGTAATCAGCGACTTAATCGCTTTACACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATNNNNNGCGGTAAGGCATCTTCCGGCAATGTCAGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGAt                                                                                                            >  2:222792/1‑150 (MQ=255)
                                         cACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAGGCGGAt                                                                                                            >  2:1618928/1‑150 (MQ=255)
                                         cACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGAt                                                                                                            >  1:1423527/1‑150 (MQ=255)
                                           cGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGAttt                                                                                                          >  1:824864/1‑150 (MQ=255)
                                             tCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCg                                                                                                        <  2:494206/150‑1 (MQ=255)
                                              cAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                                                                                       <  2:657292/150‑1 (MQ=255)
                                              cAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                                                                                       <  1:1420187/150‑1 (MQ=255)
                                              cAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                                                                                       >  2:1382227/1‑150 (MQ=255)
                                                gTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTa                                                                                                     >  1:1332092/1‑150 (MQ=255)
                                                           cTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAcc                                                                                          <  2:573945/150‑1 (MQ=255)
                                                             tAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGttt                                                                                                                            >  1:110588/1‑114 (MQ=255)
                                                             tAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGttt                                                                                                                            <  2:110588/114‑1 (MQ=255)
                                                             tAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCa                                                                                        <  1:1500152/150‑1 (MQ=255)
                                                              aaTCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGa                                                                                                <  2:1361729/141‑1 (MQ=255)
                                                              aaTCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGa                                                                                                >  1:1361729/1‑141 (MQ=255)
                                                               aTCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAgg                                                                                      >  2:1285684/1‑150 (MQ=255)
                                                               aTCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCACTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAgg                                                                                      >  2:1484571/1‑150 (MQ=255)
                                                               aTCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGCTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCACTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAgg                                                                                      >  2:1478566/1‑150 (MQ=255)
                                                                  gCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTg                                                                                   <  2:266151/150‑1 (MQ=255)
                                                                    cttAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTc                                                                                 >  1:250752/2‑150 (MQ=255)
                                                                       taaCGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAGCGGATAAACCCAGGCTGTCCAg                                                                              <  2:16508/149‑1 (MQ=255)
                                                                        aaCGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt                                                                             <  1:548932/150‑1 (MQ=255)
                                                                        aaCGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt                                                                             >  1:267499/1‑150 (MQ=255)
                                                                          cGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTc                                                                           >  2:1062281/1‑150 (MQ=255)
                                                                           gCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGATGTCCAGTTCg                                                                          <  1:946605/150‑1 (MQ=255)
                                                                             tttACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGcc                                                                        <  1:1155976/150‑1 (MQ=255)
                                                                             tttACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGcc                                                                        >  2:775509/1‑150 (MQ=255)
                                                                             tttACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGcc                                                                        <  2:775328/150‑1 (MQ=255)
                                                                               tACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTa                                                                      >  2:881498/1‑150 (MQ=255)
                                                                                 cGATAACCCGGCATATCCATCACGTTAAGCAACGGCTACAGAAAGGCATCATCCGGCAAAGTATACACATTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTaaa                                                                    >  2:945078/1‑150 (MQ=255)
                                                                                  gaaaaCCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTaa                                                                     >  2:1204771/4‑148 (MQ=255)
                                                                                  gATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTaa                                                                     <  1:1204771/148‑1 (MQ=255)
                                                                                   aTAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACt                                                                  <  1:413657/150‑1 (MQ=255)
                                                                                   aTAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACt                                                                  <  2:893482/150‑1 (MQ=255)
                                                                                     aaCCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTc                                                                <  2:322343/150‑1 (MQ=255)
                                                                                     aaCCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTc                                                                >  2:317309/1‑150 (MQ=255)
                                                                                          ggCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGc                                                           <  2:1406160/150‑1 (MQ=255)
                                                                                           gCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCa                                                          <  1:43321/150‑1 (MQ=255)
                                                                                              tGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCatat                                                       >  2:546678/1‑150 (MQ=255)
                                                                                              tGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCatat                                                       >  2:1448763/1‑150 (MQ=255)
                                                                                                   aTCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAAcgc                                                  <  2:1279887/150‑1 (MQ=255)
                                                                                                     cAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACgcgc                                                <  1:53172/150‑1 (MQ=255)
                                                                                                      aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCt                                               >  2:1192522/1‑150 (MQ=255)
                                                                                                        gTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTAc                                             >  2:825037/1‑150 (MQ=255)
                                                                                                             gCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGAcc                                        <  1:1394190/150‑1 (MQ=255)
                                                                                                             gCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGAcc                                        <  1:1244450/150‑1 (MQ=255)
                                                                                                              cATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCCCCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACca                                       >  1:736600/1‑150 (MQ=255)
                                                                                                               aTCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACTCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACcac                                      >  2:972296/1‑150 (MQ=255)
                                                                                                                tCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACcaca                                     >  2:1118994/1‑150 (MQ=255)
                                                                                                                tCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACcaca                                     >  2:1232047/1‑150 (MQ=255)
                                                                                                                 cGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACaa                                    >  2:574546/1‑150 (MQ=255)
                                                                                                                     tGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGc                                <  2:1644958/150‑1 (MQ=255)
                                                                                                                      gCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCCTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGca                               >  1:1523734/1‑149 (MQ=255)
                                                                                                                           aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTca                          >  2:144129/1‑150 (MQ=255)
                                                                                                                           aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTca                          <  1:1285684/150‑1 (MQ=255)
                                                                                                                           aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTca                          <  2:747015/150‑1 (MQ=255)
                                                                                                                             ggCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTcacg                        >  2:1668196/1‑150 (MQ=255)
                                                                                                                             ggCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTcacg                        >  1:1361055/1‑150 (MQ=255)
                                                                                                                              gCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTcacgc                       <  1:1317890/150‑1 (MQ=255)
                                                                                                                               cATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTcacgca                      >  2:467009/1‑150 (MQ=255)
                                                                                                                                 tCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTcacgcacg                    <  1:467009/150‑1 (MQ=255)
                                                                                                                                   ttCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGAc                  <  2:332952/150‑1 (MQ=255)
                                                                                                                                   ttCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGAc                  >  2:922511/1‑150 (MQ=255)
                                                                                                                                        gCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGt             <  2:1423527/150‑1 (MQ=255)
                                                                                                                                        gCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGt             <  1:881498/150‑1 (MQ=255)
                                                                                                                                        gCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGt             <  2:824864/150‑1 (MQ=255)
                                                                                                                                          aaTGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTc           <  1:1062281/150‑1 (MQ=255)
                                                                                                                                          aaTGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTc           <  1:1200094/150‑1 (MQ=255)
                                                                                                                                           aTGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCa          <  1:1378530/150‑1 (MQ=255)
                                                                                                                                            tGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCAc         <  2:1420694/150‑1 (MQ=255)
                                                                                                                                            tGTCTGCGCTTTTTACACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCAc         >  2:1643932/1‑150 (MQ=255)
                                                                                                                                             gTCTGCACTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCAcc        <  1:1478566/150‑1 (MQ=255)
                                                                                                                                             gTCTGCACTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCAcc        <  1:1484571/150‑1 (MQ=255)
                                                                                                                                               cTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGc      >  2:1388253/1‑150 (MQ=255)
                                                                                                                                                tGCGCTTTTTCCACCAGCGCTAGCGGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATACACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACTCACGACAAAGTTCACCGcc     >  2:1067192/1‑150 (MQ=255)
                                                                                                                                                   gcTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTAGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGa  >  2:779558/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCGATGCCAGCAATTCGGCGCTGATCTTATGCGTTTCGCTCACGTCAGTAATCAGCGACTTAATCGCTTTAAACGCTTTACGATAACCCGGCATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTAAACTTCCCGGCATATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGA  >  NC_000913/1887000‑1887296

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: ‑