Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,657,426 G→A 6.2% L504L (CTG→TTG)  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,657,4260GA 6.2% 201.6 / 2.5 64L504L (CTG→TTG) hscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base A (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CAGCAGCCAGCGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCG  >  NC_000913/2657288‑2657572
                                                                                                                                          |                                                                                                                                                  
cagcagCCAGCGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGAc                                                                                                                                          <  1:166170/149‑1 (MQ=255)
       cagcGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGa                                                                                                                                   <  2:484496/149‑1 (MQ=255)
        agcGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGAc                                                                                                                                  <  2:1004855/149‑1 (MQ=255)
        agcGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGAc                                                                                                                                  <  1:170593/149‑1 (MQ=255)
         gcgcgcCGTGCAGACGTTCCAGCACACGCGCCGGTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGAcc                                                                                                                                 <  2:270336/149‑1 (MQ=255)
          cgcgcCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCt                                                                                                                                <  2:730292/149‑1 (MQ=255)
           gcgcCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTg                                                                                                                               <  2:1276265/149‑1 (MQ=255)
           gcgcCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTg                                                                                                                               >  1:1337408/1‑149 (MQ=255)
             gcCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAAACCGTCTTACGGTCTGAccgata                                                                                                                             <  2:699353/149‑5 (MQ=255)
              ccGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGGCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAAACCGTCTTACGGTCTGAccgatag                                                                                                                            <  1:1023669/149‑6 (MQ=255)
                  gCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGAc                                                                                                                        <  1:1026968/149‑1 (MQ=255)
                      aCTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCt                                                                                                                    <  2:1447540/149‑1 (MQ=255)
                          tCCAGCTCTCGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTGGTGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGGCTTTGATCATCTAAGCTATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAAc                                                                                                                <  1:354064/149‑1 (MQ=255)
                            cAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGc                                                                                                              <  2:1409127/149‑1 (MQ=255)
                              gCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCg                                                                                                            <  1:225755/149‑1 (MQ=255)
                                 ggcgcgCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTg                                                                                                         >  2:1326762/3‑149 (MQ=255)
                                    gcgcCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGAt                                                                                                      <  2:720555/149‑1 (MQ=255)
                                      gccgcTTGTACTTTTTGTTCTTCCAGCATTCGGGCTTTTACGGCCTGCTCGGCATAGCTCAGTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGAttt                                                                                                    <  1:162660/149‑1 (MQ=255)
                                              tACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCg                                                                                            >  2:575024/1‑149 (MQ=255)
                                                cTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGcc                                                                                          >  1:1054362/1‑149 (MQ=255)
                                                 tttttGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCg                                                                                         >  1:975291/1‑149 (MQ=255)
                                                      gTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCTCTATCGGTCAGACCGTAAGACGGTTTGACCGGAATAGACGCCTCAACGCCGGTCGATTTCTCCATCGCCGTCAcc                                                                                    >  2:1533800/1‑148 (MQ=255)
                                                          tGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCCCGCTCa                                                                                >  1:1191280/1‑149 (MQ=255)
                                                             cAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACGGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCaaaa                                                                             <  1:898347/149‑1 (MQ=255)
                                                             cAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCaaaa                                                                             >  2:658934/1‑149 (MQ=255)
                                                                  ttCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCg                                                                        >  1:542549/1‑149 (MQ=255)
                                                                     gggCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCg                                                                     >  1:1262168/1‑149 (MQ=255)
                                                                        cTTTTACGTCCTGCGCGGCATAGCTCATTGAGTCTTTTATCATCGAAGAGATTTCGCTATCGGTCAGACCGTAAGACGGGTTGACCTGAATAGACGCCTCAAAGCCGGTGGATTTCTCCCTCGCCTTCACGCACAAAAGACAGTAGGCa                                                                  >  2:856158/1‑149 (MQ=255)
                                                                         ttttACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCAt                                                                 >  1:879512/1‑149 (MQ=255)
                                                                            tACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGa                                                              <  1:715082/149‑1 (MQ=255)
                                                                                  cTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGa                                                        >  1:1271972/1‑149 (MQ=255)
                                                                                        ggCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGCCCGTAAGACGGTATGACCTGAATAGACGCCTCAACGGCGGTGGATTTCTACATAGCCGTCACGCTAAAAATCCAGTCGGCATCCCCATGGAAAGTCa                                                  >  2:1155485/1‑149 (MQ=255)
                                                                                          cATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCAcg                                                >  1:70190/1‑149 (MQ=255)
                                                                                           aTAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCAcgc                                               >  2:1320386/1‑149 (MQ=255)
                                                                                              gCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGaa                                            >  2:731429/1‑149 (MQ=255)
                                                                                              gCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGaa                                            >  1:345119/1‑149 (MQ=255)
                                                                                                tCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTAACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAata                                          >  1:1534817/1‑149 (MQ=255)
                                                                                                     gAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATgcgc                                     >  2:196819/1‑149 (MQ=255)
                                                                                                          ttttATCATCGAAGCGATGTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTGCTCCATCGCCGTCACGCTCAACAGACCGCCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGc                                <  2:1021266/149‑1 (MQ=255)
                                                                                                          tttGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGc                                <  2:993413/149‑1 (MQ=255)
                                                                                                           ttGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGcc                               <  1:1292191/149‑1 (MQ=255)
                                                                                                            tGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCa                              <  2:975291/149‑1 (MQ=255)
                                                                                                             gATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAg                             >  2:1523452/1‑149 (MQ=255)
                                                                                                              atcatcGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTc                                                                                       <  1:479168/90‑1 (MQ=255)
                                                                                                              atcatcGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTc                                                                                       >  2:479168/1‑90 (MQ=255)
                                                                                                                catcGAAGCGATTTCGCTATCGGTCAGACGGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATATCTCCATCGACGTCACGCTCAAAAGACCGACCGCATCGCCCTGGAACGTCAGGCGAATAAGCGCACCGCAAGCCg                          >  2:278481/1‑149 (MQ=255)
                                                                                                                 atcGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCgg                         >  2:948354/1‑149 (MQ=255)
                                                                                                                 atcGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCgg                         <  1:1326762/149‑1 (MQ=255)
                                                                                                                  tcGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGt                        >  2:1230002/1‑149 (MQ=255)
                                                                                                                    gAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAg                      >  1:1147981/1‑149 (MQ=255)
                                                                                                                     aaGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGCCCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAgc                     >  1:705138/1‑149 (MQ=255)
                                                                                                                     aaGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAgc                     >  1:1456142/1‑149 (MQ=255)
                                                                                                                       gCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAgcgc                   >  1:1571841/1‑149 (MQ=255)
                                                                                                                        cGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGcc                  <  2:1337408/149‑1 (MQ=255)
                                                                                                                         gATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCg                 >  1:1029718/1‑149 (MQ=255)
                                                                                                                          aTTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCgg                >  2:429695/1‑149 (MQ=255)
                                                                                                                          aTTTCGCAATCGGTCAGACCGTAAGACTGTTTGACCGGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCTCGCGAATATGCGCACCGCCAGCCGGTAGCGCCgg                <  1:1180710/149‑1 (MQ=255)
                                                                                                                            ttCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGaa              >  2:1313911/1‑149 (MQ=255)
                                                                                                                               tattcaggTCAAACCGTCTTACGGTCTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATAc           >  2:781908/7‑149 (MQ=255)
                                                                                                                                 tATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCggccgccac         >  2:1367674/1‑142 (MQ=255)
                                                                                                                                   tCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACg       >  1:1166473/1‑149 (MQ=255)
                                                                                                                                     ggTCAAACCGTCTTACGGTCTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCACAAAAGACCGTCGGCCTCGACCTGGAACGTCACGCGAATAAGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGcc     >  1:1051252/1‑148 (MQ=255)
                                                                                                                                      gTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAg    <  1:196819/149‑1 (MQ=255)
                                                                                                                                        cAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAgcg  <  2:1507504/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CAGCAGCCAGCGCGCCGTGCAGACTTTCCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCG  >  NC_000913/2657288‑2657572

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