Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,795,630 (A)6→5 100% coding (1378/1401 nt) gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,795,6250A.75.6% 424.1 / 214.4 197coding (1373/1401 nt)gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base . (71/78);  total (96/101)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 209 total aligned reads displayed.

CGGCGCAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGATATTAT  >  NC_000913/2795526‑2795723
                                                                                                   |                                                                                                  
cgGCGCAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCa                                                                                                    >  2:1133792/1‑100 (MQ=255)
 ggCGCAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGGGTATCTGTACCGTGCCAATTTTGGCTCGCTGGAAATAGCTGGTATTGTGGCa                                                                                                    >  2:2485901/1‑99 (MQ=255)
  gcgcAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCaaa                                                                                                  <  2:251752/100‑1 (MQ=255)
    gcAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCaaaaa                                                                                                <  2:2134485/100‑1 (MQ=255)
      aaCAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCaaaaaca                                                                                              >  2:4123263/1‑98 (MQ=255)
        cAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCAGATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCaaaaacacc                                                                                            <  2:3548307/100‑5 (MQ=255)
          gTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACAccc                                                                                          <  2:812112/100‑1 (MQ=255)
          gTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCg                                                                                         >  1:2312194/1‑100 (MQ=255)
           ttAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGt                                                                                        <  1:3075009/100‑1 (MQ=255)
             aGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTc                                                                                      <  2:1534780/100‑1 (MQ=255)
              gAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCa                                                                                     >  2:880278/1‑100 (MQ=255)
              gAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCa                                                                                     <  1:2351572/100‑1 (MQ=255)
                aGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACTGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCata                                                                                   <  2:2917880/100‑1 (MQ=255)
                aGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCata                                                                                   <  1:2671728/100‑1 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataa                                                                                  <  2:167647/100‑1 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataa                                                                                  >  2:2108742/1‑100 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑‑AAAACACCCGTTCataat                                                                                 <  1:3106581/100‑1 (MQ=255)
                  tGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTGGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataat                                                                                 >  2:2675056/1‑100 (MQ=255)
                  tGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataat                                                                                 >  2:1393861/1‑100 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCataat                                                                                 >  1:552403/1‑100 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataata                                                                                <  1:1748740/100‑1 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCataata                                                                                >  1:395783/1‑100 (MQ=255)
                    aTCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGACAAAACTGGTATTGTGGCAAAAAACACCCGTTCataata                                                                                >  2:1055311/1‑100 (MQ=255)
                    aTCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATAc                                                                               <  1:3524896/100‑1 (MQ=255)
                     tctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTGTTGTGGC‑AAAAACGCCCGTTCATAATAcg                                                                              >  1:4107991/1‑100 (MQ=255)
                     tctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATAcg                                                                              >  1:3297390/1‑100 (MQ=255)
                      ctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATAcgc                                                                             >  2:4432/1‑100 (MQ=255)
                        ccACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATAcgcg                                                                            <  1:1050076/100‑1 (MQ=255)
                          aCCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCt                                                                          <  1:3202616/100‑1 (MQ=255)
                           ccGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtga                                                                        <  1:3625535/100‑1 (MQ=255)
                             ggCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  2:742014/100‑1 (MQ=255)
                             ggCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  1:2494353/100‑1 (MQ=255)
                             ggCTAATTAGCGATAGGGATTATCTGTACCGTGCAAATTATGGCTCGCTGGAAAAAGCTGGTATAGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  2:2605910/100‑1 (MQ=255)
                              gCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  1:3867345/100‑1 (MQ=255)
                              gCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGAAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  1:1125560/100‑1 (MQ=255)
                               atctttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  2:581408/96‑1 (MQ=255)
                                ttattaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCGGGTGTTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatcat                                                                    >  1:3367992/1‑99 (MQ=255)
                                  attaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTAGCTGGCAAAAGCTGGTAATGTGGC‑AAAAACACCCGTTCATAATACGCGCtgagcatgat                                                                 >  2:3497566/1‑100 (MQ=255)
                                   ttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCtgatcaagat                                                                 >  1:581412/1‑100 (MQ=255)
                                   ttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatcatgatc                                                                <  2:3297367/100‑1 (MQ=255)
                                    tccgatctgggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCtgatcatgatca                                                               <  2:1882887/92‑1 (MQ=255)
                                     aGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAg                                                              >  2:4009236/1‑100 (MQ=255)
                                      gggATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCT‑AAAAAAAAAAAAAAAACACCCGTTCATAATACGCGCtgatcatga                                                                  <  2:165299/93‑1 (MQ=255)
                                      gCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGATAAAGCTGGTATTGTGGC‑AAAAACACCCTTTCATAATACGCGATTATCATGATCAg                                                              >  2:328595/1‑99 (MQ=255)
                                      gCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAgg                                                             <  1:1245505/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  1:48959/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  1:769103/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  2:2281076/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           >  2:1274266/1‑100 (MQ=255)
                                         aTAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGt                                                          >  2:3696000/1‑100 (MQ=255)
                                          tAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTc                                                         <  1:97093/100‑1 (MQ=255)
                                            gggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcac                                                       >  1:1176254/1‑100 (MQ=255)
                                            gggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcac                                                       >  2:845843/1‑100 (MQ=255)
                                             ggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAgcgcgtatt                                                      <  1:3578503/100‑9 (MQ=255)
                                             agATTATCTGTACCGTGCCAAGTAGAGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  2:1822576/94‑1 (MQ=255)
                                              gATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcaca                                                      >  1:1569279/1‑100 (MQ=255)
                                                  aTCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCAt                                                  >  2:3462229/1‑100 (MQ=255)
                                                  aTCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATg                                                 <  1:3278894/100‑1 (MQ=255)
                                                   tCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATg                                                 <  2:3035642/100‑1 (MQ=255)
                                                    cTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGAc                                               >  2:322580/1‑100 (MQ=255)
                                                       tACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAt                                            >  2:1983501/1‑100 (MQ=255)
                                                        aCCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAtt                                           <  1:1974276/100‑1 (MQ=255)
                                                        aCCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGG‑GAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAtt                                           <  2:4028753/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          <  2:2820417/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          <  1:954126/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          >  1:3237249/1‑100 (MQ=255)
                                                          cGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTAc                                         >  1:2867426/1‑100 (MQ=255)
                                                             gCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGtc                                      >  1:1435219/1‑100 (MQ=255)
                                                              ccAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGtc                                      <  1:2935325/100‑1 (MQ=255)
                                                                aaTTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTctct                                   <  1:2145007/100‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTTTGGC‑AAAAACTCCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTg                                  <  2:1336393/98‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGa                                <  2:1435210/100‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGa                                >  1:2888872/1‑100 (MQ=255)
                                                                     tGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATg                              <  2:2958753/100‑1 (MQ=255)
                                                                     tGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATg                              <  2:2321964/100‑1 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATgg                             >  1:238268/1‑100 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGc                            >  2:721698/1‑100 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGc                            <  1:4019649/100‑1 (MQ=255)
                                                                        cTCGCTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCt                           >  1:16917/1‑100 (MQ=255)
                                                                            cTGGAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCg                       >  1:4047770/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgc                      >  1:80425/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgct                     >  1:3861470/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgct                     >  2:3370553/1‑100 (MQ=255)
                                                                                 aaaaGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctgg                   >  1:3664767/1‑100 (MQ=255)
                                                                                 aaaaGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggc                  <  1:702525/100‑1 (MQ=255)
                                                                                   aaGCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggctg                >  1:2015004/1‑100 (MQ=255)
                                                                                    aGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggctg                <  1:1188148/100‑1 (MQ=255)
                                                                                     gCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggcttan              <  2:3169116/100‑4 (MQ=255)
                                                                                     gCTGGTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTgaa              <  2:1723532/100‑1 (MQ=255)
                                                                                       tGGTATTGTGGC‑AAAAACCCCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTgaaga            <  1:80195/100‑1 (MQ=255)
                                                                                         gTATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAAc          >  2:3776313/1‑100 (MQ=255)
                                                                                          tATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACg         <  2:2072443/100‑1 (MQ=255)
                                                                                          tATTGTGGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACg         >  1:1506412/1‑100 (MQ=255)
                                                                                            ttGTGGC‑AAAAACACCCGTTCATAATACGGGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGat       <  2:660103/100‑1 (MQ=255)
                                                                                             tgtgGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGat       >  1:1513476/1‑100 (MQ=255)
                                                                                              gtgGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  2:3432750/1‑100 (MQ=255)
                                                                                              gtgGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  1:1708107/1‑100 (MQ=255)
                                                                                              gtgGC‑AAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  1:3925370/1‑100 (MQ=255)
                                                                                               gtggaaaacaCACCCGTTCATAATACGCGCTGATCATGATCAGTAGTCACACCATTACCATGACGTGTCTGGATGGCTATCTCTGGCGGAAGAACGacat     >  2:1802943/5‑100 (MQ=255)
                                                                                                 gCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGAtatta   <  2:2945672/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGGCGCAACAGTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGATATTAT  >  NC_000913/2795526‑2795723

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