Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,087 T→A 8.6% intergenic (+72/‑166) gabT → / → gabP 4‑aminobutyrate aminotransferase, PLP‑dependent/gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,0870TA 8.6% 346.9 / 18.7 105intergenic (+72/‑166)gabT/gabP4‑aminobutyrate aminotransferase, PLP‑dependent/gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base T (50/46);  new base A (3/6);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 110 total aligned reads displayed.

CGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACAC  >  NC_000913/2793939‑2794235
                                                                                                                                                    |                                                                                                                                                    
cGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGtt                                                                                                                                                      >  1:348334/1‑149 (MQ=255)
cGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGtt                                                                                                                                                      >  1:27701/1‑149 (MQ=255)
 gCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGttt                                                                                                                                                     <  2:311746/149‑1 (MQ=255)
   tCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAc                                                                                                                                                   <  3:467625/149‑1 (MQ=255)
     aCCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgc                                                                                                                                                 <  2:774687/149‑1 (MQ=255)
     aCCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAagc                                                                                                                                                 >  3:585883/1‑149 (MQ=255)
       cATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccg                                                                                                                                               >  3:878767/1‑149 (MQ=255)
        aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                              >  3:91343/1‑149 (MQ=255)
        aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                              <  4:744846/149‑1 (MQ=255)
          tGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAAcgccttat                                                                                                                                            <  1:654546/149‑5 (MQ=255)
            aaGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                          >  1:528296/1‑149 (MQ=255)
              gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                        <  4:351366/149‑1 (MQ=255)
              gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                        >  3:867562/1‑149 (MQ=255)
               aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                       <  2:307854/149‑1 (MQ=255)
               aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTCTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                       <  2:264886/149‑1 (MQ=255)
                cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                                      <  3:59431/149‑2 (MQ=255)
                cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                                      <  2:655601/149‑2 (MQ=255)
                 gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACaa                                                                                                                                     >  1:653618/1‑149 (MQ=255)
                 gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACaa                                                                                                                                     >  2:822832/1‑149 (MQ=255)
                 gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACaa                                                                                                                                     <  4:867562/149‑1 (MQ=255)
                   tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                   >  1:476891/1‑149 (MQ=255)
                   tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                   >  2:370500/1‑149 (MQ=255)
                   tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                   <  3:96791/149‑1 (MQ=255)
                   tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                   >  2:202787/1‑149 (MQ=255)
                       aTCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                                                                                                               >  4:695921/1‑149 (MQ=255)
                         ccGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCGTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACg                                                                                                                             >  2:751664/1‑149 (MQ=255)
                          cGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGt                                                                                                                            <  4:42602/149‑1 (MQ=255)
                               gggTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACa                                                                                                                       >  1:836184/1‑149 (MQ=255)
                                   cTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                                   <  2:836184/149‑1 (MQ=255)
                                   cTGGAGATCAGCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTACGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGCATTCAGGTCAGATAAGGCGTTTACGCCGCATCCGGCAATCTGTACGTGAACAGgaa                                                                                                                   <  2:862357/149‑1 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                                                                                                           <  2:406622/140‑1 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                                                                                                           >  1:406622/1‑140 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                                                                                                                  <  3:213720/149‑1 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                                                                                                                  >  2:230949/1‑149 (MQ=255)
                                       agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                               <  2:136981/149‑1 (MQ=255)
                                       agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                               >  1:706849/1‑149 (MQ=255)
                                         atcatcAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATc                                                                                                             >  2:872851/1‑149 (MQ=255)
                                            atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAAcgcctta                                                                                                                                             <  4:842074/114‑4 (MQ=255)
                                            atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAAcgcctta                                                                                                                                             >  3:842074/1‑111 (MQ=255)
                                            atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                          <  2:848187/149‑1 (MQ=255)
                                            atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                          <  3:429214/149‑1 (MQ=255)
                                            atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                          <  4:259574/149‑1 (MQ=255)
                                             tcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATg                                                                                                         <  1:503042/149‑1 (MQ=255)
                                              caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGt                                                                                                        >  4:336162/1‑149 (MQ=255)
                                                   aGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAAt                                                                                                              >  1:340568/1‑138 (MQ=255)
                                                   aGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAAt                                                                                                              <  2:340568/138‑1 (MQ=255)
                                                     tgtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCgg                                                                                                 >  1:548584/1‑149 (MQ=255)
                                                      gtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCggg                                                                                                <  1:872851/149‑1 (MQ=255)
                                                       ttttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt                                                                                               >  4:253638/1‑149 (MQ=255)
                                                           gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTg                                                                                                      >  4:155054/1‑138 (MQ=255)
                                                           gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTg                                                                                                      <  3:155054/138‑1 (MQ=255)
                                                            atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                          <  3:573571/149‑1 (MQ=255)
                                                            atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                          <  1:550346/149‑1 (MQ=255)
                                                            atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                          >  1:79500/1‑149 (MQ=255)
                                                                ggCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAg                                                                                      >  1:502491/1‑149 (MQ=255)
                                                                 gCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGc                                                                                     >  2:836558/1‑149 (MQ=255)
                                                                   gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATATGTCCGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCg                                                                                   >  3:94301/1‑149 (MQ=255)
                                                                      gCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGctc                                                                                <  1:230949/149‑1 (MQ=255)
                                                                            gcgcCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGc                                                                          >  4:46082/1‑149 (MQ=255)
                                                                                   tCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGa                                                                   >  2:275178/1‑149 (MQ=255)
                                                                                    ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                  <  4:259361/149‑1 (MQ=255)
                                                                                    ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                  <  2:116704/149‑1 (MQ=255)
                                                                                    ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                  <  2:354112/149‑1 (MQ=255)
                                                                                     cTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAAt                                                                 <  3:703664/149‑1 (MQ=255)
                                                                                     cTACGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTGCGGGGATCGAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGCACGTGAACAGGAAGAAAACTATGTTGGCCGGGTAAGGCGGAGCCGCTTTCCGGCAAAAAGAAt                                                                 <  2:383394/149‑1 (MQ=255)
                                                                                          ccGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATa                                                            <  1:61283/149‑1 (MQ=255)
                                                                                          ccGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATa                                                            >  3:533717/1‑149 (MQ=255)
                                                                                           cGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATaa                                                           <  2:348334/149‑1 (MQ=255)
                                                                                              aGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACaa                                                        <  4:878767/149‑1 (MQ=255)
                                                                                              aGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACaa                                                        <  1:244582/149‑1 (MQ=255)
                                                                                                 cGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATta                                                     <  3:253638/149‑1 (MQ=255)
                                                                                                  gACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATtat                                                    >  3:879138/1‑149 (MQ=255)
                                                                                                    cGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATAc                                                  >  4:661176/1‑149 (MQ=255)
                                                                                                             tCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGAccc                                         >  3:42947/1‑149 (MQ=255)
                                                                                                                  tCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                >  1:743030/1‑121 (MQ=255)
                                                                                                                  tCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                <  2:743030/121‑1 (MQ=255)
                                                                                                                            ggggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTc                          >  1:77574/1‑149 (MQ=255)
                                                                                                                              ggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCgg                        <  3:695921/149‑1 (MQ=255)
                                                                                                                                aTCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCgggg                      >  1:21157/1‑149 (MQ=255)
                                                                                                                                   cAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGgcgc                   <  1:369920/149‑1 (MQ=255)
                                                                                                                                   attgTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGgcg                    >  3:741777/4‑148 (MQ=255)
                                                                                                                                   attgTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGgcg                    <  4:741777/145‑1 (MQ=255)
                                                                                                                                    aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                                   >  4:171461/1‑52 (MQ=255)
                                                                                                                                    aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                                   <  3:171461/52‑1 (MQ=255)
                                                                                                                                    aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGct                  >  2:534904/1‑149 (MQ=255)
                                                                                                                                       tCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCtctc               >  1:700602/1‑149 (MQ=255)
                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa            >  4:626863/1‑149 (MQ=255)
                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa            >  4:367680/1‑149 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTg        <  1:411838/149‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa       >  4:568308/1‑149 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa       <  2:27701/149‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCACAAGTGa       <  2:476891/149‑1 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCtctc               <  4:181570/140‑1 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCtctc               >  3:181570/1‑140 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGac      >  1:460613/1‑149 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGac      >  1:709939/1‑149 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGaca     <  4:565243/149‑1 (MQ=255)
                                                                                                                                                  gTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAGTAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacac    >  4:758913/1‑149 (MQ=255)
                                                                                                                                                   tttACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacaca   >  1:475094/1‑149 (MQ=255)
                                                                                                                                                    ttACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacacac  <  2:653618/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACAC  >  NC_000913/2793939‑2794235

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