Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,917 A→T 5.9% intergenic (‑94/+37) yddW ← / ← gadC liprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,9170AT 5.9% 273.8 / 5.6 85intergenic (‑94/+37)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base T (2/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGA  >  NC_000913/1568769‑1569060
                                                                                                                                                    |                                                                                                                                               
cTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                 >  2:281930/1‑149 (MQ=255)
cTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                 >  2:31303/1‑149 (MQ=255)
    tATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTAt                                                                                                                                             >  4:298743/1‑149 (MQ=255)
         cTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTc                                                                                                                                        >  8:98833/1‑149 (MQ=255)
            gTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTcc                                                                                                                                       <  6:358677/147‑1 (MQ=255)
                 ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTc                                                                                                                                 >  3:142338/1‑148 (MQ=255)
                 ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTc                                                                                                                                 <  4:142338/148‑1 (MQ=255)
                 ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCt                                                                                                                                >  1:58069/1‑149 (MQ=255)
                    aTCGTTAATATCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCCTACTTATCCCCGTTTAGTCAATACGACAGCAAGCACGAAAATGGGAGCGATTAATTATCGCTCCCTTGTCtt                                                                                                                               <  8:270664/147‑1 (MQ=255)
                     tCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtat                                                                                                                             >  8:421754/1‑148 (MQ=255)
                        ttAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAAc                                                                                                                          >  5:211605/1‑148 (MQ=255)
                         tAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTTTCCCCGTTAAGTCAATACGACAGCAAGAACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCa                                                                                                                        >  4:304954/1‑149 (MQ=255)
                         tAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCGCCCTTGTCTTATAACCa                                                                                                                        >  1:402517/1‑149 (MQ=255)
                           aTTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCAtt                                                                                                                      <  8:232319/149‑1 (MQ=255)
                                     tCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCACTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATgg                                                                                                             <  7:344370/148‑1 (MQ=255)
                                            cAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGtt                                                                                                     <  7:421754/149‑1 (MQ=255)
                                               atatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTCCTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                  <  1:176431/149‑1 (MQ=255)
                                                tatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCtt                                                                                                 >  5:49906/1‑149 (MQ=255)
                                                  tCCATATGTGCTGGTTTCCTGTAAACAGATGTGCGCTCATTCACGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGTGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                <  2:279325/148‑1 (MQ=255)
                                                       atGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt                                                                                           >  3:392194/1‑148 (MQ=255)
                                                          tgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatc                                                                                       <  7:410272/149‑1 (MQ=255)
                                                            cTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                      >  3:349385/1‑148 (MQ=255)
                                                             tGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATcac                                                                                     >  7:316587/1‑148 (MQ=255)
                                                               gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTcc                                                                                                                                       <  2:274047/96‑1 (MQ=255)
                                                               gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTcc                                                                                                                                       >  1:274047/1‑96 (MQ=255)
                                                                tttCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCAGTTAAGTCAATACGACAGCAAGCACGAAAAAGTGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAata                                                                                 <  3:60279/149‑1 (MQ=255)
                                                                  tCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAata                                                                                 >  1:348059/1‑147 (MQ=255)
                                                                   ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGTGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAtata                                                                               <  1:381588/148‑1 (MQ=255)
                                                                   ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAtata                                                                               <  8:316587/148‑1 (MQ=255)
                                                                     ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                            >  2:209093/1‑149 (MQ=255)
                                                                         aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                        <  4:349385/149‑1 (MQ=255)
                                                                         aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                        >  7:33950/1‑149 (MQ=255)
                                                                                   cgcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCa                                                               <  3:166040/148‑1 (MQ=255)
                                                                                    gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc                                                                                          >  7:407707/1‑120 (MQ=255)
                                                                                    gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc                                                                                          <  8:407707/120‑1 (MQ=255)
                                                                                    gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcg                                                             >  7:211878/1‑149 (MQ=255)
                                                                                    gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcg                                                             >  4:302901/1‑149 (MQ=255)
                                                                                     cTCATTCTCGTACTTATCCCAGTTAATTCAATACGACAGAAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgc                                                            <  2:348059/149‑1 (MQ=255)
                                                                                        aTTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCggg                                                         >  7:331719/1‑149 (MQ=255)
                                                                                         ttCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGt                                                        >  3:368962/1‑149 (MQ=255)
                                                                                            tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                            >  1:124006/1‑126 (MQ=255)
                                                                                            tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                            <  2:124006/126‑1 (MQ=255)
                                                                                              gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGc                                                      <  7:98833/146‑1 (MQ=255)
                                                                                              gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCa                                                     <  3:298743/147‑1 (MQ=255)
                                                                                               tACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAgg                                                   >  7:98722/1‑148 (MQ=255)
                                                                                                 cTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGgaag                                                <  6:49906/149‑1 (MQ=255)
                                                                                                  ttATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGgaag                                                >  4:152826/1‑148 (MQ=255)
                                                                                                   tATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGgaag                                                >  7:254431/1‑147 (MQ=255)
                                                                                                    aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                   >  8:10776/1‑127 (MQ=255)
                                                                                                    aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                   <  7:10776/127‑1 (MQ=255)
                                                                                                    aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGgaagaa                                              <  4:368962/148‑1 (MQ=255)
                                                                                                       cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGa                                          <  4:392194/149‑1 (MQ=255)
                                                                                                        ccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaag                                             <  7:104703/145‑1 (MQ=255)
                                                                                                        ccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaag                                             >  8:104703/1‑145 (MQ=255)
                                                                                                        ccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGAc                                         >  8:86873/1‑149 (MQ=255)
                                                                                                           ttAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCtt                                      >  3:77846/1‑149 (MQ=255)
                                                                                                             aaGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGCCATGGGTAGGGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCttt                                     <  1:281930/148‑1 (MQ=255)
                                                                                                             aaGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCtttt                                    >  7:48843/1‑149 (MQ=255)
                                                                                                              aGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAg                                                    <  2:169803/132‑1 (MQ=255)
                                                                                                              aGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTg                                   <  8:331719/149‑1 (MQ=255)
                                                                                                                   aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                              >  1:47795/1‑149 (MQ=255)
                                                                                                                    tACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCCCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTAATTCATCACAATATAGTGTGTTTAACGTGAACGAGGGTGCAGGAAGAAGTGACATTTTGGTGCGt                             >  7:186378/1‑149 (MQ=255)
                                                                                                                      cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGtt                            <  6:211605/148‑1 (MQ=255)
                                                                                                                      cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGtt                            <  4:77846/148‑1 (MQ=255)
                                                                                                                       gACAGCAATCACGAAAAAGGGAGCGATGAATTATCGATCCCTTGTCTTATAACGATTCAGACATTGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTTGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                          <  8:186378/149‑1 (MQ=255)
                                                                                                                       gACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                          <  2:361818/117‑1 (MQ=255)
                                                                                                                       gACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                          >  1:361818/1‑117 (MQ=255)
                                                                                                                         cAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGa                        <  8:211878/149‑1 (MQ=255)
                                                                                                                              agcaCTAAAAAGGGAGCGATGAATTCTCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATCTTGTTGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGtt                   <  3:304954/149‑1 (MQ=255)
                                                                                                                                caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGa                 <  8:33950/149‑1 (MQ=255)
                                                                                                                                 aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                <  3:152826/149‑1 (MQ=255)
                                                                                                                                 aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                <  7:392833/149‑1 (MQ=255)
                                                                                                                                 aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                <  1:31303/149‑1 (MQ=255)
                                                                                                                                 aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                <  8:254431/149‑1 (MQ=255)
                                                                                                                                   gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTTTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg              >  6:96409/1‑149 (MQ=255)
                                                                                                                                    aaaagGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTgg             >  8:37976/1‑149 (MQ=255)
                                                                                                                                         gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                <  1:131722/141‑1 (MQ=255)
                                                                                                                                         gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                >  2:131722/1‑141 (MQ=255)
                                                                                                                                              cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaag                                             <  8:184062/107‑1 (MQ=255)
                                                                                                                                              cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaag                                             >  7:184062/1‑107 (MQ=255)
                                                                                                                                               gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc          >  8:151246/1‑141 (MQ=255)
                                                                                                                                               gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc          <  7:151246/141‑1 (MQ=255)
                                                                                                                                               gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa  <  3:302901/149‑1 (MQ=255)
                                                                                                                                               gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa  <  1:209093/149‑1 (MQ=255)
                                                                                                                                               gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa  <  8:98722/149‑1 (MQ=255)
                                                                                                                                                aTGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa  >  2:222145/1‑148 (MQ=255)
                                                                                                                                                aTGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa  <  1:222145/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGA  >  NC_000913/1568769‑1569060

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