Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,231,659 T→G intergenic (‑28/+225) yjbD ← / ← lysC DUF3811 family protein/lysine‑sensitive aspartokinase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,231,6590TG90.7% 154.7 / 6.6 54intergenic (‑28/+225)yjbD/lysCDUF3811 family protein/lysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (25/24);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

TGCTTTTTCTTCAACTGGCGGGCTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATACAGTATTT  >  NC_000913/4231429‑4231885
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
tGCTTTTTCTTCAACTGGCGGGCTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATCCGGGGGAGTGCCATGATGGTGCCCTCAGTAAGCCAGAGAGGTCAggg                                                                                                                                                                                                                               >  2:591676/1‑234 (MQ=255)
tGCTTTTTCTTCAACTGGCGGGCTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAgtg                                                                                                                                                                                                                               >  1:510405/1‑236 (MQ=255)
                 gCGGGCTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCGTGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTg                                                                                                                                                                                                               >  2:74922/1‑235 (MQ=255)
                      cTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCa                                                                                                                                                                                                         <  2:719243/236‑1 (MQ=255)
                              gCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGATGGTTTCGGAGGTGGTCAATGCGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGat                                                                                                                                                                                                 <  2:444767/236‑1 (MQ=255)
                               cTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGa                                                                                                                                                                                                  <  2:428614/234‑1 (MQ=255)
                               cTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGa                                                                                                                                                                                                  >  1:428614/1‑234 (MQ=255)
                                              catccATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCt                                                                                                                                                                                 <  2:142993/233‑1 (MQ=255)
                                               aaCCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTg                                                                                                                                                                                    >  2:401684/1‑232 (MQ=255)
                                                                     tctttcttAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAAc                                                                                                                                                                                                     >  1:318358/1‑193 (MQ=255)
                                                                     tctttcttAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAAc                                                                                                                                                                                                     <  2:318358/193‑1 (MQ=255)
                                                                                   gTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATGCATTTTTCCTTTCCTTAACTATTTATCTCTTTCGATAATTACTTaga                                                                                                                                            >  2:677514/1‑236 (MQ=255)
                                                                                   gTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTCCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGGACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGCTActtaattaca                                                                                                                                            >  2:663735/1‑234 (MQ=255)
                                                                                          tCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAg                                                                                                                                                                                                           <  2:192572/166‑1 (MQ=255)
                                                                                          tCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAg                                                                                                                                                                                                           >  1:192572/1‑166 (MQ=255)
                                                                                                  ggTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATaaaaaa                                                                                                                             >  1:735458/1‑236 (MQ=255)
                                                                                                             gACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGa                                                                                                                                                                                                       <  1:246758/151‑1 (MQ=255)
                                                                                                             gACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGa                                                                                                                                                                                                       >  2:246758/1‑151 (MQ=255)
                                                                                                                         tCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGAt                                                                                                      >  1:474921/1‑236 (MQ=255)
                                                                                                                             cgcCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATaaa                                                                                                   >  1:584386/1‑235 (MQ=255)
                                                                                                                              gcCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGa                                                                                                 >  1:51688/1‑236 (MQ=255)
                                                                                                                                     ccAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGt                                                                                          <  2:735458/236‑1 (MQ=255)
                                                                                                                                          agTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCa                                                                                                                                                                                                                                  >  1:120752/1‑95 (MQ=255)
                                                                                                                                          agTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCa                                                                                                                                                                                                                                  <  2:120752/95‑1 (MQ=255)
                                                                                                                                             gtTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACa                                                                                     <  2:474921/233‑1 (MQ=255)
                                                                                                                                             gtTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACa                                                                                     <  1:401684/233‑1 (MQ=255)
                                                                                                                                              ttttAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATCAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACACAGt                                                                                 >  2:104897/1‑236 (MQ=255)
                                                                                                                                                      tCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTc                                                                         >  2:17162/1‑236 (MQ=255)
                                                                                                                                                        aCGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCtt                                                                                                                                                                           >  1:660834/1‑136 (MQ=255)
                                                                                                                                                        aCGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCtt                                                                                                                                                                           <  2:660834/136‑1 (MQ=255)
                                                                                                                                                        aCGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGt                                                                       >  1:716674/1‑236 (MQ=255)
                                                                                                                                                                cgcgTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATtgat                                                               <  2:584386/236‑1 (MQ=255)
                                                                                                                                                                             ctctTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACa                                                  >  2:187018/1‑236 (MQ=255)
                                                                                                                                                                                         aaTGCGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATa                                      <  2:313272/236‑1 (MQ=255)
                                                                                                                                                                                            aCGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACg                                    <  1:591676/235‑1 (MQ=255)
                                                                                                                                                                                             cGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGc                                  <  1:17162/236‑1 (MQ=255)
                                                                                                                                                                                                  gAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGc                                                                                                                                                                                   >  1:598372/1‑86 (MQ=255)
                                                                                                                                                                                                  gAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGc                                                                                                                                                                                   <  2:598372/86‑1 (MQ=255)
                                                                                                                                                                                                  gAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGaa                             >  1:461524/1‑236 (MQ=255)
                                                                                                                                                                                                     tGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATtgat                                                               <  2:232673/199‑1 (MQ=255)
                                                                                                                                                                                                     tGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATtgat                                                               >  1:232673/1‑199 (MQ=255)
                                                                                                                                                                                                     tGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGc                                  <  2:371810/228‑1 (MQ=255)
                                                                                                                                                                                                     tGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGc                                  >  1:371810/1‑228 (MQ=255)
                                                                                                                                                                                                         atgatgGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAAtata                        <  2:461524/234‑1 (MQ=255)
                                                                                                                                                                                                           gatgGTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCtt                    <  1:74922/236‑1 (MQ=255)
                                                                                                                                                                                                               gTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATa                                      >  2:715445/1‑214 (MQ=255)
                                                                                                                                                                                                               gTGCCCTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATa                                      <  1:715445/214‑1 (MQ=255)
                                                                                                                                                                                                                  cccTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAg                                                >  2:151396/1‑201 (MQ=255)
                                                                                                                                                                                                                  cccTCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAg                                                <  1:151396/201‑1 (MQ=255)
                                                                                                                                                                                                                   catCAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGcc             <  2:716674/233‑1 (MQ=255)
                                                                                                                                                                                                                      cAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATAc         <  1:187018/236‑1 (MQ=255)
                                                                                                                                                                                                                      cAGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATa          <  2:510405/235‑1 (MQ=255)
                                                                                                                                                                                                                       aGTGAGCCAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATACa        >  1:436201/1‑236 (MQ=255)
                                                                                                                                                                                                                             ccAGAGAGTGCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATACAGTAttt  >  1:406651/1‑236 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
TGCTTTTTCTTCAACTGGCGGGCTTTTTTCGCTTCAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCTTCTCTTAACTATTTATCTCTTTCGATAATTAATTAGAGATGCAGATAAAAAAATTAAAGGCAATTATTCTCCGATAAAGACGCTGGTTAACAGAGTACAGGCTCGTTTATTGATGAGCATAGTGACAAGAAAATCAATACGGCCCGAAATATAGCTTCCAGGCCATACAGTATTT  >  NC_000913/4231429‑4231885

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