Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,094 G→T 6.1% intergenic (+79/‑159) 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,0940GT 6.1% 342.4 / 10.0 98intergenic (+79/‑159)gabT/gabP4‑aminobutyrate aminotransferase, PLP‑dependent/gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base G (48/44);  new base T (2/4);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 102 total aligned reads displayed.

CATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTTTCG  >  NC_000913/2793946‑2794240
                                                                                                                                                    |                                                                                                                                                  
cATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccg                                                                                                                                                    >  3:878767/1‑149 (MQ=255)
 aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                   <  4:744846/149‑1 (MQ=255)
 aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                   >  3:91343/1‑149 (MQ=255)
     aaGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                               >  1:528296/1‑149 (MQ=255)
       gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                             <  4:351366/149‑1 (MQ=255)
       gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                             >  3:867562/1‑149 (MQ=255)
       gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                             >  4:566532/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                                                                                                                                          <  4:867562/149‑1 (MQ=255)
          gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACaa                                                                                                                                          >  1:653618/1‑149 (MQ=255)
            tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                        >  2:202787/1‑149 (MQ=255)
            tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                        >  2:370500/1‑149 (MQ=255)
            tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                        >  1:476891/1‑149 (MQ=255)
            tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                        <  3:96791/149‑1 (MQ=255)
                aTCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                                                                                                                    >  4:695921/1‑149 (MQ=255)
                  ccGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCGTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACg                                                                                                                                  >  2:751664/1‑149 (MQ=255)
                  ccGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACg                                                                                                                                  >  1:699105/1‑149 (MQ=255)
                   cGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGt                                                                                                                                 <  4:42602/149‑1 (MQ=255)
                       aGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                                                                                                                             <  4:178981/149‑1 (MQ=255)
                       aGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                                                                                                                             >  3:237147/1‑149 (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                                                                                                                    >  1:706849/1‑149 (MQ=255)
                                agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                                    <  2:136981/149‑1 (MQ=255)
                                agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                                    <  2:782051/149‑1 (MQ=255)
                                  atcatcAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATc                                                                                                                  >  2:872851/1‑149 (MQ=255)
                                     atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                               <  3:429214/149‑1 (MQ=255)
                                     atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                               <  2:848187/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                                                                                                                   <  2:340568/138‑1 (MQ=255)
                                            aGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAAt                                                                                                                   >  1:340568/1‑138 (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:79500/1‑149 (MQ=255)
                                                     atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                               <  1:550346/149‑1 (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)
                                                                  gTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTcc                                                                                  >  2:369920/1‑149 (MQ=255)
                                                                     gcgcCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGc                                                                               >  4:46082/1‑149 (MQ=255)
                                                                            tCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGa                                                                        >  2:275178/1‑149 (MQ=255)
                                                                             ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                       <  2:116704/149‑1 (MQ=255)
                                                                             ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                       <  4:259361/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                                                                <  1:461782/149‑1 (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)
                                                                                           gACGCTGCGCGTCTTGTCCGGTCTACGGGGATCCAGGTCGGGTAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATtat                                                         >  2:150283/1‑149 (MQ=255)
                                                                                           gACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATtat                                                         >  3:879138/1‑149 (MQ=255)
                                                                                             cGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATAc                                                       >  4:661176/1‑149 (MQ=255)
                                                                                                      tCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGAccc                                              >  3:42947/1‑149 (MQ=255)
                                                                                                           tCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                     <  2:743030/121‑1 (MQ=255)
                                                                                                           tCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                     >  1:743030/1‑121 (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                         <  4:741777/145‑1 (MQ=255)
                                                                                                                            attgTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGgcg                         >  3:741777/4‑148 (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)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGTGAAATGTCGGGGCGCTCTCCCa                 >  3:728700/1‑149 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa                 >  4:367680/1‑149 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa                 >  4:626863/1‑149 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTg             <  1:411838/149‑1 (MQ=255)
                                                                                                                                        ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa            <  2:27701/149‑1 (MQ=255)
                                                                                                                                        ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa            >  4:568308/1‑149 (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)
                                                                                                                                                  ccgcATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTTTCg  <  4:334458/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTTTCG  >  NC_000913/2793946‑2794240

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: