Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,712,518 (C)6→5 100% intergenic (+360/‑251) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,5130C.99.0% 474.7 / ‑2.1 103intergenic (+355/‑256)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (58/44);  total (58/45)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 116 total aligned reads displayed.

TATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCA  >  NC_000913/1712414‑1712608
                                                                                                   |                                                                                                
tatTGAGCTACAAAACATTACACTGGCTATTTTTCAGAAACTGTACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAc                                                                                                  <  2:2395192/100‑1 (MQ=255)
tatTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAc                                                                                                  >  2:1148837/1‑100 (MQ=255)
    gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                              >  1:3987537/1‑100 (MQ=255)
    gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                              <  2:348347/100‑1 (MQ=255)
    gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                              <  2:3296618/100‑1 (MQ=255)
    gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcccac                                                                                              >  2:3869918/1‑98 (MQ=255)
     aCCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccct                                                                                             >  1:905572/1‑99 (MQ=255)
     aCCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcccccc                                                                                             <  1:3696070/100‑1 (MQ=255)
      ccTACAAAACATTACACTGGCTATTTTTTTTTTTTTTTACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                              >  2:697512/1‑98 (MQ=255)
      ccTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccctc                                                                                            <  1:1588989/100‑3 (MQ=255)
       cTACAAAACATTACACTGGCTATTTTTTTTTTTTTTTTCTATCTCAATAGTAAAAACGCTAAATATCAGAACATATCGCCGTTAAGCGATCAccccctc                                                                                            >  2:3654428/1‑97 (MQ=255)
       cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccctct                                                                                           <  1:2939720/100‑4 (MQ=255)
         aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATATCAGAACATATCGCCATTAAGCCATTActccctctc                                                                                          >  2:2040698/1‑95 (MQ=255)
         aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTca                                                                                        >  2:3890549/1‑100 (MQ=255)
         aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTca                                                                                        >  1:3125122/1‑100 (MQ=255)
         aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTca                                                                                        >  1:1440401/1‑100 (MQ=255)
           aaaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagt                                                                                      <  2:3960395/100‑1 (MQ=255)
           aaaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagt                                                                                      <  2:3369092/100‑1 (MQ=255)
            aaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtg                                                                                     <  2:3791678/100‑1 (MQ=255)
              aCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgca                                                                                   <  2:358852/100‑1 (MQ=255)
              aCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgc                                                                                    <  2:513788/99‑1 (MQ=255)
               cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcag                                                                                  <  2:3259317/100‑1 (MQ=255)
               cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcag                                                                                  <  1:1689994/100‑1 (MQ=255)
                aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcagt                                                                                 >  1:881600/1‑100 (MQ=255)
                   acacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑aa                                                                             >  1:3009144/1‑100 (MQ=255)
                    cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑aaa                                                                            <  1:2459279/100‑1 (MQ=255)
                          cTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATc                                                                       <  2:3879040/99‑1 (MQ=255)
                          cTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCt                                                                      >  1:3241995/1‑100 (MQ=255)
                           tATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTg                                                                     <  2:4049751/100‑1 (MQ=255)
                             tttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGc                                                                    >  2:2570863/1‑99 (MQ=255)
                                ttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGtt                                                                >  1:1433449/1‑100 (MQ=255)
                                  cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTAc                                                              >  2:3639477/1‑100 (MQ=255)
                                  cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTAc                                                              >  2:271895/1‑100 (MQ=255)
                                  cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTAc                                                              >  1:1543454/1‑100 (MQ=255)
                                   aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACg                                                             >  2:2802865/1‑100 (MQ=255)
                                    gAAACTGGACTATCTCACTAGTAAAAACGCGAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CACACACTCACTGCAGTG‑AAAAAATCTGCCGTTACg                                                             >  2:3779605/1‑99 (MQ=255)
                                      aaCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCCTCA‑CCCCCTCTCAATGCATTG‑AAAAAATCTGACGTTACAtt                                                           >  2:1139678/1‑99 (MQ=255)
                                        cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAAATCTGCCGTTACGtttt                                                         >  2:1448/1‑100 (MQ=255)
                                          ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGtttttt                                                       >  1:53797/1‑99 (MQ=255)
                                          ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTg                                                      <  1:1109514/100‑1 (MQ=255)
                                          ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTg                                                      <  1:52978/100‑1 (MQ=255)
                                          ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTg                                                      <  1:885747/100‑1 (MQ=255)
                                           gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTg                                                      <  2:3987511/99‑1 (MQ=255)
                                            aCTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CACACTCTCAGTGCAGTG‑AAAAAATCTACCGTTAAGTTTTTTg                                                      <  2:623259/98‑1 (MQ=255)
                                             cTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGaaa                                                   <  1:4098489/100‑1 (MQ=255)
                                               aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGaaaa                                                  >  2:3423709/1‑99 (MQ=255)
                                                 ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAATAAtt                                               >  2:2615561/1‑100 (MQ=255)
                                                 ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAAtt                                               >  1:3298432/1‑100 (MQ=255)
                                                  tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGAAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                              >  2:2890426/1‑100 (MQ=255)
                                                   cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                             >  1:3422072/1‑100 (MQ=255)
                                                    aatAGTAAAAACGCTAAATAGCAGAACATAGCGCCGTTAAGCCATCA‑CCCCCTCTAAAAACAGAA‑AAAAAATCTGCCGTTACGTTTTTTGAAAAAt                                                <  2:838008/94‑1 (MQ=255)
                                                      tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAAc                                           >  2:1693214/1‑99 (MQ=255)
                                                        gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                        >  2:1120708/1‑100 (MQ=255)
                                                        gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                        <  2:790280/100‑1 (MQ=255)
                                                          aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                      <  1:875783/100‑1 (MQ=255)
                                                          aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                      >  1:1803441/1‑100 (MQ=255)
                                                          aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                      <  2:1619412/100‑1 (MQ=255)
                                                          aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                      >  1:3730906/1‑100 (MQ=255)
                                                          aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGACATCA‑CCCCCTCTCAGTGCAGTG‑CAAAAATCTGCCGTTACGTTTTTTGCAAAATTTAACGCTgg                                      <  2:1219950/100‑1 (MQ=255)
                                                           aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                     >  1:4106947/1‑100 (MQ=255)
                                                           aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                     <  1:3402971/100‑1 (MQ=255)
                                                             aaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                    <  1:1148844/100‑1 (MQ=255)
                                                               cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc                                 <  1:2322071/100‑1 (MQ=255)
                                                               cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc                                 >  1:3738092/1‑100 (MQ=255)
                                                                gCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                <  2:3749203/100‑1 (MQ=255)
                                                                 cTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAACTTTAACGCTGGATAACAt                               >  2:1592813/1‑100 (MQ=255)
                                                                  tAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAtt                              >  1:3998499/1‑100 (MQ=255)
                                                                   aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                             >  1:1050767/1‑100 (MQ=255)
                                                                   aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                             >  1:2104364/1‑100 (MQ=255)
                                                                   aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                             >  1:1947982/1‑100 (MQ=255)
                                                                   aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                             >  2:3377826/1‑100 (MQ=255)
                                                                      tAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAAGCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTccc                          <  2:2295477/100‑1 (MQ=255)
                                                                       aGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                         >  1:595684/1‑100 (MQ=255)
                                                                        gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATATCATTTCCCgg                        >  2:2705784/1‑100 (MQ=255)
                                                                        gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                        <  1:3377846/100‑1 (MQ=255)
                                                                        gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                        >  1:2885229/1‑100 (MQ=255)
                                                                        gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                        >  1:2356382/1‑100 (MQ=255)
                                                                        gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                        >  1:3819180/1‑100 (MQ=255)
                                                                         cAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                       <  1:1184204/100‑1 (MQ=255)
                                                                           gAACATATCGCCGTTAAGCCCTCA‑CCCCCTCTCCGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                        <  2:2404173/97‑1 (MQ=255)
                                                                            aaCATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                    >  2:1774220/1‑100 (MQ=255)
                                                                              gatctCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGt                  <  2:3202248/96‑1 (MQ=255)
                                                                              cATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGt                  >  1:2968401/1‑100 (MQ=255)
                                                                               atatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtt                 <  1:1345458/100‑1 (MQ=255)
                                                                                tatTGCCGTTAAGCCATCA‑ACCCCTGTCAGTGAAATG‑GAAAAATCTGCCGTTACGTTGTTTGAAAAATTGAACGCTGGATAACATTTCCCGGAATGGTTg                <  2:77844/100‑1 (MQ=255)
                                                                                tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                <  2:1114760/100‑1 (MQ=255)
                                                                                tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                <  2:1114831/100‑1 (MQ=255)
                                                                                tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                <  1:577099/100‑1 (MQ=255)
                                                                                tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                >  1:3123619/1‑100 (MQ=255)
                                                                                 atCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTACCGCTGGATAACATTTCCCGGAATGGTTg                >  2:1747230/1‑99 (MQ=255)
                                                                                 atCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGa               >  1:1242393/1‑100 (MQ=255)
                                                                                   ngCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt             >  1:1709919/2‑100 (MQ=255)
                                                                                   caccattAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt             <  2:1235797/95‑1 (MQ=255)
                                                                                   cGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtagatc             >  1:3202267/1‑95 (MQ=255)
                                                                                    gCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt            >  2:2952267/1‑100 (MQ=255)
                                                                                      cGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcc          <  1:531445/100‑1 (MQ=255)
                                                                                           aGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTTTTTTTTTTTACACTGGAAAACATTTCCCGGAATGGTTTAATTcc          >  2:2342759/1‑95 (MQ=255)
                                                                                           aGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc     <  1:3260850/100‑1 (MQ=255)
                                                                                              cATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCa  >  1:554152/1‑100 (MQ=255)
                                                                                              cATCA‑CCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCa  >  1:1067578/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTG‑AAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCA  >  NC_000913/1712414‑1712608

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: ‑