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)
          gTTAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCg                                                                                         >  1:2312194/1‑100 (MQ=255)
           ttAGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGt                                                                                        <  1:3075009/100‑1 (MQ=255)
             aGAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTc                                                                                      <  2:1534780/100‑1 (MQ=255)
              gAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCa                                                                                     >  2:880278/1‑100 (MQ=255)
              gAAGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCa                                                                                     <  1:2351572/100‑1 (MQ=255)
                aGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACTGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCata                                                                                   <  2:2917880/100‑1 (MQ=255)
                aGTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCata                                                                                   <  1:2671728/100‑1 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataa                                                                                  <  2:167647/100‑1 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataa                                                                                  >  2:2108742/1‑100 (MQ=255)
                 gTGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGC‑‑AAAACACCCGTTCataat                                                                                 <  1:3106581/100‑1 (MQ=255)
                  tGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTGGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataat                                                                                 >  2:2675056/1‑100 (MQ=255)
                  tGATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataat                                                                                 >  2:1393861/1‑100 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCataat                                                                                 >  1:552403/1‑100 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataata                                                                                <  1:1748740/100‑1 (MQ=255)
                   gATCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCataata                                                                                >  1:395783/1‑100 (MQ=255)
                    aTCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGACAAAACTGGTATTGTGGCAAAAAACACCCGTTCataata                                                                                >  2:1055311/1‑100 (MQ=255)
                    aTCTCCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATAc                                                                               <  1:3524896/100‑1 (MQ=255)
                     tctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTGTTGTGGCAAAAACGCCCGTTCATAATAcg                                                                              >  1:4107991/1‑100 (MQ=255)
                     tctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATAcg                                                                              >  1:3297390/1‑100 (MQ=255)
                      ctcCACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATAcgc                                                                             >  2:4432/1‑100 (MQ=255)
                        ccACCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATAcgcg                                                                            <  1:1050076/100‑1 (MQ=255)
                          aCCGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCt                                                                          <  1:3202616/100‑1 (MQ=255)
                           ccGGCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtga                                                                        <  1:3625535/100‑1 (MQ=255)
                             ggCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  2:742014/100‑1 (MQ=255)
                             ggCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  1:2494353/100‑1 (MQ=255)
                             ggCTAATTAGCGATAGGGATTATCTGTACCGTGCAAATTATGGCTCGCTGGAAAAAGCTGGTATAGTGGCAAAAACACCCGTTCATAATACGCGCtgatc                                                                      <  2:2605910/100‑1 (MQ=255)
                              gCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  1:3867345/100‑1 (MQ=255)
                              gCTTATTAGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGAAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  1:1125560/100‑1 (MQ=255)
                               atctttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCtgatca                                                                     <  2:581408/96‑1 (MQ=255)
                                ttattaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCGGGTGTTGTGGCAAAAACACCCGTTCATAATACGCGCtgatcat                                                                    >  1:3367992/1‑99 (MQ=255)
                                  attaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTAGCTGGCAAAAGCTGGTAATGTGGCAAAAACACCCGTTCATAATACGCGCtgagcatgat                                                                 >  2:3497566/1‑100 (MQ=255)
                                   ttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCtgatcaagat                                                                 >  1:581412/1‑100 (MQ=255)
                                   ttaGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatcatgatc                                                                <  2:3297367/100‑1 (MQ=255)
                                    tccgatctgggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCtgatcatgatca                                                               <  2:1882887/92‑1 (MQ=255)
                                     aGCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAg                                                              >  2:4009236/1‑100 (MQ=255)
                                      gggATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTAAAAAAAAAAAAAAAACACCCGTTCATAATACGCGCtgatcatga                                                                  <  2:165299/93‑1 (MQ=255)
                                      gCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGATAAAGCTGGTATTGTGGCAAAAACACCCTTTCATAATACGCGATTATCATGATCAg                                                              >  2:328595/1‑99 (MQ=255)
                                      gCGATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAgg                                                             <  1:1245505/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  1:48959/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  1:769103/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  2:2281076/100‑1 (MQ=255)
                                        gATAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           >  2:1274266/1‑100 (MQ=255)
                                         aTAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGt                                                          >  2:3696000/1‑100 (MQ=255)
                                          tAGGGATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTc                                                         <  1:97093/100‑1 (MQ=255)
                                            gggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcac                                                       >  1:1176254/1‑100 (MQ=255)
                                            gggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcac                                                       >  2:845843/1‑100 (MQ=255)
                                             ggATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAgcgcgtatt                                                      <  1:3578503/100‑9 (MQ=255)
                                             agATTATCTGTACCGTGCCAAGTAGAGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAg                                                           <  2:1822576/94‑1 (MQ=255)
                                              gATTATCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTcaca                                                      >  1:1569279/1‑100 (MQ=255)
                                                  aTCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCAt                                                  >  2:3462229/1‑100 (MQ=255)
                                                  aTCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATg                                                 <  1:3278894/100‑1 (MQ=255)
                                                   tCTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATg                                                 <  2:3035642/100‑1 (MQ=255)
                                                    cTGTACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGAc                                               >  2:322580/1‑100 (MQ=255)
                                                       tACCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAt                                            >  2:1983501/1‑100 (MQ=255)
                                                        aCCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAtt                                           <  1:1974276/100‑1 (MQ=255)
                                                        aCCGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGGAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCAtt                                           <  2:4028753/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          <  2:2820417/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          <  1:954126/100‑1 (MQ=255)
                                                         ccGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTa                                          >  1:3237249/1‑100 (MQ=255)
                                                          cGTGCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTAc                                         >  1:2867426/1‑100 (MQ=255)
                                                             gCCAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGtc                                      >  1:1435219/1‑100 (MQ=255)
                                                              ccAATTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGtc                                      <  1:2935325/100‑1 (MQ=255)
                                                                aaTTATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTctct                                   <  1:2145007/100‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTTTGGCAAAAACTCCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTg                                  <  2:1336393/98‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGa                                <  2:1435210/100‑1 (MQ=255)
                                                                   tATGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGa                                >  1:2888872/1‑100 (MQ=255)
                                                                     tGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATg                              <  2:2958753/100‑1 (MQ=255)
                                                                     tGGCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATg                              <  2:2321964/100‑1 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATgg                             >  1:238268/1‑100 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGc                            >  2:721698/1‑100 (MQ=255)
                                                                       gCTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGc                            <  1:4019649/100‑1 (MQ=255)
                                                                        cTCGCTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCt                           >  1:16917/1‑100 (MQ=255)
                                                                            cTGGAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCg                       >  1:4047770/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgc                      >  1:80425/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgct                     >  1:3861470/1‑100 (MQ=255)
                                                                              ggAAAAAGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgct                     >  2:3370553/1‑100 (MQ=255)
                                                                                 aaaaGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctgg                   >  1:3664767/1‑100 (MQ=255)
                                                                                 aaaaGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggc                  <  1:702525/100‑1 (MQ=255)
                                                                                   aaGCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggctg                >  1:2015004/1‑100 (MQ=255)
                                                                                    aGCTGGTATTGTGGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggctg                <  1:1188148/100‑1 (MQ=255)
                                                                                     gCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCgctggcttan              <  2:3169116/100‑4 (MQ=255)
                                                                                     gCTGGTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTgaa              <  2:1723532/100‑1 (MQ=255)
                                                                                       tGGTATTGTGGCAAAAACCCCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTgaaga            <  1:80195/100‑1 (MQ=255)
                                                                                         gTATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAAc          >  2:3776313/1‑100 (MQ=255)
                                                                                          tATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACg         <  2:2072443/100‑1 (MQ=255)
                                                                                          tATTGTGGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACg         >  1:1506412/1‑100 (MQ=255)
                                                                                            ttGTGGCAAAAACACCCGTTCATAATACGGGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGat       <  2:660103/100‑1 (MQ=255)
                                                                                             tgtgGCAAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGat       >  1:1513476/1‑100 (MQ=255)
                                                                                              gtgGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  2:3432750/1‑100 (MQ=255)
                                                                                              gtgGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  1:1708107/1‑100 (MQ=255)
                                                                                              gtgGCAAAAACACCCGTTCATAATACGCGCTGATCATGATCAGGAGTCACACCATGACCATTACGTCTCTGGATGGCTATCGCTGGCTGAAGAACGatat     >  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: