Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,658,448 G→C 5.7% A163G (GCC→GGC) ‡ 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,4480GC 5.7% 220.8 / 2.4 66A163G (GCC→GGC) ‡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 = 1.00e+00

GTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCAG  >  NC_000913/2658301‑2658596
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                         <  2:1302505/110‑6 (MQ=255)
                                                   aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                         <  2:1296147/110‑6 (MQ=255)
                                                   aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                         >  1:1296147/1‑105 (MQ=255)
                                                   aGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGATGCCCAGCGTCagggca                                                                                                                                         >  1:1302505/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                            <  2:949914/149‑1 (MQ=255)
                                                                                                                         ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCagcagc                            <  1:1230099/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                                      <  2:693761/124‑1 (MQ=255)
                                                                                                                                        ccTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg                                      >  1:693761/1‑124 (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)
                                                                                                                                                   gCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCAg  >  2:495941/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCCCCGCCGCCGTTTCAATCATCGGCAG  >  NC_000913/2658301‑2658596

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