Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 11.7% intergenic (+39/‑181) metH → / → yjbB cobalamin‑dependent methionine synthase/putative inorganic phosphate export protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,5500GA11.7% 373.9 / 29.7 138intergenic (+39/‑181)metH/yjbBcobalamin‑dependent methionine synthase/putative inorganic phosphate export protein
Reads supporting (aligned to +/- strand):  ref base G (65/55);  new base A (6/10);  total (72/66)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCA  >  NC_000913/4227405‑4227698
                                                                                                                                                 |                                                                                                                                                    
cacaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa                                                                                                                                                  >  1:1180289/1‑150 (MQ=255)
 acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                                                                 >  2:2481328/1‑150 (MQ=255)
 acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                                                                 <  2:2177115/150‑1 (MQ=255)
     aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCCCAAATCTGTCACTTTTCGTTACAACAAACAGGGCGCTCAATGAg                                                                                                                                             >  1:2090929/1‑150 (MQ=255)
       tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                                                           <  1:1156953/150‑1 (MQ=255)
        cAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGc                                                                                                                                          >  1:1005639/1‑150 (MQ=255)
         aGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                                                                         >  1:1042985/1‑150 (MQ=255)
          gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                                                        <  1:1847065/150‑1 (MQ=255)
           cgcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCCGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCt                                                                                                                                       <  1:1891524/150‑1 (MQ=255)
            gcgATCAGGCTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                                                                      <  2:2067460/150‑1 (MQ=255)
             cgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                                                                     <  1:1000930/150‑1 (MQ=255)
              gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                                                                    <  1:1978759/150‑2 (MQ=255)
                tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttg                                                                                                                                  <  1:1213081/150‑4 (MQ=255)
                  aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                                                                <  1:32115/150‑1 (MQ=255)
                  aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCCCAATGAGTGCCCTGTCTCtt                                                                                                                                <  1:1553269/150‑1 (MQ=255)
                   ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                                                               <  2:1867127/150‑7 (MQ=255)
                   ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                                                               <  2:1875939/150‑7 (MQ=255)
                    gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                                                              >  1:2451484/1‑150 (MQ=255)
                      tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                                                            >  2:1489689/1‑150 (MQ=255)
                      tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                                                            <  2:1613333/150‑1 (MQ=255)
                      tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                                                            >  1:583670/1‑150 (MQ=255)
                       gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACAGAATCTGGGGTATGACGCGGTCTGATTTATATTTATGTGACTTTTCCTTTATAAAAAGAGGGCTCTAAATGAGTGCACTTTGTCTttatta                                                                                                                           <  2:814560/150‑1 (MQ=255)
                         acataaTGCCCGCCGTAAAGGTATCCGCGTTAACGAAGTTTAGCGCTGGCTGGCACCGAATATGGGGTATCACGCGGACTGATTCAAAAATCTGTCACTTTTCCTAACAACAAACAGGGCGCTCAATGAGTGCACTGTCTCttaattaaa                                                                                                                         >  2:1950955/6‑149 (MQ=255)
                          gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                                                        <  1:2157088/150‑1 (MQ=255)
                          gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                                                        <  2:944708/150‑1 (MQ=255)
                           tttATGGGCGGCGTAAAGGTATGAGCGTTACAGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGAGTGATTCAGTAATCTTTCACTTTTCCTTACAACAAACAGGTCGCTCAATGAGTGCCCTGTATCTTTATTAattt                                                                                                                       <  1:1950955/149‑3 (MQ=255)
                              aTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                                                                    <  1:2475225/150‑1 (MQ=255)
                               tGCCCGCCGTAAAGGTATGAGAGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCATTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTAATTAATATGaaa                                                                                                                   >  2:564282/1‑150 (MQ=255)
                                gCCCGCTGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                                                                  >  2:2318303/1‑150 (MQ=255)
                                 cccGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                 >  2:1307744/1‑150 (MQ=255)
                                   cgccgTAAATGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                                                               <  2:1542846/150‑1 (MQ=255)
                                   cgccgTAAAGGTATTAGCGGTTCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGAGTCGGACTGATTCACAAATCTTTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGATTCACt                                                                                                               <  1:1567328/150‑1 (MQ=255)
                                   cgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                                                               >  1:1454632/1‑150 (MQ=255)
                                       gTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtata                                                                                                           >  1:676453/1‑150 (MQ=255)
                                       gTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtata                                                                                                           <  1:2417339/150‑1 (MQ=255)
                                        tAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                                                          >  1:351177/1‑150 (MQ=255)
                                          aaGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTTCAACAAACAGGGCGGTCAATGAGTGCGCTGTCTCTTTATTTATATGAAGCGCTTATACTg                                                                                                        >  2:2370976/1‑150 (MQ=255)
                                            ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                                                                      <  1:949616/150‑1 (MQ=255)
                                            ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                                                                      >  1:905344/1‑150 (MQ=255)
                                            ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                                                                      <  2:2316016/150‑1 (MQ=255)
                                               aTGAGCGTTACCGAAGTTGAGTGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                                                                   <  1:1190744/150‑1 (MQ=255)
                                                tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                                                                  >  2:508679/1‑150 (MQ=255)
                                                   gCGTTACCGATGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCATTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                                                                               <  1:2276709/150‑1 (MQ=255)
                                                     gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCGCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                                                                             <  1:1806084/150‑1 (MQ=255)
                                                        aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                                                                          >  1:725533/1‑150 (MQ=255)
                                                        aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTAACAAATATGTCACTATTCCTTAAAAAAAACAGGGCGCTCAATGAGTCGCCTGTCTCTTTATTAATATGAAACAATTATACTCGAAACAGGCTGGaa                                                                                          >  2:2316134/1‑150 (MQ=255)
                                                              gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                                                                    >  1:2490861/1‑150 (MQ=255)
                                                              gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCCATCAATCATTCCCATGTCGATTTATTAATTTGAAAAAATAATAATGGAAACAGGCTGGAATAATCt                                                                                    >  2:1851379/1‑150 (MQ=255)
                                                                tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                                                                  >  2:1811544/1‑150 (MQ=255)
                                                                   gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                                                               <  1:2254999/150‑1 (MQ=255)
                                                                    cgcTGGCTGGCACCGAATCTGGGGTAGGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTAAGCCg                                                                              <  2:1791281/150‑1 (MQ=255)
                                                                       tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                                                                           >  2:1039084/1‑150 (MQ=255)
                                                                            ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                                                                      <  1:465762/150‑1 (MQ=255)
                                                                              cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                                                                    >  2:443279/1‑150 (MQ=255)
                                                                               aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                                                                   <  1:1262364/150‑1 (MQ=255)
                                                                                ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                                                                  <  1:742355/150‑1 (MQ=255)
                                                                                ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                                                                  <  1:459474/150‑1 (MQ=255)
                                                                                 cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                                                                 >  2:1552911/1‑150 (MQ=255)
                                                                                   aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                                                               <  1:1535309/150‑1 (MQ=255)
                                                                                       tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                                                           >  2:2199766/1‑150 (MQ=255)
                                                                                       tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                                                           >  1:1959876/1‑150 (MQ=255)
                                                                                          ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca                                                        >  1:1487803/1‑150 (MQ=255)
                                                                                           gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                                                       <  2:726948/150‑1 (MQ=255)
                                                                                           gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                                                       <  1:2287436/150‑1 (MQ=255)
                                                                                                     gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg                                             <  1:2043516/150‑1 (MQ=255)
                                                                                                       cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc                                           <  2:1693756/150‑1 (MQ=255)
                                                                                                         gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg                                         >  1:191071/1‑150 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                     >  2:1434658/1‑150 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                     >  1:599821/1‑150 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                     >  1:174271/1‑150 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                     >  1:748840/1‑150 (MQ=255)
                                                                                                              acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTCTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGc                                    >  1:1017002/1‑150 (MQ=255)
                                                                                                               caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGATACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCg                                   <  2:1127158/150‑1 (MQ=255)
                                                                                                               caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCg                                   >  1:1389059/1‑150 (MQ=255)
                                                                                                                        cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTcctg                          <  1:772722/150‑3 (MQ=255)
                                                                                                                         aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGt                         <  2:2349032/150‑1 (MQ=255)
                                                                                                                          cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                        >  1:217437/1‑150 (MQ=255)
                                                                                                                                cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                  >  1:1107548/1‑150 (MQ=255)
                                                                                                                                cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                  <  2:351177/150‑1 (MQ=255)
                                                                                                                                  tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATACTCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTc                >  1:554380/1‑150 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt              >  2:2100610/1‑150 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt              >  1:1085374/1‑150 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt              >  1:1702881/1‑150 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt             >  2:1232995/1‑150 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGAGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt              >  1:2140997/1‑150 (MQ=255)
                                                                                                                                     aacaaacaGGTCGCTCAATGAGTGCGCTGTGTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATGTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTTGTCAtt             <  1:1433134/150‑1 (MQ=255)
                                                                                                                                      acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc            >  1:876605/1‑150 (MQ=255)
                                                                                                                                        agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCaa          >  1:165913/3‑150 (MQ=255)
                                                                                                                                         aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAt         <  2:1072955/150‑1 (MQ=255)
                                                                                                                                          acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAtt        >  1:1221171/1‑150 (MQ=255)
                                                                                                                                          acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAtt        <  2:1005639/150‑1 (MQ=255)
                                                                                                                                          acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAtt        <  2:599821/150‑1 (MQ=255)
                                                                                                                                           caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg       >  1:1488044/1‑150 (MQ=255)
                                                                                                                                           caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg       >  1:1741204/1‑150 (MQ=255)
                                                                                                                                             gggCGCTCAATGAGTGCCCTGTCTCTCTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt     <  2:1017002/150‑1 (MQ=255)
                                                                                                                                             gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt     >  2:942095/1‑150 (MQ=255)
                                                                                                                                               gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTc   <  2:671236/150‑1 (MQ=255)
                                                                                                                                               gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTc   <  1:1307744/150‑1 (MQ=255)
                                                                                                                                                cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCgg                          >  1:112635/1‑126 (MQ=255)
                                                                                                                                                cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCgg                          <  2:112635/126‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCA  >  NC_000913/4227405‑4227698

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