Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,944,114	(G)6→7	noncoding (411/2904 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,944,108	1	.	G	89.1%	130.2 / 13.2	55	noncoding (405/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (6/0);  new base G (20/29);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 7.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGG  >  NC_000913/3943975‑3944197                                                                                                                                        |                                                                                              agCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATCCAGGGTGACAGCGCCGTCCACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑g                                                                                     <  2:874820/140‑1 (MQ=255) agCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑g                                                                                     >  2:1355257/1‑140 (MQ=37) agCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑g                                                                                     >  1:2947/1‑140 (MQ=37) agCCTGAATCAGTATGTGTGATAGTTGAACCGCCTGGAAAGGCGCGCGATACAGGGAGACAGCCTCGTAAAAAAAAATGCTTATATTGAGAAAT‑GATGAGCAGGCCGGGAAACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggagtg                                                                                      >  2:1233489/1‑136 (MQ=2)   ccTGAATCAGTATGTGTTTTAGAGGAAGCGTCTGGAAAGTCGCGCGATCCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑gga                                                                                   <  2:403198/140‑3 (MQ=11)   ccTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑G‑GGG‑G‑GAc                                                                                   >  1:803077/1‑140 (MQ=38)      gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑ggacca                                                                                >  1:879163/1‑135 (MQ=255)      gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑ggacca                                                                                >  2:209883/1‑135 (MQ=255)      gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑g‑ggg‑g‑ggacca                                                                                >  1:1234258/1‑135 (MQ=255)           aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCAtcctc                                                                            <  2:1427001/140‑1 (MQ=255)           aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCAtcctc                                                                            <  1:343118/140‑1 (MQ=255)           aGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGAAACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCAtcctc                                                                            <  2:533969/140‑1 (MQ=255)             tATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          <  2:2947/140‑1 (MQ=255)             tATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          <  1:847715/140‑1 (MQ=255)             tATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          >  1:242886/1‑140 (MQ=255)             tATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          >  2:1269378/1‑140 (MQ=255)                gtgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAgg                                                                       <  1:1269378/140‑1 (MQ=21)                 tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGc                                                                      >  1:1179734/1‑140 (MQ=21)                   tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTTAATT‑T‑‑GTGGG‑G‑GATCATCCTCCAAGGCTa                                                                    >  1:474871/1‑140 (MQ=11)                   tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACA‑GTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTaa                                                                   >  1:805360/1‑140 (MQ=255)                   tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGAACACGTGGTATCGTGTATGAATA‑T‑‑GTGGG‑G‑CACCATCCTCCAAGGCTa                                                                    >  2:415660/1‑140 (MQ=11)                    gtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTaa                                                                   <  2:991209/140‑1 (MQ=21)                    gtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTaa                                                                   <  1:71934/140‑1 (MQ=21)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          <  1:1348974/130‑1 (MQ=25)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCa                                                                          >  2:1348974/1‑130 (MQ=25)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAgg                                                                       >  1:1391159/1‑133 (MQ=21)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAgg                                                                       <  2:1391159/133‑1 (MQ=21)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATa                                                                >  1:342219/1‑140 (MQ=21)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACAAGTGGTATCCTGTATGAATA‑TG‑G‑GGG‑G‑GACCGTCCTCCAAGGCTAAATa                                                                >  2:556044/1‑140 (MQ=11)                       aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACAAGTGGTATCATGTTTGATTA‑T‑‑G‑GGG‑GAGACCATCCTCCAAGGCTAAATa                                                                >  2:494988/1‑140 (MQ=11)                         tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTTTGAATATT‑‑G‑GGT‑G‑GGCCATCCACCAAGGCTAAATACt                                                              >  1:687116/1‑140 (MQ=11)                         tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTTTGAATA‑T‑‑G‑GGGCG‑GACCATCCTCCAAGGCTAAATACt                                                              >  1:675447/1‑140 (MQ=255)                           gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCAtcctcc                                                                           <  2:9892/125‑1 (MQ=25)                           gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGCTATCCTGTATGACTG‑TG‑G‑GGG‑G‑GACAATCCCCCAAGGCTAAATACTcc                                                            >  2:695937/1‑140 (MQ=11)                           gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAACTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCAtcctcc                                                                           >  1:9892/1‑125 (MQ=255)                              gCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCctga                                                         <  2:193804/140‑1 (MQ=21)                              gCGTCTGGAAAGGCGCGAGATAAAGGGTGACAGCCACGTACACAAAAACGCACATATTGTGAGCTCGAGGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCctga                                                         <  2:474871/140‑1 (MQ=11)                               cGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCctgac                                                        <  2:342219/140‑1 (MQ=21)                                 tCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCctgactg                                                      >  2:51450/1‑140 (MQ=21)                                   tGGAAAGGCGCGCGATATAGGGAGACAGCCCTGTAGACAAAAATGCACCTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCtgactgac                                                    <  2:570890/140‑1 (MQ=11)                                   tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCtgactgac                                                    <  2:1234258/140‑1 (MQ=21)                                    ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTTCACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                   <  1:1463411/140‑1 (MQ=255)                                    ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                   <  2:1368815/140‑1 (MQ=21)                                        aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGA‑GAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                              <  1:1300211/140‑1 (MQ=255)                                           cgcgcgATACAGGGTTACAGCCTCGCACAAAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                            <  1:413660/140‑1 (MQ=11)                                             cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TGGG‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGa                                           >  2:1166296/1‑140 (MQ=255)                                              gcgATACAGGGTGACCGCCCCGCACCCAACAATGCACTTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                         <  1:494988/140‑1 (MQ=11)                                              gcgATACAGGGTGACATCCCCGTACACAAAATTCCGCATTTTGTGAGCCCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                         <  2:688968/140‑1 (MQ=11)                                              gcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                         <  2:342729/140‑1 (MQ=21)                                                            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTAc                                   >  2:1342269/1‑132 (MQ=255)                                                            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTAc                                   <  1:1342269/132‑1 (MQ=21)                                                            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAggg                           >  2:1259251/1‑140 (MQ=255)                                                            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAggg                           >  2:1050333/1‑140 (MQ=21)                                                             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                          <  2:242886/140‑1 (MQ=21)                                                                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCg                    >  1:1225021/1‑140 (MQ=255)                                                                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGa              >  1:1520707/1‑140 (MQ=255)                                                                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAccc          >  1:44822/1‑140 (MQ=255)                                                                              tGCACATATTGTGAGGTCGATGCGTAGGGTGGGACACCTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcccc         <  1:424237/140‑1 (MQ=11)                                                                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGc      <  1:1050333/140‑1 (MQ=21)                                                                                  caTATTGTGAGCTCGGTGAGTATGGGGGGGCACGTGGGATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATGCTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCg     <  1:756363/140‑1 (MQ=11)                                                                                    tatTGTGAGCTCGATGAGTAGGGCGGGACTCGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGGCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg   <  1:315279/140‑1 (MQ=255)                                                                                     atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAgg  >  1:32105/1‑140 (MQ=14)                                                                                     atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑TG‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAgg  >  1:958842/1‑140 (MQ=14)                                                                                                                                        |                                                                                              AGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATA‑T‑‑G‑GGG‑G‑GACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGG  >  NC_000913/3943975‑3944197

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