Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,416 G→T 8.4% L272M (CTG→ATG) ‡ gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,4160GT 8.4% 248.3 / 11.4 71L272M (CTG→ATG) ‡gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base T (2/4);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

GCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATC  >  NC_000913/3963273‑3963564
                                                                                                                                               |                                                                                                                                                    
gcggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCg                                                                                                                                                 >  4:98030/1‑149 (MQ=255)
   gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                                                                                                                                              <  4:181871/149‑1 (MQ=255)
    cgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCat                                                                                                                                             <  3:182509/149‑2 (MQ=255)
       tGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttt                                                                                                                                          <  1:754552/149‑5 (MQ=255)
           tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc                                                                                                                                      >  1:281813/1‑149 (MQ=255)
           tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc                                                                                                                                      <  2:755712/149‑1 (MQ=255)
           tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc                                                                                                                                      <  2:281813/149‑1 (MQ=255)
            atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCg                                                                                                                                     >  3:285994/1‑149 (MQ=255)
                cGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa                                                                                                                                 <  3:98030/149‑1 (MQ=255)
                   aGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCattt                                                                                                                                           >  3:589591/1‑133 (MQ=255)
                   aGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCattt                                                                                                                                           <  4:589591/136‑4 (MQ=255)
                      gTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTaa                                                                                                                           <  2:619576/149‑1 (MQ=255)
                       tACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc                                                                                                                          <  4:123660/149‑1 (MQ=255)
                             tGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGt                                                                                                                    <  4:790084/149‑1 (MQ=255)
                                         gCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGc                                                                                                        >  1:729692/1‑149 (MQ=255)
                                                  ccAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                               >  1:486548/1‑149 (MQ=255)
                                                  ccAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                               >  1:343771/1‑149 (MQ=255)
                                                      aTGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCa                                                                                           <  3:528525/149‑1 (MQ=255)
                                                       tGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAg                                                                                          >  4:259106/1‑149 (MQ=255)
                                                          aaCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtc                                                                                       >  4:262266/1‑149 (MQ=255)
                                                          aaCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtc                                                                                       >  4:177465/1‑149 (MQ=255)
                                                           aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGt                                                                                                            >  3:503265/1‑127 (MQ=255)
                                                           aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGt                                                                                                            <  4:503265/127‑1 (MQ=255)
                                                           aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcc                                                                                      >  3:532980/1‑149 (MQ=255)
                                                               gCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCa                                                                                  >  1:321280/1‑149 (MQ=255)
                                                                ccGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTcc                                                                                                                 >  4:396059/1‑117 (MQ=255)
                                                                ccGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTcc                                                                                                                 <  3:396059/117‑1 (MQ=255)
                                                                     gACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGAc                                                                            >  3:479401/1‑149 (MQ=255)
                                                                        cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGc                                                                         >  4:405714/1‑149 (MQ=255)
                                                                         aGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTCCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCa                                                                                                          <  1:28808/115‑1 (MQ=255)
                                                                         aGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTCCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCa                                                                                                          >  2:28808/1‑115 (MQ=255)
                                                                           gCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaat                                                                      <  3:390447/149‑1 (MQ=255)
                                                                              ttCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaa                                                                   >  3:165720/1‑149 (MQ=255)
                                                                              ttCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaa                                                                   >  1:587400/1‑149 (MQ=255)
                                                                               tCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaat                                                                  <  4:285994/149‑1 (MQ=255)
                                                                               tCACGCAACGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaat                                                                  <  3:689409/149‑1 (MQ=255)
                                                                                     aGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCg                                                            >  2:466306/1‑149 (MQ=255)
                                                                                        gcACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTg                                                         >  2:539240/1‑149 (MQ=255)
                                                                                          aCCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGtt                                                       >  1:398726/1‑149 (MQ=255)
                                                                                                    aGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCAGGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCTCAATGAGTGGCTCGTTGTTTCAATtgctg                                             >  1:587043/1‑149 (MQ=255)
                                                                                                       gTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCaa                                          <  1:652750/149‑1 (MQ=255)
                                                                                                          aTACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                                   >  2:782514/1‑137 (MQ=255)
                                                                                                          aTACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                                   >  2:45012/1‑137 (MQ=255)
                                                                                                          aTACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                                   <  1:782514/137‑1 (MQ=255)
                                                                                                          aTACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                                   <  1:45012/137‑1 (MQ=255)
                                                                                                            acacTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgct                                              >  3:641248/1‑140 (MQ=255)
                                                                                                            acacTGAATATGCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCACCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgct                                              <  4:641248/140‑1 (MQ=255)
                                                                                                                 gAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAgg                                <  4:532980/149‑1 (MQ=255)
                                                                                                                 gAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAgg                                <  3:262266/149‑1 (MQ=255)
                                                                                                                 gAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAgg                                <  1:466306/149‑1 (MQ=255)
                                                                                                                    tatTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTa                             >  2:604823/1‑149 (MQ=255)
                                                                                                                     atTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTaa                            >  4:613709/1‑149 (MQ=255)
                                                                                                                          gTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGt                       >  4:176250/1‑149 (MQ=255)
                                                                                                                          gTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGt                       >  1:371296/1‑149 (MQ=255)
                                                                                                                            tCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTc                     <  2:343771/149‑1 (MQ=255)
                                                                                                                              ggTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcat                   >  3:668601/1‑149 (MQ=255)
                                                                                                                              ggTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcat                   <  4:479401/149‑1 (MQ=255)
                                                                                                                              ggTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcat                   <  4:406594/149‑1 (MQ=255)
                                                                                                                                tAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc                 >  3:510651/1‑149 (MQ=255)
                                                                                                                                tAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc                 >  1:394243/1‑149 (MQ=255)
                                                                                                                                 aaaaaTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGttt                                                      <  2:108850/111‑1 (MQ=255)
                                                                                                                                 aaaaaTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGttt                                                      >  1:108850/1‑111 (MQ=255)
                                                                                                                                 aaaaaTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcca                <  2:587400/149‑1 (MQ=255)
                                                                                                                                    tctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtg                                                >  4:190201/3‑114 (MQ=255)
                                                                                                                                    tctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtg                                                <  3:190201/112‑1 (MQ=255)
                                                                                                                                     aTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc            <  2:729692/149‑1 (MQ=255)
                                                                                                                                          gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg       <  4:668601/149‑1 (MQ=255)
                                                                                                                                          gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg       >  1:241627/1‑149 (MQ=255)
                                                                                                                                          gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg       >  4:469936/1‑149 (MQ=255)
                                                                                                                                           aTCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGc      >  2:21960/1‑149 (MQ=255)
                                                                                                                                              aGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCcat   >  1:477004/1‑149 (MQ=255)
                                                                                                                                               ggATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCcatc  <  3:405714/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATC  >  NC_000913/3963273‑3963564

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