Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,089 A→G 5.6% intergenic (+74/‑164) 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,0890AG 5.6% 2641.7 / 95.1 803intergenic (+74/‑164)gabT/gabP4‑aminobutyrate aminotransferase, PLP‑dependent/gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base A (406/349);  new base G (24/21);  total (432/371)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 817 total aligned reads displayed.

CTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTT  >  NC_000913/2793941‑2794237
                                                                                                                                                    |                                                                                                                                                    
cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTa                                                                                                                                                      >  3:3176917/1‑149 (MQ=255)
   aCCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTgcgc                                                                                                                                                   <  3:3834514/149‑1 (MQ=255)
   aCCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgc                                                                                                                                                   <  4:5527550/149‑1 (MQ=255)
     cATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccg                                                                                                                                                 >  3:4377119/1‑149 (MQ=255)
     cATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                                                                                                                  <  1:5439471/148‑1 (MQ=255)
      atcgAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                <  3:3637423/146‑1 (MQ=255)
      aTTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                >  2:5289076/1‑149 (MQ=255)
         gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATc                                                                                                                                             <  3:5990393/149‑1 (MQ=255)
         gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATc                                                                                                                                             <  3:5708582/149‑1 (MQ=255)
          aaGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCAccc                                                                                                                                            >  2:699194/1‑149 (MQ=255)
              cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                                        <  3:2230396/149‑2 (MQ=255)
                 tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGGTCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     >  3:6470069/1‑149 (MQ=255)
                 tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     >  4:5200771/1‑149 (MQ=255)
                 tCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     <  1:454317/149‑1 (MQ=255)
                    gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAc                                                                                                                                                     <  1:1803073/130‑1 (MQ=255)
                    gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGt                                                                                                                                  >  4:6337807/1‑149 (MQ=255)
                        cGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGt                                                                                                                              <  1:6639322/149‑1 (MQ=255)
                         gTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                                                                                                             >  3:4633360/1‑149 (MQ=255)
                          tCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGa                                                                                                                            >  4:4486687/1‑149 (MQ=255)
                             gggTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACa                                                                                                                         <  2:4190308/149‑1 (MQ=255)
                               gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                                                                                                                       >  1:5835022/1‑149 (MQ=255)
                               gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                                                                                                                       >  3:6314687/1‑149 (MQ=255)
                                 cTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                                                                                                                     <  1:2193358/149‑1 (MQ=255)
                                      gaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAAt                                                                                                                <  2:2607892/149‑1 (MQ=255)
                                         catcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                                                                                                             >  4:4906043/1‑149 (MQ=255)
                                          atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt                                                                                                            <  1:3784689/149‑1 (MQ=255)
                                            caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGt                                                                                                          <  2:2859881/149‑1 (MQ=255)
                                             aGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGtt                                                                                                         <  4:3219972/149‑1 (MQ=255)
                                              gCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTg                                                                                                        <  2:1809599/149‑1 (MQ=255)
                                                cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                          >  3:3032229/1‑113 (MQ=255)
                                                cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                          <  4:3032229/113‑1 (MQ=255)
                                                cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGc                                                                                                      >  3:2905029/1‑149 (MQ=255)
                                                cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGc                                                                                                      <  4:3486877/149‑1 (MQ=255)
                                                     ttttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTgcgc                                                                                                                                                   <  2:21279/99‑1 (MQ=255)
                                                     ttttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                                                                                                                  <  2:4678468/100‑1 (MQ=255)
                                                     ttttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGt                                                                                                                                  <  2:3857447/116‑1 (MQ=255)
                                                         gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGTGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAgg                                                                                             >  1:1411522/1‑149 (MQ=255)
                                                         gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAgg                                                                                             <  3:1132602/149‑1 (MQ=255)
                                                             aGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGa                                                                                         <  4:3621607/149‑1 (MQ=255)
                                                             aGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGa                                                                                         <  2:5835022/149‑1 (MQ=255)
                                                               gCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGc                                                                                       >  4:1831935/1‑149 (MQ=255)
                                                               gCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGc                                                                                       >  2:3114032/1‑149 (MQ=255)
                                                                 gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTAc                                                                                                                                >  3:5432879/1‑106 (MQ=255)
                                                                 gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCg                                                                                     >  1:4998524/1‑149 (MQ=255)
                                                                 gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCg                                                                                     <  1:3281389/149‑1 (MQ=255)
                                                                     cAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGctct                                                                                 >  4:6137536/1‑149 (MQ=255)
                                                                       gTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTcc                                                                               >  2:4437096/1‑149 (MQ=255)
                                                                       gTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTcc                                                                               >  3:6829562/1‑149 (MQ=255)
                                                                          gcgcCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGc                                                                            >  1:4202365/1‑149 (MQ=255)
                                                                           cgccgcTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCa                                                                           >  1:2092955/1‑149 (MQ=255)
                                                                               gcTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCaaaaa                                                                       <  2:6557354/149‑1 (MQ=255)
                                                                                cTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAg                                                                      >  2:1250710/1‑149 (MQ=255)
                                                                                 tCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAg                                                                                              >  4:5855491/1‑124 (MQ=255)
                                                                                  ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                    <  4:2754590/149‑1 (MQ=255)
                                                                                   cTATGCCGGAGGCGACGCTTCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAAt                                                                   >  3:5613852/1‑149 (MQ=255)
                                                                                    tATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                  >  4:1285117/1‑149 (MQ=255)
                                                                                    tATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATc                                                                  <  4:1953375/149‑1 (MQ=255)
                                                                                     aTGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCa                                                                 >  2:6129186/1‑149 (MQ=255)
                                                                                         cGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATaa                                                             <  1:2801863/149‑1 (MQ=255)
                                                                                          ggAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAAc                                                            <  2:779070/149‑1 (MQ=255)
                                                                                             ggCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTGCGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAAt                                                         <  3:410365/149‑1 (MQ=255)
                                                                                               cGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATta                                                       <  2:3301100/149‑1 (MQ=255)
                                                                                                gACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAAt                                                                   >  4:5640601/1‑136 (MQ=255)
                                                                                                 aCGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATtata                                                     >  4:6982948/1‑149 (MQ=255)
                                                                                                    cTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATAcgc                                                  <  2:6810388/149‑1 (MQ=255)
                                                                                                      gcgcGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGt                                                <  1:7150863/149‑1 (MQ=255)
                                                                                                           tCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGAccc                                           >  4:2265345/1‑149 (MQ=255)
                                                                                                           tCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGAccc                                           >  1:4361610/1‑149 (MQ=255)
                                                                                                               gTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTGCGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGgcg                                       >  4:1456757/1‑149 (MQ=255)
                                                                                                               gTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGgcg                                       >  1:389971/1‑149 (MQ=255)
                                                                                                                  cGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCgg                                    <  1:1012771/149‑1 (MQ=255)
                                                                                                                  cGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCgg                                    <  2:2295002/149‑1 (MQ=255)
                                                                                                                   ggCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCggg                                   >  2:713001/1‑149 (MQ=255)
                                                                                                                    gCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                                                                                                                          <  2:3547971/61‑1 (MQ=255)
                                                                                                                    gCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                                                                                                                          >  1:3547971/1‑61 (MQ=255)
                                                                                                                       tACGGGGATCCAGGTCGGATAAGGCGTTTGCGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAAt                               >  4:395265/1‑149 (MQ=255)
                                                                                                                         cGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGt                             <  2:4848471/149‑1 (MQ=255)
                                                                                                                          ggtgATCCAGGTCGGATAAGGCGTTTACGCCGCATTCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTc                            >  4:6702455/4‑149 (MQ=255)
                                                                                                                          ggggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                    >  1:946690/1‑109 (MQ=255)
                                                                                                                          ggggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGAAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTc                            >  2:5674600/1‑149 (MQ=255)
                                                                                                                           gggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCg                           <  1:1250710/149‑1 (MQ=255)
                                                                                                                             gATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCggg                         <  3:5200771/149‑1 (MQ=255)
                                                                                                                              aTCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCgggg                        <  1:7010060/149‑1 (MQ=255)
                                                                                                                                  aGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGct                    >  3:99729/1‑149 (MQ=255)
                                                                                                                                  aGGGCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGct                    >  1:6260659/1‑149 (MQ=255)
                                                                                                                                        gATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa              >  3:1977780/1‑149 (MQ=255)
                                                                                                                                         aTAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCaa             <  3:6982948/149‑1 (MQ=255)
                                                                                                                                         aTAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCaa             <  1:2991817/149‑1 (MQ=255)
                                                                                                                                          tAAGGCGTTTGCGCCGCATCCGACAATCTGCACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAg            <  1:4710420/149‑1 (MQ=255)
                                                                                                                                          tAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAg            <  4:1647347/149‑1 (MQ=255)
                                                                                                                                           aaGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGt           >  3:823633/1‑149 (MQ=255)
                                                                                                                                            aGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTg          <  3:1428698/149‑1 (MQ=255)
                                                                                                                                             ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa         <  3:1858774/149‑1 (MQ=255)
                                                                                                                                             ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa         <  3:15802/149‑1 (MQ=255)
                                                                                                                                             ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGa         <  2:1405386/149‑1 (MQ=255)
                                                                                                                                              gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCg                           >  1:2494263/1‑130 (MQ=255)
                                                                                                                                               cGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGaca       <  3:4906043/149‑1 (MQ=255)
                                                                                                                                                gTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacac      <  3:3479048/149‑1 (MQ=255)
                                                                                                                                                 tttACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacaca     <  3:2780486/149‑1 (MQ=255)
                                                                                                                                                 tttACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGacaca     <  4:823633/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTT  >  NC_000913/2793941‑2794237

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