Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 18.5% intergenic (‑150/+17) ygfI ← / ← yggE putative DNA‑binding transcriptional regulator/oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3230AT18.5% 481.7 / 96.5 205intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (68/99);  new base T (17/21);  total (85/120)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

Only 100 of 208 total aligned reads displayed.

TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAA  >  NC_000913/3067224‑3067422
                                                                                                   |                                                                                                   
tgatgaACCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                     <  1:1522905/100‑1 (MQ=255)
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                     <  1:1370746/100‑1 (MQ=255)
  atgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGCGTAAAGGCATTAGCCCGACGAATACg                                                                                                   <  1:2549585/100‑1 (MQ=255)
   tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACg                                                                                                   >  2:1944273/1‑99 (MQ=255)
   tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCAGTAGCCCGACGTATTCg                                                                                                   >  2:1436993/1‑99 (MQ=255)
   tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGACATTAGCCCGACGAATACGt                                                                                                  >  1:1479153/1‑100 (MQ=255)
      aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                               >  1:3226547/1‑100 (MQ=255)
       aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                              <  2:3290534/100‑1 (MQ=255)
       aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCggg                                                                                              <  1:697085/100‑1 (MQ=255)
        ccGCGAGCAGAGATAGTTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                             <  1:3498345/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  2:127780/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  1:2713219/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  1:945800/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                             <  1:1769619/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                             <  1:544097/100‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                            <  2:1313113/100‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                            <  1:2426248/100‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                            <  1:1078906/100‑1 (MQ=255)
          gcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTa                                                                                           >  2:3936912/1‑100 (MQ=255)
           cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                          >  1:4006032/1‑100 (MQ=255)
            gAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACa                                                                                         <  1:2198414/100‑1 (MQ=255)
              gCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACaaa                                                                                       <  1:2338277/100‑1 (MQ=255)
               cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                      <  1:1982876/100‑1 (MQ=255)
                  aGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtatt                                                                                   <  2:787063/100‑1 (MQ=255)
                     tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                                <  2:2992260/100‑1 (MQ=255)
                      aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                               >  2:86584/1‑100 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCACTATCGGCTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                            >  2:3614029/1‑100 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                            <  2:1113772/100‑1 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                            <  1:3936526/100‑1 (MQ=255)
                              aCGTCTTGGTTTTACAGAAAATTACAAAAAAAGGAGGCAATATAGGATAAAGGAATAAGCACGACAAATACGTCATGATACAAATAGTATTGTGCTGAAg                                                                       >  2:1752850/1‑100 (MQ=255)
                                gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTACAgg                                                                      >  1:3799688/1‑99 (MQ=255)
                                 tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                                                    >  2:1526076/1‑100 (MQ=255)
                                   aTGGTTTTACAGAAAATTACAGAAAAAGGAGGAAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                                  <  2:1075913/100‑1 (MQ=255)
                                     ggTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                                                <  2:551033/100‑1 (MQ=255)
                                       ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                                                              >  1:1001403/1‑100 (MQ=255)
                                         ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                            <  2:4104109/100‑1 (MQ=255)
                                         ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                            <  2:1297885/100‑1 (MQ=255)
                                           acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                          <  1:1482251/100‑1 (MQ=255)
                                           acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                          >  2:2638527/1‑100 (MQ=255)
                                           acagaaaattacagaaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                         >  1:1743912/1‑100 (MQ=255)
                                            cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                         <  1:2643905/100‑1 (MQ=255)
                                            cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                         >  2:2731554/1‑100 (MQ=255)
                                                aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                     <  1:1627311/100‑1 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgttn                                                    <  1:3462397/100‑2 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                    <  1:3781364/100‑1 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                    >  1:2477150/1‑100 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                    >  1:408262/1‑100 (MQ=255)
                                                  attacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGa                                                   <  1:1080161/100‑1 (MQ=255)
                                                  agtacagaaaaAGGAGGCAAGATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATGAGTGTGCTGCAGGTGTTTTAGCGGGTTGTTGa                                                   <  1:1950949/98‑1 (MQ=255)
                                                    tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                                                 <  1:3659288/100‑1 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  2:422348/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  2:3035358/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  1:2018166/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               >  1:889005/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               <  1:3745230/100‑1 (MQ=255)
                                                       agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcccc                                              >  2:1258664/1‑98 (MQ=255)
                                                       agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                               <  2:3684524/99‑1 (MQ=255)
                                                        gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                             <  2:4006005/100‑1 (MQ=255)
                                                        gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                             >  2:3821950/1‑100 (MQ=255)
                                                        gaaaaAGGAGGAAATATCGGGTAAAGGCATTATCCAGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                             <  2:769122/100‑1 (MQ=255)
                                                         aaaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAg                                            <  1:2937609/100‑1 (MQ=255)
                                                           aaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATGATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGt                                          <  2:2181890/100‑1 (MQ=255)
                                                             aggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTc                                        <  1:668994/100‑1 (MQ=255)
                                                              ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                       >  2:2334328/1‑100 (MQ=255)
                                                              ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                       <  1:1664977/100‑1 (MQ=255)
                                                                 ggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                                    >  1:3363512/1‑100 (MQ=255)
                                                                 ggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                                    <  1:995255/100‑1 (MQ=255)
                                                                  gCAATATCGGGTAAAGCCATTAGCCCGCCGAATACGTCGGGCTACAAATATCATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                   <  1:189516/100‑1 (MQ=255)
                                                                   cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                  >  2:3170590/1‑100 (MQ=255)
                                                                    aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTgg                                 <  1:3237667/100‑1 (MQ=255)
                                                                    aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTgg                                 <  1:2882110/100‑1 (MQ=255)
                                                                      tatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                               <  2:2900668/100‑1 (MQ=255)
                                                                       atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                              >  1:3846263/1‑100 (MQ=255)
                                                                       atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                              <  1:1168198/100‑1 (MQ=255)
                                                                         cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                            >  1:3912982/1‑100 (MQ=255)
                                                                         cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                            >  2:2024506/1‑100 (MQ=255)
                                                                         cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                            >  1:753081/1‑100 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                           <  1:2138709/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                           <  1:4048121/100‑1 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                          >  1:775388/1‑100 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                          <  1:4108667/100‑1 (MQ=255)
                                                                                aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCg                     <  2:3088297/100‑1 (MQ=255)
                                                                                   cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                  <  1:4066446/100‑1 (MQ=255)
                                                                                   cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                  <  1:3341299/100‑1 (MQ=255)
                                                                                     ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg                <  1:3685858/100‑1 (MQ=255)
                                                                                       aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat              >  1:2093528/1‑100 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGata             >  1:3102675/1‑99 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTATTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  2:1797545/1‑100 (MQ=255)
                                                                                        gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTCTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  2:3769989/1‑100 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  1:487551/1‑100 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             <  1:2585546/100‑1 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             <  2:4083805/100‑1 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc             >  2:1478171/1‑100 (MQ=255)
                                                                                         cccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca            >  1:479152/1‑100 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTACGGTGCTTCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTGCTAACTAGAAGACCACATCGACCTGatca            <  2:491451/100‑1 (MQ=255)
                                                                                            gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGAtcatca         <  1:2811203/100‑1 (MQ=255)
                                                                                              cgaATACGTCGGGCTACAAGTATTATTGTGCTGCAGGTGTTTTCGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa       <  2:2715595/100‑1 (MQ=255)
                                                                                                aaTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACt     >  1:1679514/1‑100 (MQ=255)
                                                                                                  tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa   >  1:303573/1‑100 (MQ=255)
                                                                                                  tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa   <  1:388483/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAA  >  NC_000913/3067224‑3067422

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