Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,124 G→A 6.0% G465D (GGC→GAC)  yjbB → putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,229,1240GA 6.0% 257.3 / 5.9 83G465D (GGC→GAC) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base A (3/2);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

GATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTGGGA  >  NC_000913/4228976‑4229272
                                                                                                                                                    |                                                                                                                                                    
gATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                                                                                                                                                      <  4:665038/149‑1 (MQ=255)
    tcggtgAGCGCATGGGCAGCGAAATTGCTGACAAATCACGGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAATGGATGCGCTCTATGTGCAATCGCTCAGTAATTTAAAGCTGGCAATGTCGGTGGTCTTCTCTGGCGAt                                                                                                                                                  <  2:840359/144‑1 (MQ=255)
       tcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                                                                                                                                               <  2:239438/149‑1 (MQ=255)
         gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                                                                                                                                             >  2:80827/1‑149 (MQ=255)
         gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                                                                                                                                             >  1:115866/1‑149 (MQ=255)
          aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAcc                                                                                                                                            >  1:816951/1‑149 (MQ=255)
                  ggCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                                                                                                                                    <  2:587881/149‑1 (MQ=255)
                     aGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtcg                                                                                                                                 >  2:148785/1‑149 (MQ=255)
                     aGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtcg                                                                                                                                 >  1:23512/1‑149 (MQ=255)
                       cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGtt                                                                                                                               >  4:448932/1‑149 (MQ=255)
                        gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                                                                                                                              <  1:285920/149‑1 (MQ=255)
                            ttGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGt                                                                                                                          <  3:85643/149‑1 (MQ=255)
                                tGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgca                                                                                                                      <  4:566552/149‑1 (MQ=255)
                                tGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgca                                                                                                                      <  1:307650/149‑1 (MQ=255)
                                 gACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcag                                                                                                                     >  1:132453/1‑149 (MQ=255)
                                      aTCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAc                                                                                                                >  3:416430/1‑149 (MQ=255)
                                       tCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACa                                                                                                               >  1:147022/1‑149 (MQ=255)
                                          cTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCg                                                                                                            >  1:771236/1‑149 (MQ=255)
                                           tGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCTGCGACGCAAAcgacga                                                                                                           <  1:557289/149‑6 (MQ=255)
                                               agcaCGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACCTCGTTTTc                                                                                                       <  3:814484/149‑1 (MQ=255)
                                                    ggcggGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCAtt                                                                                                  <  3:448932/149‑1 (MQ=255)
                                                     gcggGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTc                                                                                                 <  4:457923/149‑1 (MQ=255)
                                                     gcggGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTc                                                                                                 <  2:307863/149‑1 (MQ=255)
                                                          cATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAAt                                                                                            <  4:284866/149‑1 (MQ=255)
                                                          cATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAAt                                                                                            <  3:186410/149‑1 (MQ=255)
                                                             tttCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATcgc                                                                                         <  1:148797/149‑1 (MQ=255)
                                                                  cTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTa                                                                                    >  2:84364/1‑149 (MQ=255)
                                                                  cTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTa                                                                                    >  4:521026/1‑149 (MQ=255)
                                                                       tGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTccc                                                                               >  3:327442/1‑149 (MQ=255)
                                                                             gAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTcccacgccc                                                                         >  3:578508/1‑149 (MQ=255)
                                                                               aGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTcccacgcccac                                                                       >  3:230840/1‑149 (MQ=255)
                                                                               aGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTcccacgcccac                                                                       >  1:714721/1‑149 (MQ=255)
                                                                                ggAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTc                                                                                                       <  4:738618/116‑1 (MQ=255)
                                                                                ggAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTc                                                                                                       >  3:738618/1‑116 (MQ=255)
                                                                                ggAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTcccacgcccacg                                                                      >  2:390408/1‑149 (MQ=255)
                                                                                      ggATGCGCTCTATGAGCAAGTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATc                                                                >  3:517189/1‑149 (MQ=255)
                                                                                       gATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATACCCACGCCCACGTCGATCg                                                               >  4:601390/1‑149 (MQ=255)
                                                                                         tGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGcc                                                             >  2:10897/1‑149 (MQ=255)
                                                                                            gcTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGc                                                          >  1:656941/1‑149 (MQ=255)
                                                                                            gcTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGc                                                          >  2:529896/1‑149 (MQ=255)
                                                                                              tctATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCACCGCTATTCCCACGCCCACGTCGATCGCCTGCAt                                                        <  3:588481/149‑1 (MQ=255)
                                                                                                tATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCa                                                      <  1:231184/149‑1 (MQ=255)
                                                                                                tATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCa                                                      >  1:609372/1‑149 (MQ=255)
                                                                                                tATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGGGTTGTTCTCTGGCGATGCCACCGGCGCTCGTCGTTTGAGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCTATCGCCTACATct                                                      >  2:498756/1‑148 (MQ=255)
                                                                                                  tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGc                                                    >  3:690770/1‑149 (MQ=255)
                                                                                                  tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGc                                                    >  1:845523/1‑149 (MQ=255)
                                                                                                   gAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCa                                                   >  3:525099/1‑149 (MQ=255)
                                                                                                   gAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTACGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCa                                                   >  4:539453/1‑149 (MQ=255)
                                                                                                    aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCaa                                                  >  4:484351/1‑149 (MQ=255)
                                                                                                       aaTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAAc                                               <  1:80827/149‑1 (MQ=255)
                                                                                                           gCTCAGTAATTTAAAGCTGGTAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCACCGCTATTCCCGCGCCCACGTCGATCGCCTGCATCAGCAAAACGTGc                                           <  1:473972/149‑1 (MQ=255)
                                                                                                            cTCAGTAATTTAAAGCTGGCAATGTCGGTGTTATTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCa                                          >  3:135862/1‑149 (MQ=255)
                                                                                                             tCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCag                                         <  3:103846/149‑1 (MQ=255)
                                                                                                             tCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGACGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCag                                         >  1:675039/1‑149 (MQ=255)
                                                                                                               aGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCagag                                       >  1:701414/1‑149 (MQ=255)
                                                                                                               aGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCagag                                       >  1:817932/1‑149 (MQ=255)
                                                                                                                gTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGc                                      >  1:766808/1‑149 (MQ=255)
                                                                                                                 tAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCGAAACGTGCAGAGCa                                     >  4:260923/1‑149 (MQ=255)
                                                                                                                    tttAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTg                                  <  3:12392/149‑1 (MQ=255)
                                                                                                                       aaaGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAGACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGaaa                               <  4:499368/149‑1 (MQ=255)
                                                                                                                         aGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCATGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAAcc                             <  3:539453/149‑1 (MQ=255)
                                                                                                                         aGCTGGCAATGGCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTCTTCGCATTCTTAATCGCCGATATTCCCACGCCCCCGTCGATGGCCTGCATCAGCAAAACGTGCAGAGCATCGAAGcc                             <  2:634269/149‑1 (MQ=255)
                                                                                                                               cAATGTCGGTGTCCTTCCCTGGCGATGTCATCAGCGCTCGTCGTTGGCGTCGCAGCAAACATCGTTTTCGCAGTCTTAATTGCCGCTATTCCCGCTCCCACTTCGATCGCCTGCAGCAGCAAAACGTGCAGAGCATTGAAACCAGTTcc                       <  2:771236/149‑2 (MQ=255)
                                                                                                                                aaTGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGc                      >  3:147640/1‑149 (MQ=255)
                                                                                                                                 aTGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCt                     >  3:389277/1‑149 (MQ=255)
                                                                                                                                      ggTGTTCTTCTCTGACGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCt                                                            <  4:163759/105‑1 (MQ=255)
                                                                                                                                      ggTGTTCTTCTCTGACGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCt                                                            <  2:515590/105‑1 (MQ=255)
                                                                                                                                      ggTGTTCTTCTCTGACGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCt                                                            >  3:163759/1‑105 (MQ=255)
                                                                                                                                      ggTGTTCTTCTCTGACGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCt                                                            >  1:515590/1‑105 (MQ=255)
                                                                                                                                       gtgtTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACAttt               >  4:245030/1‑149 (MQ=255)
                                                                                                                                         gtTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAg             <  1:822118/149‑1 (MQ=255)
                                                                                                                                         gtTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAg             <  1:529896/149‑1 (MQ=255)
                                                                                                                                          ttcttcTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACAt                 >  2:831081/1‑144 (MQ=255)
                                                                                                                                          ttcttcTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACAt                 <  1:831081/144‑1 (MQ=255)
                                                                                                                                          ttcttcTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAgg            >  4:259474/1‑149 (MQ=255)
                                                                                                                                           tcttctCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGc           <  1:808483/149‑1 (MQ=255)
                                                                                                                                             ttctCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCAtt                                   <  3:638590/123‑1 (MQ=255)
                                                                                                                                             ttctCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCAtt                                   >  4:638590/1‑123 (MQ=255)
                                                                                                                                               ctctGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTAc       <  1:603055/149‑1 (MQ=255)
                                                                                                                                                tctGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACt      >  4:835764/1‑149 (MQ=255)
                                                                                                                                                 ctGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGCGCATTGAAACCAGTTCGCTACATTTAGGCTTACTg     >  4:168602/1‑149 (MQ=255)
                                                                                                                                                  tGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTgg    <  4:636477/149‑1 (MQ=255)
                                                                                                                                                  tGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTgg    <  2:180865/149‑1 (MQ=255)
                                                                                                                                                  tGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTgg    >  1:95151/1‑149 (MQ=255)
                                                                                                                                                   ggCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTggg   <  4:554809/149‑1 (MQ=255)
                                                                                                                                                    gCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTGGga  <  4:422967/149‑1 (MQ=255)
                                                                                                                                                    gCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTGGga  >  3:10621/1‑149 (MQ=255)
                                                                                                                                                    gCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTGGga  <  2:609372/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAACATCGTTTTCGCATTCTTAATCGCCGCTATTCCCACGCCCACGTCGATCGCCTGCATCAGCAAAACGTGCAGAGCATTGAAACCAGTTCGCTACATTTAGGCTTACTGGGA  >  NC_000913/4228976‑4229272

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