Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,658,447 G→C 5.7% A163A (GCC→GCG) ‡ 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,4470GC 5.7% 222.0 / 2.7 66A163A (GCC→GCG) ‡hscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (2/2);  total (35/31)
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.88e-01

ACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCA  >  NC_000913/2658299‑2658595
                                                                                                                                                    |                                                                                                                                                    
aCGTCCCGCCACCGAGGTCATAACCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg                                                                                                                                                      <  1:1493416/149‑1 (MQ=255)
  gTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTACCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCTCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc                                                                                                                                                    >  1:613007/1‑149 (MQ=255)
    cccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCAt                                                                                                                                                  <  2:291529/149‑1 (MQ=255)
    cccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCAt                                                                                                                                                  >  1:51655/1‑149 (MQ=255)
     ccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCCGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc                                                                                                                                                 <  1:1581886/149‑1 (MQ=255)
     ccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc                                                                                                                                                 <  1:388795/149‑1 (MQ=255)
          aCCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCaa                                                                                                                                            <  1:1041537/149‑1 (MQ=255)
           ccGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCaaa                                                                                                                                           <  2:876226/149‑1 (MQ=255)
            cGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAg                                                                                                                                          >  2:854450/1‑149 (MQ=255)
               ggTCATAACCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGAGAGCCGCAGCGGGCGGTTCGTTAAGTAAGCGCAGGACGGGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTAc                                                                                                                                       <  1:449325/149‑1 (MQ=255)
               ggTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTAc                                                                                                                                       <  2:664299/149‑1 (MQ=255)
                 tCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                                                     >  1:1353394/1‑149 (MQ=255)
                     aaaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAACCGCAGGACCTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGa                                                                                                                                 >  2:1230099/1‑149 (MQ=255)
                        cGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCCTCGTCAAAGTACGCCGGAACg                                                                                                                              >  1:1478043/1‑149 (MQ=255)
                        cGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACg                                                                                                                              >  1:696380/1‑149 (MQ=255)
                            gcgcACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCCGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                                          <  1:1154019/146‑1 (MQ=255)
                            gATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTACGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                                          <  2:1127217/149‑1 (MQ=255)
                              tCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGCACGGGGATa                                                                                                                        >  1:1128117/1‑149 (MQ=255)
                                    cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCTCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                                                  <  1:1178867/149‑1 (MQ=255)
                                    cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                                                  <  1:1127648/149‑1 (MQ=255)
                                          gACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGGCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATcc                                                                                                            <  2:1157918/149‑1 (MQ=255)
                                            ccGGAATCCAGCCAGTAGGCGATAGCCGCAGCGGTCGGTTCGAGAAGTAAGCGCAGGACGTGAAGGCCCGCCAGCCGCGCCGCGTCTTTGGTGCCCTGCCGCTGGGCATCGTCAAAGTACGCCGGAACGGTGCTACCTCCCCCATCCAg                                                                                                          >  1:263162/1‑149 (MQ=255)
                                             cGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGc                                                                                                         >  2:1416168/1‑149 (MQ=255)
                                             cGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGc                                                                                                         >  2:490098/1‑149 (MQ=255)
                                                  tCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGcc                                                                                                    >  1:1127295/1‑149 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                                          >  2:203904/1‑124 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                                          <  1:203904/124‑1 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                                                 >  2:769524/1‑149 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                                                 >  2:838326/1‑149 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                        >  1:1302505/1‑105 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                        <  2:1296147/110‑6 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                        <  2:1302505/110‑6 (MQ=255)
                                                     aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                        >  1:1296147/1‑105 (MQ=255)
                                                         cGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGACCTGACGCTGGGAATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAgg                                                                                             >  1:1211761/1‑149 (MQ=255)
                                                           gaGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGc                                                                                           <  1:130171/148‑1 (MQ=255)
                                                             ggCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCtt                                                                                         >  1:1156548/1‑149 (MQ=255)
                                                              gCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTc                                                                                        >  1:1453020/1‑149 (MQ=255)
                                                                gATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAg                                                                                      <  1:41296/149‑1 (MQ=255)
                                                                gATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGCGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAg                                                                                      <  1:1537119/149‑1 (MQ=255)
                                                                     ccGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACAACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGcc                                                                                 >  2:1107391/1‑149 (MQ=255)
                                                                       gCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCcg                                                                               >  2:611006/1‑149 (MQ=255)
                                                                              tcggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCgccgcc                                                                        >  2:634788/1‑149 (MQ=255)
                                                                                  ttCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTg                                                                    <  1:101086/149‑1 (MQ=255)
                                                                                         agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGGCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                                                             >  1:593408/1‑149 (MQ=255)
                                                                                           taagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGa                                                           >  2:782086/1‑149 (MQ=255)
                                                                                            atgcgcAGGAGGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGTCGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACTCCTTCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGAt                                                          <  2:613007/147‑1 (MQ=255)
                                                                                                 cAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCg                                                     <  2:51655/149‑1 (MQ=255)
                                                                                                         gAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcg                                             >  1:748448/1‑149 (MQ=255)
                                                                                                         gAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcg                                             >  2:353162/1‑149 (MQ=255)
                                                                                                               ccGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc                                       >  1:120080/1‑149 (MQ=255)
                                                                                                               accaaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc                                       >  1:641752/6‑149 (MQ=255)
                                                                                                                  ccAGACGCGCCGCGTCTTTGGTGCCGTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCggg                                    <  2:980089/149‑1 (MQ=255)
                                                                                                                         cgccgcGTCGTTGGTGCCCTGACGCTGGGCATCGTCTAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTcagca                             <  2:1449658/149‑1 (MQ=255)
                                                                                                                           ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCagcagc                           <  1:1230099/149‑1 (MQ=255)
                                                                                                                           ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCagcagc                           <  2:949914/149‑1 (MQ=255)
                                                                                                                              cgTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGcccc                        >  1:1171769/1‑149 (MQ=255)
                                                                                                                                  tttGGTGCCCTGACGCTGGGCAGCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCccgccg                    <  1:78359/149‑1 (MQ=255)
                                                                                                                                       tGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGATTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGtt               >  2:1195103/1‑149 (MQ=255)
                                                                                                                                        gCCCTGACGCTGGGCATCGTCATAGTACGCCGGAACGGTGATAACTACATCATCCTGGTCGCTTGCGAGGGCTTCATTTGCCAGCTCAGCCAGTGCTTTGAGGATGTCCGCAGAAACGCTTATCGGTTTGAGCTGCCCCGCCGCCGttt              <  2:1211761/149‑1 (MQ=255)
                                                                                                                                        gCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGttt              <  1:244644/149‑1 (MQ=255)
                                                                                                                                         cccTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTc             >  1:315877/1‑149 (MQ=255)
                                                                                                                                          ccTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCCCCGGGTTCAGCAGCCCCGCCGCCGTTTCa            >  1:1538828/1‑149 (MQ=255)
                                                                                                                                          ccTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg                                     >  1:693761/1‑124 (MQ=255)
                                                                                                                                          ccTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg                                     <  2:693761/124‑1 (MQ=255)
                                                                                                                                            tGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACAATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAat          >  1:1579486/1‑149 (MQ=255)
                                                                                                                                             gACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAatc         <  1:611006/149‑1 (MQ=255)
                                                                                                                                                  tGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCgg    <  1:448123/149‑1 (MQ=255)
                                                                                                                                                    ggCATCGTCAAAGTCCGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCa  <  1:854450/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCA  >  NC_000913/2658299‑2658595

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