Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,634,087	C→A	8.8%	T97T (ACC→ACA)	ytjC →	phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,634,087	0	C	A	8.8%	230.0 / 7.8	68	T97T (ACC→ACA)	ytjC	phosphatase
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base A (3/3);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAG  >  NC_000913/4633948‑4634224                                                                                                                                            |                                                                                                                                          tcatcaGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGAcc                                                                                                                                           <  2:917947/140‑1 (MQ=255)  catcaGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCg                                                                                                                                          <  1:1859387/140‑1 (MQ=255)   atcaGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCga                                                                                                                                         >  2:1125323/1‑140 (MQ=255)     caGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCgaag                                                                                                                                       <  2:1722370/140‑1 (MQ=255)        tAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCgaagaag                                                                                                                                    >  2:461341/1‑140 (MQ=255)        tAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCgaagaag                                                                                                                                    >  2:262300/1‑140 (MQ=255)           cGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGaagaagaag                                                                                                                                 >  1:1207840/1‑140 (MQ=255)             aTTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGTCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAagag                                                                                                                               <  1:1455735/140‑1 (MQ=255)             aTTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAagag                                                                                                                               >  2:1600912/1‑140 (MQ=255)             aTTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAagag                                                                                                                               >  1:1723109/1‑140 (MQ=255)                       cACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTc                                                                                                                     <  1:1463993/140‑1 (MQ=255)                        aCCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCg                                                                                                                    >  1:784404/1‑140 (MQ=255)                        aCCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGG‑AAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCgg                                                                                                                   >  1:1096218/1‑140 (MQ=255)                             gCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGc                                                                                                               <  2:606369/140‑1 (MQ=255)                                  cggcggAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAATATACATATCGATTCTCTGACCGAAGAAGTAGAGAACTGGCGTCGGCAGCTGGTc                                                                                                          >  1:535865/1‑140 (MQ=255)                                  cggcggAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTc                                                                                                          >  2:1213401/1‑140 (MQ=255)                                       gAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     >  1:1500544/1‑140 (MQ=255)                                       gAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     >  1:1500412/1‑140 (MQ=255)                                         aaTCATCGCCCTGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCa                                                                                                   <  1:202735/140‑1 (MQ=255)                                         aaTCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCa                                                                                                   <  1:1030617/140‑1 (MQ=255)                                           tcatcGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCAcc                                                                                                 <  2:206683/140‑1 (MQ=255)                                             atcGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGt                                                                                               >  2:1584072/1‑140 (MQ=255)                                                gCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGa                                                                                            <  1:994697/140‑1 (MQ=255)                                                gCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGa                                                                                            <  2:1647295/140‑1 (MQ=255)                                                 cccAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGAc                                                                                           <  2:1326622/140‑1 (MQ=255)                                                 cccAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGAc                                                                                           <  2:1326435/140‑1 (MQ=255)                                                     ggCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGc                                                                                       >  2:1520542/1‑140 (MQ=255)                                                      gCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCg                                                                                      <  1:1607100/140‑1 (MQ=255)                                                        cTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTa                                                                                    <  1:805679/140‑1 (MQ=255)                                                        cTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTa                                                                                    <  1:1590884/140‑1 (MQ=255)                                                           cGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAATAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTc                                                                                 >  1:1728178/1‑140 (MQ=255)                                                           cGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTc                                                                                 >  1:423712/1‑140 (MQ=255)                                                                  gACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAgg                                                                          <  1:1217490/140‑1 (MQ=255)                                                                  gACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAgg                                                                          <  1:76979/140‑1 (MQ=255)                                                                  gACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAgg                                                                          <  2:502744/140‑1 (MQ=255)                                                                  gACATCAACTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGCTTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAgg                                                                          <  2:1580742/140‑1 (MQ=255)                                                                     atcatcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGa                                                                       >  1:1291356/1‑140 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTc                                                                    <  2:235826/140‑1 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     <  2:1787837/107‑1 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     <  2:1787725/107‑1 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     <  2:1374056/107‑1 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     >  1:1787725/1‑107 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     >  1:1374056/1‑107 (MQ=255)                                                                        atcTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACAGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATgg                                                                                                     >  1:1787837/1‑107 (MQ=255)                                                                          cTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGTCCGAAGTAGAAGTGAACTGGCGTCGGCAGCTGGTCAATGGCTCCGTTGACGGGCGTATTCCTGAAGGCGTGTCaa                                                                  <  1:1394059/140‑1 (MQ=255)                                                                             tGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGc                                                               >  2:1788213/1‑140 (MQ=255)                                                                             tGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGc                                                               >  2:1788187/1‑140 (MQ=255)                                                                               aTTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAg                                                             <  1:1737133/140‑1 (MQ=255)                                                                                tcctcGCCCGCGCGAATTAAACATCGGAGTGCTGGAAAAAAGACCTATCGCGTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAgg                                                            <  2:557798/138‑1 (MQ=255)                                                                                    cGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCt                                                        <  1:1644143/140‑1 (MQ=255)                                                                                    cGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCt                                                        >  1:497180/1‑140 (MQ=255)                                                                                      ccTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATCTCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCa                                                      <  1:770284/140‑1 (MQ=255)                                                                                          cGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGa                                                  <  1:1075119/140‑1 (MQ=255)                                                                                          cGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGa                                                  <  1:1543543/140‑1 (MQ=255)                                                                                                ttAAACATGGGTGTGCTGGAAAAATGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGt                                            <  2:1258557/140‑1 (MQ=255)                                                                                                     cATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATg                                       <  2:726111/140‑1 (MQ=255)                                                                                                        gggTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCgg                                    >  1:511938/1‑140 (MQ=255)                                                                                                         ggTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGc                                   >  2:1213055/1‑140 (MQ=255)                                                                                                             tgCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTg                               >  1:1575532/1‑140 (MQ=255)                                                                                                             tgCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTg                               >  1:43179/1‑140 (MQ=255)                                                                                                                tGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGaa                            >  2:535854/1‑140 (MQ=255)                                                                                                                 ggAAAAAAGACCTCTCGCCTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAAt                           <  2:1900968/140‑1 (MQ=255)                                                                                                                 ggAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAAt                           <  1:234981/140‑1 (MQ=255)                                                                                                                 ggAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAAt                           <  1:846539/140‑1 (MQ=255)                                                                                                                        aGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCg                    <  1:741754/140‑1 (MQ=255)                                                                                                                             tatCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACt               <  1:1600912/140‑1 (MQ=255)                                                                                                                              atCGATTCTCTGACCGAAGAAGAAGTGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACtt              >  2:907440/1‑140 (MQ=255)                                                                                                                                cGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTAc            <  2:370946/140‑1 (MQ=255)                                                                                                                                        tGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGa    >  2:1892915/1‑140 (MQ=255)                                                                                                                                          cccGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAg  <  2:1866360/139‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          TCATCAGTAGCGATTTAGGACGCACCCGGCGTACGGCGGAAATCATCGCCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAG  >  NC_000913/4633948‑4634224

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