Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	43.3%	91.5 / 22.5	60	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/26);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.45e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACG  >  NC_000913/2337576‑2337869                                                                                                                                                    |                                                                                                                                                   gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                 <  2:252234/151‑1 (MQ=255)  cgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCtt                                                                                                                                                <  2:459027/151‑1 (MQ=255)            gAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacg                                                                                                                                      >  2:779953/1‑151 (MQ=255)                    tCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTgcg                                                                                                                              >  2:247128/1‑151 (MQ=255)                           aCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAAAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGt                                                                                                                        >  2:806569/1‑150 (MQ=255)                               aTAAACTTTCATCGAATAGACGCGACCCCGGCTGGAGAAGCAAAGAATAGGGTCGGGAGTGTTCCCCACCAGCAGTCGGTCAATAAAGTCTTCTTTTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                      <  2:498250/148‑1 (MQ=255)                                tAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                    >  1:687556/1‑149 (MQ=255)                                 aaaCTTTCATCGAATAGACGCGCCCACGGCTGGAGAAGCACAAAATATGGTCGTGGTTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                  <  2:633221/150‑1 (MQ=255)                                 aaaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                    >  1:679259/1‑148 (MQ=255)                                  aaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGAGCGGCAGATTTACGTTTCCCGCCACGACGCTGCGCTTTGTATTCAGa                                                                                                                 >  1:523294/1‑150 (MQ=255)                                      ttCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCg                                                                                                            >  2:726286/1‑151 (MQ=255)                                      ttCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCg                                                                                                            >  2:506479/1‑151 (MQ=255)                                             aaTAGCCGCGACCACGGCTGGAGAAGCACAGAATATGGTCGGGAGTTTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                        <  2:605812/148‑1 (MQ=255)                                             aaTAGACCCGCCCCCGGCGGGAGAAGAACAGAATATGGTCGTGGGTGTCCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTAAATACGTGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATa                                                                                                     <  1:322302/151‑1 (MQ=255)                                                 gACGCGCCCACGGCGGGAGAAGCACAAAATATGGTCGTGTGTGTCCCCCACAAGCAGTCGTTCAAAAAAGTCTTCTTTTTTAATACTTGCGGAAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                  <  1:903176/150‑1 (MQ=255)                                                   cgcgCCCCCGCCTGGAAAAGCACAGAATATGGTCGTGATTGTCCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTaa                                                                                                <  2:68276/150‑1 (MQ=255)                                                        ccACGGCTGGAGAACCACAGAATATGGTCGTGAGTGTCCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg                                                                                              <  2:372198/147‑1 (MQ=255)                                                           cGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt                                                                                       >  2:572527/1‑151 (MQ=255)                                                               tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                   >  1:292623/1‑151 (MQ=255)                                                                ggAGAAGCACAGAATTTGGTCGTGAGTGTTCCCCCCCAGCAGTCGGTCAAAAAATTCTTCTTTTTTAATCCGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTgg                                                                                     <  1:764776/148‑1 (MQ=255)                                                                ggAAAACCACAAAATATGGTCGTGATTTTTCCCCCCCACCTGTCGTTCAAAAAATTTTTCTTTTTTAATCCGTGCGGAAATTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTga                                                                                  <  1:592566/151‑1 (MQ=255)                                                                 gagaAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                   >  2:768812/1‑149 (MQ=255)                                                                  agaAGCACAGAATATGGTCGTGAGTTTTCCCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgaga                                                                                <  2:225658/151‑1 (MQ=255)                                                                    aaGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTTCGGCAGATTTACCTTTCCCGCCACGAAGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACCTAGCCCTGGTgcg                                                                                 >  1:563807/1‑146 (MQ=255)                                                                    aaGCACAGAATAGGGTCGTGAGTGTCCCCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGTGCGGAAGATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagaga                                                                              <  1:173565/151‑1 (MQ=255)                                                                    aaGCAAAAAATATGGTGGTGAGTGTCCGCCACCAGCAGTGGGTCGAAAATGTCTTTTTCTTTAATGCGGGGGGAGGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgaga                                                                                <  1:779330/149‑1 (MQ=255)                                                                    aaGCAAAAAATAGGGTGGTGAGTGTTCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATGCGTGGGACATATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag                                                                               <  1:37665/150‑1 (MQ=255)                                                                          aaattatGGTCTTGGGTGTCCCCCCCCAGGAGTCGCTCAATATCTTTTTCTTCTTTATTCCGTCCGTCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                        <  1:506500/147‑1 (MQ=255)                                                                                     gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCt                                                             >  2:426922/1‑151 (MQ=255)                                                                                            ttCGCCCCCAGCAGTCGGTCAATAAAGTTTTCTTTTTAAATCCGTGCGGCAGTTTCCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgg                                                       <  1:886286/150‑1 (MQ=255)                                                                                            ttCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                      >  1:376024/1‑151 (MQ=255)                                                                                            tccGCCCCCAGCAGTCGGTCGATAAATTCTTCTTCTTTAATACGGGCGCAAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                         <  1:534124/146‑1 (MQ=255)                                                                                             tCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                     >  2:748710/1‑151 (MQ=255)                                                                                               gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                    >  2:406781/1‑150 (MQ=255)                                                                                               ccccccaGCAGCCGGTCAATAAGTTCTTCTTTTTTAATCCGTGGGGCAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                    <  1:782289/146‑1 (MQ=255)                                                                                                ccaccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGAt                                                  >  2:150583/1‑151 (MQ=255)                                                                                                 caccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                 >  2:169513/1‑151 (MQ=255)                                                                                                   ccaGCAGTCGGTCGATAAATTCTTCTTTTTTAATACGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCa                                                <  1:491992/150‑1 (MQ=255)                                                                                                    cggcagTCGTTCGAAAAGGTCTTCTTCTTTAATACGTGCGGCATATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCa                                                <  1:680539/147‑1 (MQ=255)                                                                                                    cagcagTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                 >  2:708634/1‑148 (MQ=255)                                                                                                        agTCGTTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGAt                                             <  2:857144/148‑1 (MQ=255)                                                                                                        agTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                          >  1:887618/1‑151 (MQ=255)                                                                                                        agTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                          >  2:363400/1‑151 (MQ=255)                                                                                                          tcggtcgAAAAAGTCTTCTTCTTAAATCCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTc                                         <  1:748671/150‑1 (MQ=255)                                                                                                           cgatggagtaggttttcttcttTAAAGCGTCGGTGAGTTTTACGTTGCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCTTCACGACCACATCTTCCTGGGTGATCAGATc                                            <  1:186092/132‑1 (MQ=255)                                                                                                            ggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                       <  1:308339/150‑1 (MQ=255)                                                                                                            ggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                        >  2:220347/1‑149 (MQ=255)                                                                                                            ggtcgATAAAGTCTTCTTCTTTAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                       <  2:155074/150‑1 (MQ=255)                                                                                                            ggtcgATAAAGTCTTCTTCTTTAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                       <  1:708595/150‑1 (MQ=255)                                                                                                             gtcaATAAAGTCTTTTTTTTTAATACGGGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                        <  1:806529/148‑1 (MQ=255)                                                                                                                          tcttcttTAATCCGTGGGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                         <  2:11400/150‑1 (MQ=255)                                                                                                                             tgtttATTTCGTGGGCCATATTTCCCTTTTCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAGGTTTg                           <  1:350289/143‑1 (MQ=255)                                                                                                                             tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                     >  1:904354/1‑151 (MQ=255)                                                                                                                                   aTACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg               >  2:50795/1‑151 (MQ=255)                                                                                                                                     aCGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTg             >  2:153168/1‑151 (MQ=255)                                                                                                                                        tGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAtt          >  2:777112/1‑151 (MQ=255)                                                                                                                                         ggggCAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAttt         <  2:246386/149‑1 (MQ=255)                                                                                                                                              aGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt     >  2:328476/1‑150 (MQ=255)                                                                                                                                                aTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg  >  2:802891/1‑151 (MQ=255)                                                                                                                                                 tttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt     <  2:215224/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACG  >  NC_000913/2337576‑2337869

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