Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,430,144 G→C 8.5% A32G (GCC→GGC)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,430,1440GC 8.5% 316.3 / 18.3 94A32G (GCC→GGC) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (48/38);  new base C (3/5);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

AGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAA  >  NC_000913/2429998‑2430289
                                                                                                                                                  |                                                                                                                                                 
aGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCcggcg                                                                                                                                                 <  4:661425/149‑1 (MQ=255)
 gCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGTcggcgg                                                                                                                                                >  1:622576/1‑149 (MQ=255)
    gCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCAt                                                                                                                                             <  2:248818/149‑1 (MQ=255)
     ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                            >  1:588779/1‑149 (MQ=255)
     ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                            <  2:160080/149‑1 (MQ=255)
     ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATc                                                                                                                                            <  1:55128/149‑1 (MQ=255)
      cAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATca                                                                                                                                           <  3:581087/149‑2 (MQ=255)
        gCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                         <  4:434388/149‑4 (MQ=255)
        gCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                         <  1:656361/149‑4 (MQ=255)
            tGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTTCGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACcg                                                                                                                                     >  4:595026/1‑149 (MQ=255)
              gggTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACcgcg                                                                                                                                   >  4:177276/1‑149 (MQ=255)
                 tAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATg                                                                                                                                >  4:533641/1‑149 (MQ=255)
                    cGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTg                                                                                                                             >  3:254304/1‑149 (MQ=255)
                     gCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGa                                                                                                                            <  4:359155/149‑1 (MQ=255)
                      cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                           <  1:472960/149‑1 (MQ=255)
                        caTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGc                                                                                                                         <  4:490897/149‑1 (MQ=255)
                           gACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAcc                                                                                                                      >  4:433889/1‑149 (MQ=255)
                           gACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAcc                                                                                                                      >  1:458716/1‑149 (MQ=255)
                            aCCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCg                                                                                                                     <  1:160299/149‑1 (MQ=255)
                             ccAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGt                                                                                                                    >  4:525649/1‑149 (MQ=255)
                              cAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGtt                                                                                                                   >  1:707276/1‑149 (MQ=255)
                               aaTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCt                                                                                                                                          <  3:548049/125‑1 (MQ=255)
                               aaTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCt                                                                                                                                          >  4:548049/1‑125 (MQ=255)
                               aaTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                                  >  3:18290/1‑149 (MQ=255)
                                aTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCGGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTaa                                                                                                                 <  2:401362/149‑1 (MQ=255)
                                  gCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAgg                                                                                                               <  2:622576/149‑1 (MQ=255)
                                   cccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGc                                                                                                              >  4:225801/1‑149 (MQ=255)
                                    ccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGca                                                                                                             >  1:384360/1‑149 (MQ=255)
                                     cATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcat                                                                                                            >  3:631471/1‑149 (MQ=255)
                                     cATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcat                                                                                                            >  2:356634/1‑149 (MQ=255)
                                      atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                            <  2:66361/116‑1 (MQ=255)
                                      atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                            >  1:66361/1‑116 (MQ=255)
                                      atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCCCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatc                                                                                                           <  3:525649/149‑1 (MQ=255)
                                       tatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatca                                                                                                          >  1:391110/1‑149 (MQ=255)
                                        atTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                           <  4:593399/131‑1 (MQ=255)
                                        atTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                           >  3:593399/1‑131 (MQ=255)
                                         ttGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATa                                                                                                        >  3:595793/1‑149 (MQ=255)
                                           gCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATaaa                                                                                                      <  4:239895/149‑1 (MQ=255)
                                            cccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAAt                                                                                                     >  4:162183/1‑149 (MQ=255)
                                               tGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                  >  3:561188/1‑149 (MQ=255)
                                                gCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGAc                                                                                                 <  4:595793/149‑1 (MQ=255)
                                                  aaaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTg                                                                                               <  1:635594/149‑1 (MQ=255)
                                                      gCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTa                                                                                           <  1:356634/149‑1 (MQ=255)
                                                        tGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAc                                                                                         >  4:149915/1‑149 (MQ=255)
                                                        tGCATATGGCGAGCTTCAAATACATCGCTCAAAAGCCCGTTCGATTTACGCAAACGGAAGCAGATATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAc                                                                                         >  4:113736/1‑149 (MQ=255)
                                                         gCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAcc                                                                                        <  1:411006/149‑1 (MQ=255)
                                                         gCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAcc                                                                                        >  1:223258/1‑149 (MQ=255)
                                                           atatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCgg                                                                                      <  3:177276/149‑1 (MQ=255)
                                                            tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGc                                                                                     <  2:394094/149‑1 (MQ=255)
                                                                gCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATaa                                                                                 <  2:315939/149‑1 (MQ=255)
                                                                gCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATaa                                                                                 >  2:744367/1‑149 (MQ=255)
                                                                  gAGCTTCAAATACATCGCTCACCAGCCCGTTCGCGTTACGCAAACGGAAGCGGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCGGACAGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAaca                                                                               <  4:286754/149‑1 (MQ=255)
                                                                      ttCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                           >  1:258471/1‑149 (MQ=255)
                                                                       tCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGc                                                                          <  2:588779/149‑1 (MQ=255)
                                                                         aaaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGa                                                                        <  3:396072/149‑1 (MQ=255)
                                                                         aaaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGa                                                                        <  1:744367/149‑1 (MQ=255)
                                                                           aTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                      <  3:71107/149‑1 (MQ=255)
                                                                               aTCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                  >  1:718461/1‑149 (MQ=255)
                                                                                 cGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACa                                                                >  4:311715/1‑149 (MQ=255)
                                                                                  gCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACaa                                                               <  4:170392/149‑1 (MQ=255)
                                                                                             cGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATa                                                    >  1:532572/1‑149 (MQ=255)
                                                                                               ttCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACg                                                  >  4:728054/1‑149 (MQ=255)
                                                                                                tCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGt                                                 <  1:548870/149‑1 (MQ=255)
                                                                                                   cTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                           <  4:675036/72‑1 (MQ=255)
                                                                                                   cTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                           >  3:675036/1‑72 (MQ=255)
                                                                                                    tttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttt                                             <  2:384360/149‑1 (MQ=255)
                                                                                                    tttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttt                                             >  4:533352/1‑149 (MQ=255)
                                                                                                      tACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGTGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttttt                                           >  3:113220/1‑149 (MQ=255)
                                                                                                         gCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                        >  3:90369/1‑149 (MQ=255)
                                                                                                         gCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                        >  4:466960/1‑149 (MQ=255)
                                                                                                           aaaCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCg                                      >  4:446319/1‑149 (MQ=255)
                                                                                                            aaCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGt                                     >  2:461324/1‑149 (MQ=255)
                                                                                                               ggAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTaa                                  <  3:327686/149‑1 (MQ=255)
                                                                                                               ggAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTaa                                  >  3:202767/1‑149 (MQ=255)
                                                                                                                 aaGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGc                                <  1:243459/149‑1 (MQ=255)
                                                                                                                  aGCAGTTATTGGCATCTATAGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGcc                               <  1:127955/149‑1 (MQ=255)
                                                                                                                   gCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCca                              >  2:6795/1‑149 (MQ=255)
                                                                                                                   gCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCca                              <  2:458716/149‑1 (MQ=255)
                                                                                                                   gCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTCAAGCca                              >  4:746683/1‑149 (MQ=255)
                                                                                                                       ttATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACAtc                          >  4:28632/1‑149 (MQ=255)
                                                                                                                           tGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                      >  1:380376/1‑149 (MQ=255)
                                                                                                                            ggCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCaga                     <  2:707276/149‑1 (MQ=255)
                                                                                                                             gCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCagag                    >  1:156692/1‑149 (MQ=255)
                                                                                                                             gCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCagag                    >  2:744441/1‑149 (MQ=255)
                                                                                                                               aTCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACGGTTAAGGCATCACAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCa                  >  4:556037/1‑149 (MQ=255)
                                                                                                                                   aTGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa              >  2:376002/1‑149 (MQ=255)
                                                                                                                                     ggtcaggatgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACg            >  2:608891/7‑149 (MQ=255)
                                                                                                                                      gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                  >  3:144590/1‑62 (MQ=255)
                                                                                                                                      gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                  <  4:144590/62‑1 (MQ=255)
                                                                                                                                       tgatgaTGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGct          >  2:439303/1‑149 (MQ=255)
                                                                                                                                         aggatgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGctct        >  4:670448/3‑149 (MQ=255)
                                                                                                                                         aggatgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGctct        >  4:275265/3‑149 (MQ=255)
                                                                                                                                          tgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTg       <  1:689483/149‑1 (MQ=255)
                                                                                                                                           gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTgg      <  3:433889/149‑1 (MQ=255)
                                                                                                                                               ccGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCaa  <  2:81516/149‑1 (MQ=255)
                                                                                                                                               ccGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCaa  <  4:338681/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAA  >  NC_000913/2429998‑2430289

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