Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,080 A→G 8.6% intergenic (+65/‑173) 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,0800AG 8.6% 320.3 / 16.2 104intergenic (+65/‑173)gabT/gabP4‑aminobutyrate aminotransferase, PLP‑dependent/gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base A (48/47);  new base G (3/6);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 107 total aligned reads displayed.

TTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCA  >  NC_000913/2793933‑2794225
                                                                                                                                                   |                                                                                                                                                 
ttGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                                                                                                                  >  2:461782/1‑149 (MQ=255)
ttGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                                                                                                                  >  1:782051/1‑149 (MQ=255)
ttGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                                                                                                                  >  1:307854/1‑149 (MQ=255)
 tGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAg                                                                                                                                                 <  1:616873/149‑1 (MQ=255)
 tGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAg                                                                                                                                                 <  1:295751/149‑1 (MQ=255)
   tACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGc                                                                                                                                               >  3:685890/1‑149 (MQ=255)
   tACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGc                                                                                                                                               <  1:458216/149‑1 (MQ=255)
     ccGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGt                                                                                                                                             <  3:309853/149‑1 (MQ=255)
     ccGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGt                                                                                                                                             <  2:332798/149‑1 (MQ=255)
      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                                                                                                                                    <  4:744846/149‑1 (MQ=255)
              aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                    >  3:91343/1‑149 (MQ=255)
                tGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAAcgccttat                                                                                                                                  <  1:654546/149‑5 (MQ=255)
                    gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                              <  4:351366/149‑1 (MQ=255)
                    gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                              >  4:566532/1‑149 (MQ=255)
                    gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                              >  3:867562/1‑149 (MQ=255)
                     aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                             <  2:307854/149‑1 (MQ=255)
                     aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTCTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                             <  2:264886/149‑1 (MQ=255)
                      cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                            <  2:655601/149‑2 (MQ=255)
                      cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                            <  3:59431/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                                                                                                                         <  3:96791/149‑1 (MQ=255)
                         tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                         >  1:476891/1‑149 (MQ=255)
                         tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                         >  2:202787/1‑149 (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                                                                                                                 >  1:406622/1‑140 (MQ=255)
                                          tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                                                                                                 <  2:406622/140‑1 (MQ=255)
                                          tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                                                                                                        >  2:230949/1‑149 (MQ=255)
                                          tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                                                                                                        <  3:213720/149‑1 (MQ=255)
                                             agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                     <  2:136981/149‑1 (MQ=255)
                                             agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGaaa                                                                                                     <  2:782051/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                                                                                                <  3:429214/149‑1 (MQ=255)
                                                  atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                <  2:848187/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)
                                                             ttttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt                                                                                     >  4:253638/1‑149 (MQ=255)
                                                                 gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTg                                                                                            >  4:155054/1‑138 (MQ=255)
                                                                  atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                >  1:79500/1‑149 (MQ=255)
                                                                  atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                <  1:550346/149‑1 (MQ=255)
                                                                  atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                <  3:573571/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                                                        <  4:259361/149‑1 (MQ=255)
                                                                                          ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                        <  2:116704/149‑1 (MQ=255)
                                                                                          ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                        <  2:354112/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          >  3:741777/4‑148 (MQ=255)
                                                                                                                                         attgTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGgcg          <  4:741777/145‑1 (MQ=255)
                                                                                                                                          aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                         <  3:171461/52‑1 (MQ=255)
                                                                                                                                          aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                         >  4:171461/1‑52 (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)
                                                                                                                                                   |                                                                                                                                                 
TTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCA  >  NC_000913/2793933‑2794225

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: