Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,692,976	T→C	5.1%	T430A (ACG→GCG)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,692,976	0	T	C	5.1%	267.1 / 3.1	79	T430A (ACG→GCG)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand):  ref base T (45/30);  new base C (2/2);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

CATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAATT  >  NC_000913/1692831‑1693124                                                                                                                                                  |                                                                                                                                                     cATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGa                                                                                                                                                    >  1:168117/1‑148 (MQ=255)   tGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCg                                                                                                                                                 <  4:724962/149‑1 (MQ=255)         cTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATaa                                                                                                                                            >  2:91749/1‑148 (MQ=255)           aTTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaa                                                                                                                                         >  3:742693/1‑148 (MQ=255)            tttttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaa                                                                                                                                        <  3:874713/149‑3 (MQ=255)               ttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaattc                                                                                                                                     <  4:742693/149‑6 (MQ=255)                  tcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGataa                                                                                                                                  >  4:311802/1‑149 (MQ=255)                   cATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGataa                                                                                                                                  >  1:450115/1‑148 (MQ=255)                     caGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGAtaataa                                                                                                                               <  3:78665/148‑1 (MQ=255)                      aGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGAtaataa                                                                                                                               >  1:558512/1‑148 (MQ=255)                                ttGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAg                                                                                                                     >  1:903582/1‑148 (MQ=255)                                       tCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTg                                                                                                                                                     <  2:375093/109‑1 (MQ=255)                                       tCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTg                                                                                                                                                     >  1:375093/1‑109 (MQ=255)                                       tCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGAtaataa                                                                                                                               >  1:391320/1‑131 (MQ=255)                                       tCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGAtaataa                                                                                                                               <  2:391320/131‑1 (MQ=255)                                           attattTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATcc                                                                                                         >  3:81020/1‑149 (MQ=255)                                                 tttATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAg                                                                                                    >  1:11217/1‑148 (MQ=255)                                                 tttATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaa                                                                                                                                         >  1:768825/1‑110 (MQ=255)                                                     ttGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACt                                                                                               >  4:28373/1‑149 (MQ=255)                                                       gAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTa                                                                                              <  2:787250/148‑1 (MQ=255)                                                        aaTATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                             <  1:405571/148‑1 (MQ=255)                                                        aaTATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                             <  2:558512/148‑1 (MQ=255)                                                         atatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAg                                                                                            >  1:190654/1‑148 (MQ=255)                                                          tatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAgg                                                                                           >  2:792450/1‑148 (MQ=255)                                                          tatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAgg                                                                                           >  2:656685/1‑148 (MQ=255)                                                                 taGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGa                                                                                    <  1:254128/148‑1 (MQ=255)                                                                   gTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtgt                                                                                 >  3:858372/1‑149 (MQ=255)                                                                   gTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtgt                                                                                 >  4:131847/1‑149 (MQ=255)                                                                       tatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCg                                                                              >  2:202447/1‑148 (MQ=255)                                                                        atCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCgg                                                                             >  2:153608/1‑148 (MQ=255)                                                                        atCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGa                                                                            <  3:311802/149‑1 (MQ=255)                                                                        atCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGa                                                                            >  3:341661/1‑149 (MQ=255)                                                                           gCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                         >  3:749206/1‑149 (MQ=255)                                                                              gATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGccc                                                                       >  1:666967/1‑148 (MQ=255)                                                                               aTTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCat                                                                     >  3:260980/1‑149 (MQ=255)                                                                                ttAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAACCCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                         <  4:522163/144‑1 (MQ=255)                                                                                ttAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                         >  3:522163/1‑144 (MQ=255)                                                                                ttAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCata                                                                    >  3:121847/1‑149 (MQ=255)                                                                                ttAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCat                                                                     >  1:122807/1‑148 (MQ=255)                                                                                ttAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCat                                                                     >  2:407457/1‑148 (MQ=255)                                                                                 tAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCata                                                                    <  2:740349/148‑1 (MQ=255)                                                                                  aaTTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataat                                                                  <  3:131847/149‑1 (MQ=255)                                                                                     tgctgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCAtaataa                                                                <  1:656685/148‑1 (MQ=255)                                                                                           tgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCgg                                                         >  4:823180/1‑149 (MQ=255)                                                                                            gcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCgg                                                         >  2:735921/1‑148 (MQ=255)                                                                                                ttttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataa                                                                   >  1:431730/1‑134 (MQ=255)                                                                                                ttttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataa                                                                   <  2:431730/134‑1 (MQ=255)                                                                                                ttttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGttt                                                    >  3:751187/1‑149 (MQ=255)                                                                                                   tttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCa                                                 <  4:858372/149‑1 (MQ=255)                                                                                                    ttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCa                                                 >  1:296267/1‑148 (MQ=255)                                                                                                        gATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGa                                            >  3:37776/1‑149 (MQ=255)                                                                                                        gATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGa                                            >  4:622360/1‑149 (MQ=255)                                                                                                        gATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGa                                            >  4:107894/1‑149 (MQ=255)                                                                                                          ttATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGAt                                           >  1:915089/1‑148 (MQ=255)                                                                                                              cAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTg                                             <  1:118990/142‑1 (MQ=255)                                                                                                              cAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTg                                             >  2:118990/1‑142 (MQ=255)                                                                                                               aaTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCg                                      <  2:168117/148‑1 (MQ=255)                                                                                                               aaTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCg                                      <  1:153608/148‑1 (MQ=255)                                                                                                               aaTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGa                                     >  3:207944/1‑149 (MQ=255)                                                                                                                aTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGat                                    >  3:368525/1‑149 (MQ=255)                                                                                                                  ttCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACC‑GATTGGCGatat                                  >  1:666964/1‑148 (MQ=255)                                                                                                                      aCCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATcc                              <  4:751187/149‑1 (MQ=255)                                                                                                                      aCCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATc                               <  2:11217/148‑1 (MQ=255)                                                                                                                        cACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGc                            <  3:622360/149‑1 (MQ=255)                                                                                                                        cACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGc                            >  3:146341/1‑149 (MQ=255)                                                                                                                            aTTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTa                         <  2:752590/148‑1 (MQ=255)                                                                                                                                  ttGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCa                   >  2:836888/1‑148 (MQ=255)                                                                                                                                     aaTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAAt               >  3:387989/1‑149 (MQ=255)                                                                                                                                        tttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaa             >  2:808178/1‑148 (MQ=255)                                                                                                                                        tttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaa             <  2:122807/148‑1 (MQ=255)                                                                                                                                         ttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaaa            <  1:202447/148‑1 (MQ=255)                                                                                                                                          tttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAg           <  2:915089/148‑1 (MQ=255)                                                                                                                                            tATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGc         <  1:407457/148‑1 (MQ=255)                                                                                                                                             aTCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAg       <  4:368525/149‑1 (MQ=255)                                                                                                                                             aTCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAg       >  4:745189/1‑149 (MQ=255)                                                                                                                                               ccGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAgg      <  2:903582/148‑1 (MQ=255)                                                                                                                                                cGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGa     <  1:792450/148‑1 (MQ=255)                                                                                                                                                 gTGAGCGGATAAAAGCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAt   >  3:67564/1‑149 (MQ=255)                                                                                                                                                 gTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAt   <  4:341661/149‑1 (MQ=255)                                                                                                                                                  tGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAtt  <  3:28373/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAATT  >  NC_000913/1692831‑1693124

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