Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,157 G→A L212L (CTG→TTG)  sthA ← soluble pyridine nucleotide transhydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,1570GA98.0% 151.7 / ‑6.8 50L212L (CTG→TTG) sthAsoluble pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/30);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGG  >  NC_000913/4159917‑4160403
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
gTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGAt                                                                                                                                                                                                                                                 <  1:194599/249‑1 (MQ=255)
     cgGCTGTCAGTTTCTAGCACCATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGCGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTCCGCCCCTGTTCCAGAAGTGGTAGGAGAGAGATTCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATAGCgggg                                                                                                                                                                                                                                         >  2:216638/1‑250 (MQ=255)
             aGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGACGTGATAGGAGAGAGAATCTGCCATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTgctc                                                                                                                                                                                                                                  >  1:495548/1‑247 (MQ=255)
             aGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCa                                                                                                                                                                                                                                 <  1:176238/251‑1 (MQ=255)
                   gagtccAATGTGCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGGCAGCTTTCAGTGTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTTACATCTCTTGATCGAGAAATGCCAACAGGCTGTCGCGGGGGTTGATCAGATCCa                                                                                                                                                                                                                           <  2:70681/247‑1 (MQ=255)
                    atcccAATGTTCTGGAACGCCAGCGAACCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                                                                                                                                                           <  2:54175/248‑1 (MQ=255)
                        cAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                                                                                                                                                      >  1:292667/1‑251 (MQ=255)
                         aaTGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                                                                                                                                                     <  1:308527/251‑1 (MQ=255)
                               cTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGAATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGACGTGATAGGAGAGAGAATCTGACATCTCTGGATCGCGAACTGCCAACAGGAGATCGCGGGTGATGATAAGATCCATTTTTACATcac                                                                                                                                                                                                               >  2:140804/1‑249 (MQ=255)
                                 gTAACGGCATCGAATCGTTATGTCCGGTGCGCCCGTTGGCATATAGCAGGCAGTAAGCTTTCAGTGTGTTACCCGACTTCAGAAGCATGAGCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCCCTTTTACATCCATa                                                                                                                                                                                                             <  2:207085/251‑1 (MQ=255)
                                   aaCGGCAGCGAATCGGTTTTTCCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTGCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTGCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACAGCCATAc                                                                                                                                                                                                            <  2:430184/250‑1 (MQ=255)
                                          gCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAg                                                                                                                                                                                                    <  1:264019/251‑1 (MQ=255)
                                                 ggTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                                                                                                                                                      <  1:524572/242‑1 (MQ=255)
                                                 ggTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACg                                                                                                                                                                                             <  1:136974/251‑1 (MQ=255)
                                                 ggTATTACCGGAGCGACCGTTGGCATAGAGCAGGCTGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCACGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGCTCgga                                                                                                                                                                                              >  2:524572/1‑248 (MQ=255)
                                                  gTATTACCGGTGCGACCGTTGGCATCGAGAAGGCAGTCAGCTTTCAGTTTTTTACCCGTCTTCAGATGCATGATCACACCATCGTCACAACCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGCAATGCCAAGAGGCGGTCGCGGGCGTTGATCGGATCCACTTTTACATCCATACCGCGGCAGATCGACGc                                                                                                                                                                                            >  1:361219/1‑251 (MQ=255)
                                                              cGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                                                                                                                                                                <  1:444952/251‑1 (MQ=255)
                                                                  cGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCCTACCGCGGAAGATCGACGCATATTCACAGCCGATc                                                                                                                                                                            >  2:443054/1‑251 (MQ=255)
                                                                        cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTCCATCCATACCGCGGACGAACGAAGCATATTCACAGCCggaccctcac                                                                                                                                                                      >  2:560352/1‑242 (MQ=255)
                                                                         aTAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAg                                                                                                                                                                     <  1:340613/251‑1 (MQ=255)
                                                                          aagagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACTCATATTCACAGCCGATCACTCCAGc                                                                                                                                                                    <  2:292667/250‑1 (MQ=255)
                                                                           agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGc                                                                                                                                                                    >  1:427024/1‑250 (MQ=255)
                                                                             agCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACTCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAcc                                                                                                                                                                 <  1:299952/251‑1 (MQ=255)
                                                                                   cagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCAt                                                                                                                                                               >  1:563858/1‑247 (MQ=255)
                                                                                   cagTCAGCTTTCAGTTTTTTACCCGACGTCAGATGCATGATCACACCATCGTCACAGCCTGCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAt                                                                                                                                                           <  2:427024/251‑1 (MQ=255)
                                                                                                 ttttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGCGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGACGATCGACGCATATTAACAGCCGATCACTCCAGCACCATAGATAAGAACATGGcgcg                                                                                                                                             >  2:115106/1‑251 (MQ=255)
                                                                                                      tACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCCCTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCTCTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCGCCATAGATAAGTACATGGCGCGGTTc                                                                                                                                         <  1:212713/250‑1 (MQ=255)
                                                                                                       aaccGACTTCAGATGCATGATCCCACCATCGGCACAGCCTTCGATCTTCTCGTACTCTGCGGTTTGACGAATCACTACGCCACTGTTCCAGAAGTGATATGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCTCGGAATAGCGACGCATATTCACAGCCGATCACTCCAGCTCCATAGATAAGGTCATGGCGCGGTTCgt                                                                                                                                       <  2:564154/249‑1 (MQ=255)
                                                                                                                 aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCt                                                                                                                             >  1:6969/1‑251 (MQ=255)
                                                                                                                                ccATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCt                                                                                                              <  1:247802/251‑1 (MQ=255)
                                                                                                                                  aTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGt                                                                                                            <  1:124447/251‑1 (MQ=255)
                                                                                                                                          aGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATg                                                                                                    <  1:351272/251‑1 (MQ=255)
                                                                                                                                               tCGATCTTATCGTACTCTTCGTTATGACGACTCACTAGGCCACTGTTACAGCAGTGATAGGAGAGCGAATCTGCCATCTCTTGATCGAGAAATGCCAACAGGCGAACACGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCAGCTCACTCCTGCACCATAGATACGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTaaagccgtgg                                                                                               >  2:525338/1‑244 (MQ=255)
                                                                                                                                                      tccagTACTCTTCGTGGTTACAAATCCCTACGCCCATGTTCCAGAAGTGAAGGGGGAGAGAATCTGTAATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGTGGTTTGATCAGTTCCACTTTTACAACCAGACCTCTGAAGATCGTCGCATAGTCACAGCCGAGCTCTCCAGAACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGt                                                                                         <  2:204521/246‑1 (MQ=255)
                                                                                                                                                       ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTg                                                                                        >  1:230099/1‑250 (MQ=255)
                                                                                                                                                          gTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAt                                                                                    <  1:436701/251‑1 (MQ=255)
                                                                                                                                                                 gctgtgtgACGAATAACGACGCCACTGTTCCAGAAGTGATAGTAGAGGGAATCTGACATCTCTTGATCGTGAATTGCCAACAGGCGACTGTGGGTGGTGATCAGATCCATTTTTACATCCATACCTCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGGGCATTCTGAGAATTGAGTCGCTTTCGTAAATGCGTGGATGGGTGTAATCAACAt                                                                              <  2:566849/246‑1 (MQ=255)
                                                                                                                                                                  cgtgttgACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGAGCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAGCGACGCATATGCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCt                                                                            <  2:280629/246‑1 (MQ=255)
                                                                                                                                                                  cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGTTAAGTACATGGCGCGGTTCGTGGAGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGATTGGGTGAAATCAAAATCt                                                                            >  2:429626/1‑251 (MQ=255)
                                                                                                                                                                          cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgatg                                                                    <  1:167368/251‑1 (MQ=255)
                                                                                                                                                                                      ccACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACgaga                                                        <  1:355418/251‑1 (MQ=255)
                                                                                                                                                                                              ccAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAg                                                 <  1:472764/250‑1 (MQ=255)
                                                                                                                                                                                                         aGGAGAGAGATTCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCATATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCGTGGATAAGTACCTGGCGCGGTTCGTGGTGCATGCTGAGTATTGCGTCGCTGTCTTAAATACGTGGATGGGGGAAATCAACATCTGTTGGATGCTATGTACGAAGGCCGCAGGCATTAACAAAt                                     >  2:560760/1‑251 (MQ=255)
                                                                                                                                                                                                                 gaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGc                             >  1:158491/1‑251 (MQ=255)
                                                                                                                                                                                                                      tgccATCTCTTGATCGAGAAATGCCCAGAGGCGATCGCGGGAGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGACGATCACGCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTAGTAACTGCGTGGAGGGGGGAAATCACCATCTGTTGGATGATATGAACGAGAGCCGCAGGCGATAACAAATTTTTCAGGGGtt                         >  2:34756/4‑250 (MQ=255)
                                                                                                                                                                                                                            cgcgtgtgcgataaaTGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTGACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACGCCAGCACCATAGATTAGTACATTGCGCGGTTCGTGGTGCATGCTGAGAATGGAGTCGCTGTCGTAAATGCGTGGAGGGGTGAAAGCAACATCTGTTGGATGATATGGACGATAGTCGCAGGCAATAACGAATTTTTCAGCTGTTAgtggt                   <  2:70805/238‑3 (MQ=255)
                                                                                                                                                                                                                                    cGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTg          <  1:7813/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gaAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCCGCACCCTTGATAAGTACATGGCGCGGTTCGTGGTGAATGCTGAGAATTGCGACGCTGTCGTAAATGCGTGGGTGGGGGAAATCAACATCTGTGTGATGATATGGACAATAGCCGTAGGACATAACAAATTTTTCTGCGAATAGTGTTTAAACGTTGCCg       >  2:308681/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCCAACAGGCGATCGCGGGTGTTGATCCGATCCACTTTTCCATCCATACCGAGGAAGATCGACGCATATTCACAGCCGATCACTCCAACACCATAGCTAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGCTATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCgg  >  1:330646/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCCAACAGGCGATCGCGGGTGTTGATCAGATCCCCTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCAGAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCg   <  2:258826/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGG  >  NC_000913/4159917‑4160403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: