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	98.8%	282.2 / ‑4.0	83	noncoding (64/195 nt)	isrC	small RNA IsrC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/54);  total (28/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

CAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTC  >  NC_000913/2071320‑2071441                                                             |                                                              cAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACAt                                                           <  1:3807356/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:1005549/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2713633/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2706849/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2649530/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2560020/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2998882/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:3197475/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:2033691/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:1820009/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:1792520/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:698376/65‑1 (MQ=255)   atatCTATTTTATCGATCGTTTACATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGc                                                         <  1:565796/65‑1 (MQ=255)     atCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcac                                                       <  1:3415897/65‑1 (MQ=255)      tCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcaca                                                      >  1:2775617/1‑65 (MQ=255)      tCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcaca                                                      <  1:3150826/65‑1 (MQ=255)      tCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcaca                                                      <  1:2722906/65‑1 (MQ=255)      tCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGcaca                                                      >  1:499918/1‑65 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATATCCTTTGTCAGTGACATGCACAGa                                                    <  1:2394256/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2170901/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2881833/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:65949/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:6922/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2508867/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:251733/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2969075/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2737492/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:2883633/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:140631/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:183595/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:1834273/65‑1 (MQ=255)        tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGa                                                    <  1:455931/65‑1 (MQ=255)          ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATa                                                  >  1:370950/1‑65 (MQ=255)          ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATa                                                  >  1:1756418/1‑65 (MQ=255)           tttATCGATCGTTTATATCGATCGATAAGCTAGTAATAACCTTTGTCAGTGACATGCAGAGATAc                                                 >  1:2855539/1‑65 (MQ=255)              atcgatcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTa                                              >  1:3967613/1‑65 (MQ=255)              atcgatcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTa                                              >  1:1253451/1‑65 (MQ=255)                cgatcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACa                                            <  1:1890211/65‑1 (MQ=255)                       ttATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGAc                                     <  1:3065056/65‑1 (MQ=255)                       ttATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGAc                                     <  1:961742/65‑1 (MQ=255)                        tataGCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACa                                    <  1:2671772/65‑1 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:1057275/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:2168295/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:2929030/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:1444644/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:1786586/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:341473/1‑65 (MQ=255)                           atcgatcgATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTc                                 >  1:3230243/1‑65 (MQ=255)                              ggtcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:994977/63‑1 (MQ=255)                              gatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:3695020/65‑1 (MQ=255)                              gatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:2750106/65‑1 (MQ=255)                              gatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:3653916/65‑1 (MQ=255)                              gatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:2296985/65‑1 (MQ=255)                              gatcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAgg                              <  1:1747611/65‑1 (MQ=255)                                tcgatAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGa                            <  1:63536/65‑1 (MQ=255)                                tcgatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGa                            <  1:2416038/65‑1 (MQ=255)                                    tAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGaacaa                        >  1:1068253/1‑65 (MQ=255)                                      aGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAacaaca                      <  1:1015485/65‑1 (MQ=255)                                      aGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAacaaca                      <  1:2862415/65‑1 (MQ=255)                                       gCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAg                     >  1:944069/1‑65 (MQ=255)                                        cTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGa                    <  1:739689/65‑1 (MQ=255)                                        cTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGa                    <  1:1246328/65‑1 (MQ=255)                                         taataataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGaa                   <  1:1267535/65‑1 (MQ=255)                                           ataataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAAcc                 <  1:1998096/65‑1 (MQ=255)                                              cagaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              <  1:3085040/62‑1 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcacg              >  1:1239337/1‑64 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              <  1:1344462/65‑1 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:3794128/1‑65 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:972528/1‑65 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:1386834/1‑65 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:2923084/1‑65 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              <  1:3657511/65‑1 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              <  1:2938579/65‑1 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:3834730/1‑65 (MQ=255)                                              ataaCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca              >  1:2290571/1‑65 (MQ=255)                                                 tccTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt           <  1:2115057/64‑1 (MQ=255)                                                 aCCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt           >  1:705629/1‑65 (MQ=255)                                                  ccTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTc          >  1:613223/1‑65 (MQ=255)                                                  ccTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTc          >  1:206814/1‑65 (MQ=255)                                                  ccTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTc          >  1:273266/1‑65 (MQ=255)                                                   cTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCa         >  1:3637284/1‑65 (MQ=255)                                                        tCAGTGACATGCACAGATACGTACAGAAAGAAATTCAGGGAACAACAGAACCACAATTCAGAAAc    >  1:3204233/1‑65 (MQ=255)                                                         cAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACt   <  1:783988/65‑1 (MQ=255)                                                         cAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACt   <  1:1098998/65‑1 (MQ=255)                                                          aGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTc  <  1:1305061/65‑1 (MQ=255)                                                          aGTGACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTc  <  1:1480414/65‑1 (MQ=255)                                                          aGTGACATGCACAGATACGTACAGAAACACATTCAGGGAACAACAGAACCACAATTCAGAAACTc  <  1:2252231/65‑1 (MQ=255)                                                             |                                                              CAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTC  >  NC_000913/2071320‑2071441

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