Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,132,434 T→A 6.5% intergenic (+167/‑182) metL → / → metF Bifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,132,4340TA 6.5% 385.2 / 11.7 132intergenic (+167/‑182)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (50/73);  new base A (5/4);  total (55/77)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

Only 100 of 132 total aligned reads displayed.

GCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACA  >  NC_000913/4132286‑4132582
                                                                                                                                                    |                                                                                                                                                    
gCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                                                                                                                                     >  2:309879/1‑150 (MQ=255)
 cccGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc                                                                                                                                                    >  1:179063/1‑150 (MQ=255)
   cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGTCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                                                                                                                                  <  1:505686/150‑1 (MQ=255)
   cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                                                                                                                                  >  2:1628525/1‑150 (MQ=255)
       aCTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTATGCAAAGCAAACCGGact                                                                                                                                              <  2:1489841/150‑3 (MQ=255)
        cTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                                                                                                                                             >  2:1028072/1‑150 (MQ=255)
         tCATGTTTTCGTGTTTATGGTTTCTAATGAAATATATTGATTTATCATAGGTTTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                            <  1:195087/150‑1 (MQ=255)
         tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAATAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                                                                                                                                            >  1:1599612/1‑150 (MQ=255)
           aTGTTTTCTGGTTTATGGTTTCTAATGAAATATATTGATTTATCCTAGGATTAGGGCGGATTTAGCGTTTACGACGTATCCGGCAAGAAGCACTAAGTACATGGTTTGCTTATATTTGCTGTCCGGTTTGCTTTGCATACCGGTttttct                                                                                                                                          <  2:1654237/150‑1 (MQ=255)
           aTGTTTTCGGGTTTATGTTTTCTAATGAAATATATTGAATTATCATAGTATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATTTTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                                                                                                                                          <  1:1156317/150‑1 (MQ=255)
           aTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                                                                                                                                          <  2:1555050/150‑1 (MQ=255)
           aTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                                                                                                                                          <  2:1076988/150‑1 (MQ=255)
                ttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                                                                                                                                     <  1:1333196/150‑1 (MQ=255)
                 tCGGGTTTATGGTTTCTTATGAAATATATTGTTTTATGATTGGATTAGGCCGGATTAATCGTTTACGACGAATCCGGCAAGTAGCAATAAGTACATGGTTAGTTTTTATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttct                                                                                                                                    <  2:1175938/150‑1 (MQ=255)
                  cGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTCCGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt                                                                                                                                   <  1:491261/150‑1 (MQ=255)
                   gggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa                                                                                                                                  <  2:570698/150‑1 (MQ=255)
                    ggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGATGCAATAAGTCCATGGTTAGTTGATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTaa                                                                                                                                 >  2:1060465/1‑149 (MQ=255)
                        tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                             >  1:1513151/1‑150 (MQ=255)
                        tATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc                                                                                                                             >  1:584879/1‑150 (MQ=255)
                           ggTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACAACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg                                                                                                                          >  2:757831/1‑150 (MQ=255)
                              ttCTAATGAAATATATTGAATTTTCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATTCCGGATTTTCTTTTTCTTACCATCCTGAAg                                                                                                                       <  1:1296835/150‑1 (MQ=255)
                                 tAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttt                                                                                                                    >  2:36077/1‑150 (MQ=255)
                                 tAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttt                                                                                                                    >  2:750818/1‑150 (MQ=255)
                                  aaTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttt                                                                                                                   >  2:1800859/1‑150 (MQ=255)
                                   aTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttttt                                                                                                                  <  1:658961/150‑1 (MQ=255)
                                     gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttttta                                                                                                                >  1:1374733/1‑149 (MQ=255)
                                     gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                                >  1:737444/1‑150 (MQ=255)
                                     gAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTc                                                                                                                >  2:1701045/1‑150 (MQ=255)
                                        atatatTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACTAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATATGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACAATCCTGAAGTTTTTTCATc                                                                                                             >  1:243822/1‑150 (MQ=255)
                                          atatTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCATGAAGCATTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTATCTTATTCTTACCATCCTGAAGTTTTTTCATCtt                                                                                                           >  1:334902/1‑150 (MQ=255)
                                            atTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                         >  1:1494316/1‑150 (MQ=255)
                                            atTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTcc                                                                                                         >  2:85497/1‑150 (MQ=255)
                                             ttGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTccc                                                                                                        <  1:548050/150‑1 (MQ=255)
                                               gAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTg                                                                                                      >  1:108918/1‑150 (MQ=255)
                                                aaTTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGa                                                                                                     <  2:1172396/150‑1 (MQ=255)
                                                     tCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGAttttt                                                                                                <  2:900410/150‑1 (MQ=255)
                                                       aTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTcc                                                                                              <  1:1257613/150‑1 (MQ=255)
                                                       aTAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTcc                                                                                              >  1:126492/1‑150 (MQ=255)
                                                                  ccGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCAtt                                                                                   >  1:482510/1‑150 (MQ=255)
                                                                  ccGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCAtt                                                                                   >  2:1123425/1‑150 (MQ=255)
                                                                   cGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTg                                                                                  <  2:1277068/150‑1 (MQ=255)
                                                                     gATTTAGCGTTTACGACGAATCCAGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGATTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGt                                                                                <  1:1379233/150‑1 (MQ=255)
                                                                     gATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGt                                                                                <  2:1620592/150‑1 (MQ=255)
                                                                      aTTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTc                                                                               <  1:66284/150‑1 (MQ=255)
                                                                             gTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTc                                                                        <  2:1550332/150‑1 (MQ=255)
                                                                              tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTTGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                       >  2:515700/1‑150 (MQ=255)
                                                                              tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                       >  1:898182/1‑150 (MQ=255)
                                                                              tttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCg                                                                       <  1:353803/150‑1 (MQ=255)
                                                                               ttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCgg                                                                      >  2:1435974/1‑150 (MQ=255)
                                                                               ttACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCgg                                                                      <  2:676723/150‑1 (MQ=255)
                                                                                tACGACGTATACGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGt                                                                     <  2:484845/150‑1 (MQ=255)
                                                                                tACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGt                                                                     >  2:495163/1‑150 (MQ=255)
                                                                                tACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGt                                                                     <  2:850730/150‑1 (MQ=255)
                                                                                     cgaATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACg                                                                <  1:889790/150‑1 (MQ=255)
                                                                                       aaTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCAACTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGcc                                                              >  2:764254/1‑150 (MQ=255)
                                                                                        aTCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGccc                                                             <  1:1535197/150‑1 (MQ=255)
                                                                                          ccGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCtt                                                           >  2:902522/1‑150 (MQ=255)
                                                                                          ccGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCtt                                                           >  2:1178280/1‑150 (MQ=255)
                                                                                           cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                          <  2:835081/150‑1 (MQ=255)
                                                                                           cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                          >  1:1281554/1‑150 (MQ=255)
                                                                                           cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTGTTCGGTTGACGCCCTTc                                                          >  2:1242389/1‑150 (MQ=255)
                                                                                           cGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTATCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTc                                                          >  2:1116634/1‑150 (MQ=255)
                                                                                            ggCAAGAAGCAATAAGTACCTGGTTAGGTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCg                                                         <  2:1279874/150‑1 (MQ=255)
                                                                                                    gCAATAAGTACATGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATACTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTGATTTTTCGGTTGACGCCCTTAGGCTTTTCCt                                                >  1:1499367/1‑150 (MQ=255)
                                                                                                     cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                                <  1:1339746/150‑1 (MQ=255)
                                                                                                     cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                                >  2:1107118/1‑150 (MQ=255)
                                                                                                     cAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCt                                                <  2:584879/150‑1 (MQ=255)
                                                                                                       aTAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTc                                              <  1:1028072/150‑1 (MQ=255)
                                                                                                        tAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCa                                             <  1:1214033/150‑1 (MQ=255)
                                                                                                        tAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCa                                             <  2:323003/150‑1 (MQ=255)
                                                                                                         aaGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCAt                                            <  2:1826697/150‑1 (MQ=255)
                                                                                                            tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                         <  1:610092/150‑1 (MQ=255)
                                                                                                            tACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCtt                                         <  2:179063/150‑1 (MQ=255)
                                                                                                               aTGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTAc                                      <  1:1628525/150‑1 (MQ=255)
                                                                                                               aTGGTTAGGTTATATTTGCAGTCCGGGTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTAc                                      <  2:1678205/150‑1 (MQ=255)
                                                                                                                  gTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATc                                   <  2:1322509/150‑1 (MQ=255)
                                                                                                                    tAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTg                                 <  1:674858/150‑1 (MQ=255)
                                                                                                                    tAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTATTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTg                                 <  1:252017/150‑1 (MQ=255)
                                                                                                                     aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                                <  1:85497/150‑1 (MQ=255)
                                                                                                                     aGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTgg                                <  1:1604829/150‑1 (MQ=255)
                                                                                                                      gTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACAATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGa                               >  1:1046299/1‑150 (MQ=255)
                                                                                                                      gTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAGTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATATTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCGTCATCTTTACATCTGGa                               >  1:55735/1‑150 (MQ=255)
                                                                                                                        ttATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACg                             <  1:1590230/150‑1 (MQ=255)
                                                                                                                          atatTTTCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATACTGAAGTTTTATCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                           <  1:668860/150‑1 (MQ=255)
                                                                                                                          atatTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTc                           <  1:1201456/150‑1 (MQ=255)
                                                                                                                             aagaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTaa                        >  1:1637193/8‑150 (MQ=255)
                                                                                                                              agaaaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTaaa                       >  2:690399/7‑150 (MQ=255)
                                                                                                                                gCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACg                     <  1:452478/150‑1 (MQ=255)
                                                                                                                                 aaatCCGGTATGCAAAGCAAACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACgg                    >  1:1789020/4‑150 (MQ=255)
                                                                                                                                  aGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGa                   <  1:5236/150‑1 (MQ=255)
                                                                                                                                     ccGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGTTTTTTCCTCACCATCATTGGTCATTTTTAGTTTGACGACCTTAGGCTTTTCCTTCATCTTTACATATGTACTTATAAACGGATaa                >  1:361315/1‑149 (MQ=255)
                                                                                                                                     ccGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAg                >  2:569557/1‑150 (MQ=255)
                                                                                                                                      cGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGa               <  1:309879/150‑1 (MQ=255)
                                                                                                                                      cGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCATTCGGATTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGa               >  1:213469/1‑150 (MQ=255)
                                                                                                                                        gTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGAtg             <  2:963251/150‑1 (MQ=255)
                                                                                                                                           tgctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGc          <  2:773014/150‑1 (MQ=255)
                                                                                                                                            gctttgcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCATGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGca         <  1:968320/150‑1 (MQ=255)
                                                                                                                                                 gcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcacaaca    <  1:750818/150‑1 (MQ=255)
                                                                                                                                                 gcATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGcacaaca    <  1:1123425/150‑1 (MQ=255)
                                                                                                                                                   aTACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTAAATCTGGACGTCTAAACGGATAGATGTGcacaacaca  >  1:1385387/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACA  >  NC_000913/4132286‑4132582

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