Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,658,422 C→G 5.2% 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 5.2% 380.5 / 7.0 110A172P (GCG→CCG) hscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base C (48/56);  new base G (2/4);  total (50/60)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 116 total aligned reads displayed.

CTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCC  >  NC_000913/2658274‑2658571
                                                                                                                                                    |                                                                                                                                                     
cTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                                                                                                                                      <  1:535352/150‑1 (MQ=255)
cTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGGCATAAACGGCGATCACGCCGTCCTGACCGGAATCCAGCCCGTAGGCGAGAGCCGCAGCGGGCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                                                                                                                                      <  1:254429/150‑1 (MQ=255)
  tAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAATTAAGCGCAGGACGTGAAGGCCCGCCAgccgcgccgc                                                                                                                                                    >  2:326107/1‑150 (MQ=255)
    aGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGt                                                                                                                                                  <  1:1546739/150‑1 (MQ=255)
      cgcAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCt                                                                                                                                                <  2:1103896/150‑1 (MQ=255)
        cAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCttt                                                                                                                                              <  2:1528305/150‑1 (MQ=255)
         aGAATGGAAATATCAAACCTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGCCCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtggc                                                                                                                                             <  2:140058/150‑4 (MQ=255)
          gAATGGAAATATCAAACGTCCCGCCCCCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGGAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTgg                                                                                                                                            <  2:111380/150‑1 (MQ=255)
           acttgAAATATCAAACCTCCCGCCCCCGAGGTCATAATCGGCGATCACGCCTTCCTGTCCGGACTCCCGCCGGTTGGCGTTAGCCGCAGCGGTCGGTTCGGTATGTAAGCGCATGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                           <  1:100104/146‑1 (MQ=255)
           aaTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                           <  1:295303/150‑1 (MQ=255)
            aTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTg                                                                                                                                          <  1:833537/150‑1 (MQ=255)
            aTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTg                                                                                                                                          <  2:440662/150‑1 (MQ=255)
               gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGATTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       >  1:1263973/1‑150 (MQ=255)
               gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       <  1:47107/150‑1 (MQ=255)
               gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       >  1:341758/1‑150 (MQ=255)
               gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       >  1:1065064/1‑150 (MQ=255)
               gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGACCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       >  2:1433220/1‑150 (MQ=255)
                  atatCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                                                                                                                                    >  1:686162/1‑150 (MQ=255)
                  atatCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                                                                                                                                    >  1:434720/1‑150 (MQ=255)
                      cAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCt                                                                                                                                >  2:889702/1‑150 (MQ=255)
                       aaaCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg                                                                                                                               >  2:624988/1‑150 (MQ=255)
                       aaaCGTCCCGCCACAGAGGTCATAAACGGCGATCACGCCTTACTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg                                                                                                                               >  1:109090/1‑150 (MQ=255)
                        aaCGTCCCGCCACCGCGGGCATACACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg                                                                                                                              <  1:1084395/150‑1 (MQ=255)
                         aCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg                                                                                                                             <  2:2203/150‑1 (MQ=255)
                         aCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg                                                                                                                             <  1:1464352/150‑1 (MQ=255)
                         aCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg                                                                                                                             <  1:1372944/150‑1 (MQ=255)
                          cGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc                                                                                                                            >  2:530435/1‑150 (MQ=255)
                             cccGCCACCGAGGTCATAAACTGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCTGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc                                                                                                                         <  2:1072046/150‑1 (MQ=255)
                                 ccACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCa                                                                                                                     >  1:156369/1‑150 (MQ=255)
                                 ccACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCa                                                                                                                     <  1:742048/150‑1 (MQ=255)
                                    ccGAGGGCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAg                                                                                                                  <  2:334553/150‑1 (MQ=255)
                                      gAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTa                                                                                                                <  1:1315777/150‑1 (MQ=255)
                                      gAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTa                                                                                                                >  2:1248340/1‑150 (MQ=255)
                                         gTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                             >  2:690136/1‑150 (MQ=255)
                                         gTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                             >  2:62880/1‑150 (MQ=255)
                                           cATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCg                                                                                                           >  1:1663187/1‑150 (MQ=255)
                                             tAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGa                                                                                                         >  1:1542983/1‑150 (MQ=255)
                                               accGGCGATCACGCCTTCCTGCCCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGGTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAAc                                                                                                       <  1:1372859/148‑1 (MQ=255)
                                               aaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCTATAGCCGCAGCGGTCGGTTCGTTAAGTACGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCTAAGTACGCCGGAAc                                                                                                       <  1:13845/150‑1 (MQ=255)
                                                aCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACg                                                                                                      >  1:47461/1‑150 (MQ=255)
                                                   gCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTg                                                                                                   <  2:1025789/150‑1 (MQ=255)
                                                    cGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                  <  1:676979/150‑1 (MQ=255)
                                                    cGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGa                                                                                                  <  1:820940/150‑1 (MQ=255)
                                                     gATCACTCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAAc                                                                                                       >  2:1677765/1‑144 (MQ=255)
                                                     gATCACTCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAAc                                                                                                       <  1:1677765/144‑1 (MQ=255)
                                                      aTCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACATGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATa                                                                                                <  2:1307154/150‑1 (MQ=255)
                                                       tCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATaa                                                                                               >  1:1406046/1‑150 (MQ=255)
                                                          cGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGCGATAACTa                                                                                            >  2:1291832/1‑150 (MQ=255)
                                                           gCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTac                                                                                           <  2:1514853/150‑1 (MQ=255)
                                                            ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                          >  2:825743/1‑150 (MQ=255)
                                                            ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                          >  2:1177605/1‑150 (MQ=255)
                                                            ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                          >  2:1148406/1‑150 (MQ=255)
                                                             cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTacac                                                                                         <  2:813145/150‑1 (MQ=255)
                                                                 cTGACCGGAATCCCGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCGTTGGTGCCCTGCCGCGGGGCATCGTCAAAGTACGCCGGCACGGTGATAACTACCCCATc                                                                                     >  1:807309/1‑150 (MQ=255)
                                                                  tGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATcc                                                                                    <  1:624988/150‑1 (MQ=255)
                                                                  tGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATcc                                                                                    <  2:525805/150‑1 (MQ=255)
                                                                   gACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATccc                                                                                   >  2:1298168/1‑149 (MQ=255)
                                                                       ggAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAg                                                                                  <  1:1553164/147‑1 (MQ=255)
                                                                       ggAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAg                                                                                  >  2:1553164/1‑147 (MQ=255)
                                                                       ggAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTc                                                                               >  2:80652/1‑150 (MQ=255)
                                                                          aTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGcc                                                                            >  2:505986/1‑150 (MQ=255)
                                                                            ccAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTTGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTg                                                                          >  2:860953/1‑150 (MQ=255)
                                                                             cAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                         >  2:1537196/1‑150 (MQ=255)
                                                                             cAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                         <  2:434720/150‑1 (MQ=255)
                                                                             cAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                         >  2:492838/1‑150 (MQ=255)
                                                                                cccGTAGGCGATAGCCGAAGCGGGCTGTTAGTTAAGTAAGCGCAGGACGTGAAGTCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAg                                                                      <  2:1055232/150‑1 (MQ=255)
                                                                                 ccGTAGGCGATAGCCGCAGCGGGCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAgg                                                                     <  2:482326/150‑1 (MQ=255)
                                                                                      ggCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCtgc                                                                >  1:357744/1‑148 (MQ=255)
                                                                                       gCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCa                                                               <  2:1031434/150‑1 (MQ=255)
                                                                                           tAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTg                                                           <  1:935098/150‑1 (MQ=255)
                                                                                           tAGCCGCAGCGGGCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTg                                                           <  2:1128990/150‑1 (MQ=255)
                                                                                             gCCGCAGCGGTCGGGGAGGTAAGTACGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCTTTGGTGCCCTGACGCTGTGCATCGTCAAAGTACGCCGTAACTGTGATAACTACACCATCCAGCTCGCCTTCCAGTGCTTCAGTTGcc                                                         <  2:292491/150‑1 (MQ=255)
                                                                                              ccGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGccc                                                        <  1:530435/150‑1 (MQ=255)
                                                                                                   gCGGTCGGTTCGTTAAGTAAGCGCCGGACGTGAAGGCCCGCCATACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGcgcc                                                   <  1:1078363/150‑1 (MQ=255)
                                                                                                    cggtcggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGcgccg                                                  <  2:1381442/150‑1 (MQ=255)
                                                                                                     ggtcggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGcgccgc                                                 >  1:1409380/1‑150 (MQ=255)
                                                                                                         ggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCGTCAGTTGCCCGCGCCGCCAGt                                             >  1:1248154/1‑150 (MQ=255)
                                                                                                          gtTCGTTAAGTAAGCGCCGGCCGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTg                                            <  1:326107/150‑1 (MQ=255)
                                                                                                            tCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACCCCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCt                                          >  1:1135093/1‑150 (MQ=255)
                                                                                                             cGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                         >  1:213427/1‑150 (MQ=255)
                                                                                                              gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCGGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCttt                                        >  2:630049/1‑150 (MQ=255)
                                                                                                              gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCttt                                        <  1:690136/150‑1 (MQ=255)
                                                                                                              gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCttt                                        <  1:62880/150‑1 (MQ=255)
                                                                                                                taagtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTGTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCAACCACCTCGCCAGCCAGGGCTTCAGTTGCCCGCGCCGCCAGGGcgatga                                      >  1:843448/1‑145 (MQ=255)
                                                                                                                taagtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGa                                      <  2:1406046/150‑1 (MQ=255)
                                                                                                                 aagtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGCGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAg                                     <  1:1291832/150‑1 (MQ=255)
                                                                                                                   gtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGa                                   <  1:1424766/150‑1 (MQ=255)
                                                                                                                     aagCGCAGGACGTGAAGGCCCGCCATACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATg                                 <  1:907768/150‑1 (MQ=255)
                                                                                                                      agCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGt                                >  1:237693/1‑150 (MQ=255)
                                                                                                                       gcgcAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTc                               <  1:1433220/150‑1 (MQ=255)
                                                                                                                            ggACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAg                          <  1:1016118/150‑1 (MQ=255)
                                                                                                                            ggACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAg                          <  1:80652/150‑1 (MQ=255)
                                                                                                                                gTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTGTGAGGATGTCCGCAGAAAc                      >  2:1439355/1‑150 (MQ=255)
                                                                                                                                 tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcg                     >  1:922780/1‑150 (MQ=255)
                                                                                                                                       cccGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc               >  2:192459/1‑150 (MQ=255)
                                                                                                                                         ccaaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg             >  2:565966/5‑150 (MQ=255)
                                                                                                                                          caaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCggg            >  1:925697/4‑150 (MQ=255)
                                                                                                                                                   gccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCagccgc   >  2:488495/1‑150 (MQ=255)
                                                                                                                                                    ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGcc  <  2:172859/150‑1 (MQ=255)
                                                                                                                                                    ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGcc  >  1:502865/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCC  >  NC_000913/2658274‑2658571

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