Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,309,870 C→A 6.3% G371C (GGC→TGC)  alsA ← D‑allose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,309,8700CA 6.3% 261.8 / 8.1 95G371C (GGC→TGC) alsAD‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (44/45);  new base A (2/4);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-02

Only 100 of 103 total aligned reads displayed.

GATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAC  >  NC_000913/4309721‑4310019
                                                                                                                                                     |                                                                                                                                                     
gaTCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                                                                       >  2:705676/1‑150 (MQ=255)
gaTCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                                                                       >  1:1303805/1‑150 (MQ=255)
  tCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGt                                                                                                                                                     >  2:230868/1‑150 (MQ=255)
   cAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGGTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTc                                                                                                                                                    >  1:946871/1‑150 (MQ=255)
     ggACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCtt                                                                                                                                                  >  1:1167428/1‑150 (MQ=255)
       aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                >  1:1223136/1‑150 (MQ=255)
       aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGAGGTACGTTTCACGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCttta                                                                                                                                                >  1:1211427/1‑149 (MQ=255)
        cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                               >  2:1326195/1‑150 (MQ=255)
          ttCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGa                                                                                                                                             >  1:114869/1‑150 (MQ=255)
            ctgctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACt                                                                                                                                           <  2:1267813/150‑1 (MQ=255)
              gctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGc                                                                                                                                         <  2:1167428/150‑1 (MQ=255)
               ctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCg                                                                                                                                        <  2:1241056/150‑1 (MQ=255)
                tgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                                                                       >  2:1520162/1‑150 (MQ=255)
                tgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                                                                       >  1:534937/1‑150 (MQ=255)
                tgATTTCCCCCGGAGAGTTCGGTGATATTCAGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                                                                       <  1:1562234/150‑1 (MQ=255)
                 gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                      <  1:370946/150‑1 (MQ=255)
                 gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                      <  1:1686292/150‑1 (MQ=255)
                  aTTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCt                                                                                                                                     <  1:1191409/150‑1 (MQ=255)
                    ttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGTAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCtat                                                                                                                                   <  2:1420317/150‑3 (MQ=255)
                     tCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGAt                                                                                                                                  <  2:1223919/150‑1 (MQ=255)
                      cccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtataa                                                                                                                                 <  2:1525471/150‑5 (MQ=255)
                       ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                                                                                <  2:667218/150‑1 (MQ=255)
                        cccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGc                                                                                                                               >  2:1641245/1‑150 (MQ=255)
                        cccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtataaag                                                                                                                               <  2:1341571/150‑7 (MQ=255)
                           ggAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCAt                                                                                                                            >  2:1233670/1‑150 (MQ=255)
                            gagaGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATg                                                                                                                           >  1:902986/1‑150 (MQ=255)
                                     gTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                                                                                  <  1:1520162/150‑1 (MQ=255)
                                     gTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                                                                                  <  1:271875/150‑1 (MQ=255)
                                      tGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGGTGGCGTTGTTTTTAAGCGGTAAGGTGCTGGTCAACTTCATGAAACAAGCCCATCGCGGCTTTATAGCCGCCGTGTTTCAGCCTGCGGCTGATCGCCATGTTCTGAGCgc                                                                                                                 >  2:710352/1‑149 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              <  1:1553670/150‑1 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              >  2:964709/1‑150 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              <  2:1583699/150‑1 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              <  2:963536/150‑1 (MQ=255)
                                            tCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaaa                                                                                                           >  1:877028/1‑150 (MQ=255)
                                            tCTGGTTTACCGAATCACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaaa                                                                                                           <  1:566209/150‑1 (MQ=255)
                                              tGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAg                                                                                                         <  1:1257947/150‑1 (MQ=255)
                                               ggTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGt                                                                                                        >  1:137597/1‑150 (MQ=255)
                                                gTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGtt                                                                                                       >  2:552636/1‑150 (MQ=255)
                                                  ttACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgg                                                                                                     <  1:1399729/150‑1 (MQ=255)
                                                    aCCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgggg                                                                                                   <  2:896508/150‑1 (MQ=255)
                                                    aCCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgggg                                                                                                   >  1:1563820/1‑150 (MQ=255)
                                                      cGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaa                                                                                                 >  1:787590/1‑150 (MQ=255)
                                                      cGAATGACATTTCAGAGCCAGCAGTTCGCGTTGAATTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGCTCGCCATGTTCGGAGCGATGGAAAAGATGGGGaa                                                                                                 >  1:307968/1‑150 (MQ=255)
                                                       gAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaa                                                                                                >  2:343486/1‑150 (MQ=255)
                                                       gAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaa                                                                                                <  2:609204/150‑1 (MQ=255)
                                                            aCATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCTCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAcc                                                                                           >  2:1563713/1‑150 (MQ=255)
                                                              aTTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGt                                                                                         <  1:1432459/150‑1 (MQ=255)
                                                              aTTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGt                                                                                         <  1:964709/150‑1 (MQ=255)
                                                               tttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGtt                                                                                        >  2:581456/1‑150 (MQ=255)
                                                                ttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTa                                                                                       <  1:1233670/150‑1 (MQ=255)
                                                                  cAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATc                                                                                     >  1:654846/1‑150 (MQ=255)
                                                                     cgcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcg                                                                                  >  2:347965/1‑150 (MQ=255)
                                                                      gcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcgg                                                                                 <  1:230868/150‑1 (MQ=255)
                                                                        cagcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcg                                                                               <  1:1566298/150‑1 (MQ=255)
                                                                         agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                                                                              <  2:1563820/150‑1 (MQ=255)
                                                                             gTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                                                                          >  2:1524193/1‑150 (MQ=255)
                                                                                cgcgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAg                                                                       <  2:1232730/150‑1 (MQ=255)
                                                                                  cgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATTAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATAGCCATGTTCTGAGCGATGGAAAATTTGGGGAAAAAAACGTTATCACGGCTGCTTTCAGTg                                                                     >  2:1319981/1‑150 (MQ=255)
                                                                                  cgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGGCTTTCAGACTGCGGCTGATCGCCATGTTCGGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTg                                                                     <  1:1558044/150‑1 (MQ=255)
                                                                                   gTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                    <  2:1223136/150‑1 (MQ=255)
                                                                                      gATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGt                                                                 >  1:483852/1‑150 (MQ=255)
                                                                                         tttCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAg                                                              <  1:1521317/150‑1 (MQ=255)
                                                                                                ggTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATccc                                                       >  1:1379793/1‑150 (MQ=255)
                                                                                                ggTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATccc                                                       <  2:877028/150‑1 (MQ=255)
                                                                                                ggTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATccc                                                       >  1:1050502/1‑150 (MQ=255)
                                                                                                    cGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTa                                                                                       >  2:295470/1‑114 (MQ=255)
                                                                                                    cGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTa                                                                                       <  1:295470/114‑1 (MQ=255)
                                                                                                     gTTGCTCGTCAACTTCATGAAACATGCCCATCGCGCCTTTATCGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGATAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCttttt                                                  <  1:1396817/150‑1 (MQ=255)
                                                                                                       tGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTc                                                >  1:683245/1‑150 (MQ=255)
                                                                                                                aCTTCATGATACAAGCCCATCTCGCCTTAATCGCCGCCGTCTTTCCGACTGCGGCTGATCGTCATGTTCTGAGCGAAGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCGTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATcc                                       <  2:299060/150‑1 (MQ=255)
                                                                                                                 cTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGCTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCa                                      >  1:175655/1‑150 (MQ=255)
                                                                                                                    cATGAAACAAGCCCCTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAACACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAggg                                   <  2:1069061/150‑1 (MQ=255)
                                                                                                                       gAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGCCTGCGGCTGATCGCCATGTTCTGCGCGATGGAAAACTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggggg                                >  1:107609/1‑148 (MQ=255)
                                                                                                                       gAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTg                                <  2:534937/150‑1 (MQ=255)
                                                                                                                       gAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTg                                >  2:4913/1‑150 (MQ=255)
                                                                                                                          aCAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAAc                             <  2:1211427/150‑1 (MQ=255)
                                                                                                                           cAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACg                            <  2:654846/150‑1 (MQ=255)
                                                                                                                            aaGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGATAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGACATCCCTTTTTTCACGGCATCCAGGGGTGAACGt                           >  2:1551694/1‑150 (MQ=255)
                                                                                                                            aaGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGACATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGt                           >  2:1509770/1‑150 (MQ=255)
                                                                                                                               cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                        >  2:1347964/1‑150 (MQ=255)
                                                                                                                                ccATCGCGCCTTTCTAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTCAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       <  2:1537782/150‑1 (MQ=255)
                                                                                                                                ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       <  2:1291257/150‑1 (MQ=255)
                                                                                                                                ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       >  2:795718/1‑150 (MQ=255)
                                                                                                                                ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGGGAACGTGgag                       <  1:343486/150‑1 (MQ=255)
                                                                                                                                ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGATGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       >  2:945560/1‑150 (MQ=255)
                                                                                                                                  aTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagag                     <  2:946871/150‑1 (MQ=255)
                                                                                                                                      cgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGatat                 >  1:671759/1‑150 (MQ=255)
                                                                                                                                      cgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGatat                 >  1:715779/1‑150 (MQ=255)
                                                                                                                                        ccTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCt               >  2:1436277/1‑150 (MQ=255)
                                                                                                                                           ttATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTg            <  1:710352/150‑1 (MQ=255)
                                                                                                                                             ataGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGcc          <  1:1709671/150‑1 (MQ=255)
                                                                                                                                              taGCCGCCGTCTTTCAGACTGCTGCTGATCGCCATGTTCTGTGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTATCACGGCATCCAGGGGGGAACGTGGAGAGATATCTTTGCCa         <  1:1319981/150‑1 (MQ=255)
                                                                                                                                              taGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCTATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCa         <  1:850958/150‑1 (MQ=255)
                                                                                                                                               aGCCGCAGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCAATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCAAt        >  2:1472276/1‑150 (MQ=255)
                                                                                                                                                gccgccGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAtt       >  2:1221882/1‑150 (MQ=255)
                                                                                                                                                gccgcaGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAtt       >  1:1243329/1‑150 (MQ=255)
                                                                                                                                                 ccgccgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTa      >  2:426433/1‑150 (MQ=255)
                                                                                                                                                    ccgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGTAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGa   <  1:1466420/150‑1 (MQ=255)
                                                                                                                                                     cgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAc  <  1:552636/150‑1 (MQ=255)
                                                                                                                                                     cgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAc  <  1:1641245/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAC  >  NC_000913/4309721‑4310019

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