Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,071,380	A→G	noncoding (64/195 nt)	isrC →	small RNA IsrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,071,380	0	A	G	100.0%	186.9 / NA	54	noncoding (64/195 nt)	isrC	small RNA IsrC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/39);  total (15/39)

ACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCC  >  NC_000913/2071315‑2071442                                                                  |                                                               aCGATCAATATCTATTTTATCGATCGTTTGTATCGATCGATAAGCTAATAATAACCTTTGTCAGTGAc                                                              >  1:740432/1‑68 (MQ=255) aCGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGAc                                                              >  1:1169648/1‑68 (MQ=255)   gATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACAt                                                            <  1:1282633/68‑1 (MQ=255)   gATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACAt                                                            <  1:1399430/68‑1 (MQ=255)   gATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACAt                                                            <  1:2773528/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:98303/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:773038/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:33576/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:2958701/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:1705052/68‑1 (MQ=255)     tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                          <  1:1543146/68‑1 (MQ=255)        atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcaca                                                       <  1:1824098/68‑1 (MQ=255)        atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCNTTGTCAGTGACATGcaca                                                       <  1:2708970/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:1562081/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:1239795/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:741447/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:1241921/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:2880498/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:2242783/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:1800674/68‑1 (MQ=255)          atCTATTTTATCGATCGTTTATATCGATCGAGAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                     <  1:3125278/68‑1 (MQ=255)               ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGTTACGt                                                <  1:1610583/68‑1 (MQ=255)               ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGt                                                >  1:1891130/1‑68 (MQ=255)               ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGt                                                >  1:546732/1‑68 (MQ=255)                tttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTa                                               >  1:803283/1‑68 (MQ=255)                  tATCGATCGTTTCTATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACa                                             >  1:784458/1‑68 (MQ=255)                   atcgatcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAg                                            >  1:75882/1‑68 (MQ=255)                        tcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGa                                       <  1:1141331/68‑1 (MQ=255)                        tcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGa                                       <  1:1427531/68‑1 (MQ=255)                         cgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGAc                                      <  1:3415188/68‑1 (MQ=255)                           tttATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACAt                                    <  1:2548646/68‑1 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               >  1:3341623/1‑68 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               <  1:20471/68‑1 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               >  1:3444979/1‑68 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               >  1:3459425/1‑68 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               >  1:570016/1‑68 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               <  1:88209/68‑1 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               >  1:3180119/1‑68 (MQ=255)                                atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                               <  1:2048402/68‑1 (MQ=255)                                  cgatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGa                             <  1:2000399/68‑1 (MQ=255)                                  cgatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGa                             <  1:3151600/68‑1 (MQ=255)                                     tcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGaaca                          <  1:3265760/68‑1 (MQ=255)                                       gatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGaacaac                        <  1:1270317/68‑1 (MQ=255)                                        atAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAacaaca                       <  1:1412108/68‑1 (MQ=255)                                        atAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAacaaca                       <  1:2657462/68‑1 (MQ=255)                                             cTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAAGAGAAcc                  <  1:1605794/68‑1 (MQ=255)                                                ataataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca               <  1:2685583/68‑1 (MQ=255)                                                   ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt            >  1:1454372/1‑68 (MQ=255)                                                      aCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAg         >  1:469193/1‑68 (MQ=255)                                                       ccTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGa        >  1:1734470/1‑68 (MQ=255)                                                         tttGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGaaa      <  1:2312592/68‑1 (MQ=255)                                                         tttGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGaaa      <  1:215979/68‑1 (MQ=255)                                                           tGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACt    <  1:3541362/68‑1 (MQ=255)                                                            gTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTc   <  1:2599394/68‑1 (MQ=255)                                                             tCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTcc  >  1:2306546/1‑68 (MQ=255)                                                                  |                                                               ACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCC  >  NC_000913/2071315‑2071442

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