Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,159,182 C→T A438A (GCC→GCT ptsG → glucose‑specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,1820CT97.0% 220.0 / ‑1.9 66A438A (GCC→GCTptsGglucose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (34/30);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGG  >  NC_000913/1159046‑1159318
                                                                                                                                        |                                                                                                                                          
cAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAggctggc                                                                                                                                         <  1:320537/140‑1 (MQ=255)
      cGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTgaag                                                                                                                                   >  1:892988/1‑140 (MQ=255)
          aGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAAc                                                                                                                               <  2:944514/140‑1 (MQ=255)
               ccAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGgc                                                                                                                          >  1:325499/1‑140 (MQ=255)
                       aTGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCggg                                                                                                                  >  2:760830/1‑140 (MQ=255)
                                   ctGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGc                                                                                                      >  2:700419/1‑140 (MQ=255)
                                         gctgcATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTc                                                                                                <  1:385697/140‑1 (MQ=255)
                                                                  ttACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACt                                                                     <  2:892988/140‑1 (MQ=255)
                                                                   tACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTa                                                                    <  2:845486/140‑1 (MQ=255)
                                                                   tACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTa                                                                    <  2:397695/140‑1 (MQ=255)
                                                                           tCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGATGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGAt                                                            >  2:798893/1‑140 (MQ=255)
                                                                             gACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATaa                                                          >  1:883236/1‑140 (MQ=255)
                                                                             gACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATaa                                                          >  2:479861/1‑140 (MQ=255)
                                                                             gACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATaa                                                          >  2:317369/1‑140 (MQ=255)
                                                                                  aTGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGa                                                     >  2:310123/1‑140 (MQ=255)
                                                                                    gTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGaaa                                                   >  1:937582/1‑140 (MQ=255)
                                                                                       ttACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAAcc                                                >  1:42023/1‑140 (MQ=255)
                                                                                        tACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCg                                               <  2:937582/140‑1 (MQ=255)
                                                                                        tACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCg                                               <  1:563032/140‑1 (MQ=255)
                                                                                         aCCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCga                                              <  1:160953/140‑1 (MQ=255)
                                                                                         aCCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCga                                              >  2:704771/1‑140 (MQ=255)
                                                                                           ccGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCgaga                                            <  1:986785/140‑1 (MQ=255)
                                                                                             gTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatg                                          >  2:831738/1‑140 (MQ=255)
                                                                                             gTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatg                                          >  1:86573/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgat                                           >  2:900247/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCga                                              <  1:566731/135‑1 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCga                                              >  2:566731/1‑135 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  2:428849/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:908265/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:998680/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         <  2:785532/140‑1 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:591936/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         <  2:998904/140‑1 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:435052/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:793437/1‑140 (MQ=255)
                                                                                              tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgg                                         >  1:913992/1‑140 (MQ=255)
                                                                                               cTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgga                                        >  2:170351/1‑140 (MQ=255)
                                                                                               cTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatgga                                        >  1:239185/1‑140 (MQ=255)
                                                                                                tGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCgaga                                            <  2:744246/135‑1 (MQ=255)
                                                                                                tGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCgaga                                            >  1:744246/1‑135 (MQ=255)
                                                                                                tGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatggat                                       <  2:903294/140‑1 (MQ=255)
                                                                                                tGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatggat                                       <  2:998680/140‑1 (MQ=255)
                                                                                                tGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAgatggat                                       <  2:591936/140‑1 (MQ=255)
                                                                                                    cgTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGt                                   >  1:1112647/1‑140 (MQ=255)
                                                                                                     gTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTa                                  <  2:913992/140‑1 (MQ=255)
                                                                                                     gTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTa                                  <  2:1112647/140‑1 (MQ=255)
                                                                                                      tCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGt                                   <  1:900247/140‑1 (MQ=255)
                                                                                                        aGCGTTGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACAt                               >  1:1163439/1‑140 (MQ=255)
                                                                                                            ttGCTGATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGt                           <  2:632024/140‑1 (MQ=255)
                                                                                                                 gATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCa                      <  1:170351/140‑1 (MQ=255)
                                                                                                                 gATGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCa                      <  2:883236/140‑1 (MQ=255)
                                                                                                                  aTGTGTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCAc                     <  1:428849/140‑1 (MQ=255)
                                                                                                                   tgtgTCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACt                    <  1:798893/140‑1 (MQ=255)
                                                                                                                    gtgtCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTa                   <  2:155388/140‑1 (MQ=255)
                                                                                                                     tgtCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTaa                  <  2:1163439/140‑1 (MQ=255)
                                                                                                                       tCTAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTaacc                >  1:314018/1‑138 (MQ=255)
                                                                                                                         tAAAGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCg              >  2:91826/1‑140 (MQ=255)
                                                                                                                           aaGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTa            <  2:42023/140‑1 (MQ=255)
                                                                                                                           aaGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTa            <  1:831738/140‑1 (MQ=255)
                                                                                                                           aaGTGGATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTa            <  1:1066218/140‑1 (MQ=255)
                                                                                                                                gATCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACg       <  2:908265/140‑1 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGt      >  2:1037569/1‑140 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGt      >  2:910144/1‑140 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGt      >  1:159797/1‑140 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCAGTAACCa                      >  2:519218/1‑124 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCAGTAACCa                      <  1:519218/124‑1 (MQ=255)
                                                                                                                                 aTCAGGCTGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGACGGATGAGTACATCCGTAACCACTAATCCGTAAGACGt      >  2:741317/1‑140 (MQ=255)
                                                                                                                                     ggctggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTggg  <  2:187390/140‑1 (MQ=255)
                                                                                                                                     ggctggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTggg  <  1:1037569/140‑1 (MQ=255)
                                                                                                                                     ggctggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTggg  <  1:91826/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
CAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTG‑‑GTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGG  >  NC_000913/1159046‑1159318

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