Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,658,422 C→G 6.5% A172P (GCG→CCG)  hscA ← DnaK‑like molecular chaperone specific for IscU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,658,4220CG 6.5% 325.6 / 11.7 92A172P (GCG→CCG) hscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base C (41/45);  new base G (4/2);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

TAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTC  >  NC_000913/2658276‑2658564
                                                                                                                                                  |                                                                                                                                              
tAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                                                                                                                               >  2:839702/1‑148 (MQ=255)
tAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                                                                                                                               >  1:318771/1‑148 (MQ=255)
  aGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcggcgc                                                                                                                                             >  2:417166/1‑148 (MQ=255)
       aGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCttt                                                                                                                                       <  3:477399/149‑1 (MQ=255)
       aGAATGGAAATATCAAACGCCCCGCCACCGAGGTCATAAAAGGCGATCACACCTTCCTGACCGGAATCCAGCCCGTAGGCGAGAGCCGCATCGGTCGGTTCGTTAAGTAACCGCAGGACGTGGAGGCCCGCCTGACGCGCCGCGTCttt                                                                                                                                       <  3:760764/149‑1 (MQ=255)
        gAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtgg                                                                                                                                       <  1:417166/148‑3 (MQ=255)
        gAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTg                                                                                                                                      <  3:715124/149‑1 (MQ=255)
          aTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                    <  3:41661/149‑1 (MQ=255)
             gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGc                                                                                                                                  >  2:144684/1‑148 (MQ=255)
               aaTATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                                                                                                                                          >  2:339694/1‑138 (MQ=255)
               aaTATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                                                                                                                                          <  1:339694/138‑1 (MQ=255)
               aaTATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                >  1:727319/1‑148 (MQ=255)
                 tatCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                                                                                                                             <  3:200544/149‑1 (MQ=255)
                      aaCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg                                                                                                                        >  4:688287/1‑149 (MQ=255)
                          tCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc                                                                                                                     >  2:838277/1‑148 (MQ=255)
                            ccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc                                                                                                                  >  3:837793/1‑149 (MQ=255)
                              gCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGt                                                                                                                <  4:796751/149‑1 (MQ=255)
                                 aCCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                    >  3:330460/1‑126 (MQ=255)
                                 aCCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                    <  4:330460/126‑1 (MQ=255)
                                     aGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGt                                                                                                          >  2:266456/1‑148 (MQ=255)
                                       gTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTAc                                                                                                        >  1:783312/1‑148 (MQ=255)
                                       gTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTAc                                                                                                        >  2:808874/1‑148 (MQ=255)
                                         cATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                      >  1:626320/1‑148 (MQ=255)
                                             aaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGa                                                                                                  >  1:182987/1‑148 (MQ=255)
                                                   gATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTg                                                                                            <  1:266456/148‑1 (MQ=255)
                                                      cACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATaa                                                                                        <  4:239431/149‑1 (MQ=255)
                                                      cACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATa                                                                                         <  1:232377/148‑1 (MQ=255)
                                                      cACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATa                                                                                         >  1:819278/1‑148 (MQ=255)
                                                       aCGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATaa                                                                                        <  2:626320/148‑1 (MQ=255)
                                                       aCGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAAc                                                                                       <  3:437112/149‑1 (MQ=255)
                                                       aCGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAAc                                                                                       >  4:141854/1‑149 (MQ=255)
                                                        cGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAAc                                                                                       <  1:501954/148‑1 (MQ=255)
                                                          ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTac                                                                                    <  3:688287/149‑1 (MQ=255)
                                                               cTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCa                                                                                >  1:11902/1‑148 (MQ=255)
                                                                   ccGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCa                                                                            <  2:819278/148‑1 (MQ=255)
                                                                       aaTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTc                                                                        >  1:209382/1‑148 (MQ=255)
                                                                        aTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCg                                                                       <  1:839702/148‑1 (MQ=255)
                                                                             gCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                  >  2:625015/1‑148 (MQ=255)
                                                                              cccGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGcc                                                                 >  2:857210/1‑148 (MQ=255)
                                                                              cccGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCa                                                                >  4:547660/1‑149 (MQ=255)
                                                                               ccGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCa                                                                <  2:318771/148‑1 (MQ=255)
                                                                               ccGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCa                                                                <  1:625015/148‑1 (MQ=255)
                                                                                  tAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAggg                                                             >  2:906691/1‑148 (MQ=255)
                                                                                    ggCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCt                                                           >  2:46618/1‑148 (MQ=255)
                                                                                      cGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTc                                                         <  1:906691/148‑1 (MQ=255)
                                                                                      cGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTc                                                         <  2:126540/148‑1 (MQ=255)
                                                                                       gATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                      >  3:500317/1‑102 (MQ=255)
                                                                                       gATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                      <  4:500317/102‑1 (MQ=255)
                                                                                         tAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGtt                                                     >  4:148981/1‑149 (MQ=255)
                                                                                             cGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGcc                                                  <  1:100768/148‑1 (MQ=255)
                                                                                                        gtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAg                                                               <  3:731487/124‑1 (MQ=255)
                                                                                                        gtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAg                                                               >  4:731487/1‑124 (MQ=255)
                                                                                                           cGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGc                                    <  1:693185/148‑1 (MQ=255)
                                                                                                           cGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGc                                    <  2:485182/148‑1 (MQ=255)
                                                                                                            gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                  <  4:830490/149‑1 (MQ=255)
                                                                                                            gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                  <  4:851961/149‑1 (MQ=255)
                                                                                                             ttAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATaa                                                                                        <  4:44742/94‑1 (MQ=255)
                                                                                                             ttAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATaa                                                                                        >  3:44742/1‑94 (MQ=255)
                                                                                                             ttAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCttt                                 >  3:354141/1‑149 (MQ=255)
                                                                                                             ttAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                  <  2:812031/148‑1 (MQ=255)
                                                                                                               aagtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTg                                <  1:144684/148‑1 (MQ=255)
                                                                                                                agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACTCCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                              <  4:46705/149‑1 (MQ=255)
                                                                                                                agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGa                               <  1:239497/148‑1 (MQ=255)
                                                                                                                agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGa                               <  1:273407/148‑1 (MQ=255)
                                                                                                                 gtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                              <  3:177661/148‑1 (MQ=255)
                                                                                                                 gtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                              >  4:177661/1‑148 (MQ=255)
                                                                                                                  taagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTa                                                                                     >  2:315763/1‑92 (MQ=255)
                                                                                                                  taagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTa                                                                                     <  1:315763/92‑1 (MQ=255)
                                                                                                                   aagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGa                            >  1:577927/1‑148 (MQ=255)
                                                                                                                   aagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGAt                           <  3:296604/149‑1 (MQ=255)
                                                                                                                    agCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTa                                                                                     <  2:46544/90‑1 (MQ=255)
                                                                                                                    agCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTa                                                                                     >  1:46544/1‑90 (MQ=255)
                                                                                                                      cgcAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATa                                                                                         >  2:270481/1‑84 (MQ=255)
                                                                                                                      cgcAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATa                                                                                         <  1:270481/84‑1 (MQ=255)
                                                                                                                        cAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTcc                       >  2:115933/1‑148 (MQ=255)
                                                                                                                        cAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCg                      >  3:524939/1‑149 (MQ=255)
                                                                                                                          ggACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGc                     >  1:488754/1‑148 (MQ=255)
                                                                                                                              gTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGaaa                <  3:272086/149‑1 (MQ=255)
                                                                                                                               tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTc                        <  2:791106/140‑1 (MQ=255)
                                                                                                                               tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTc                        >  1:791106/1‑140 (MQ=255)
                                                                                                                               tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTc                        <  2:58232/140‑1 (MQ=255)
                                                                                                                               tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTc                        >  1:58232/1‑140 (MQ=255)
                                                                                                                               tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAc               >  3:77421/1‑149 (MQ=255)
                                                                                                                                  aGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcgc             >  2:204261/1‑148 (MQ=255)
                                                                                                                                   ggCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcgcg            >  2:550326/1‑148 (MQ=255)
                                                                                                                                       ccaaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCg       >  3:729120/5‑149 (MQ=255)
                                                                                                                                       cGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc        >  2:568791/1‑148 (MQ=255)
                                                                                                                                        gCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCg       <  2:714943/148‑1 (MQ=255)
                                                                                                                                        caaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                              <  3:413581/122‑1 (MQ=255)
                                                                                                                                        caaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                              >  4:413581/4‑125 (MQ=255)
                                                                                                                                        caaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg      >  3:492870/4‑149 (MQ=255)
                                                                                                                                         ccAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg      >  2:299046/1‑148 (MQ=255)
                                                                                                                                         aaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg      >  1:630110/3‑148 (MQ=255)
                                                                                                                                            gACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGCACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTc  <  4:476024/149‑1 (MQ=255)
                                                                                                                                             aCGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTc  <  1:327649/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTC  >  NC_000913/2658276‑2658564

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: