Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,016:1 +C 16.5% intergenic (+52/+32) yfbP → / ← nuoN TPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,390,0161.C16.5% 546.3 / 113.8 158intergenic (+52/+32)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base . (61/71);  new base C (13/13);  total (74/84)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 165 total aligned reads displayed.

CAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTA  >  NC_000913/2389917‑2390115
                                                                                                    |                                                                                                     
cAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATg                                                                                                        >  2:1319947/1‑100 (MQ=255)
cAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATg                                                                                                        <  1:1359526/100‑1 (MQ=255)
cAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATg                                                                                                        <  1:1359555/100‑1 (MQ=255)
    tATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcgg                                                                                                   <  1:2829225/100‑1 (MQ=255)
     aTGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggc                                                                                                  <  2:2895970/100‑1 (MQ=255)
      tGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggcg                                                                                                 <  1:338367/100‑1 (MQ=255)
       ggttAGAGTAGGTTGGCGATAGTCGTACACAGTTGAGTTAACGCTATGCGTCTGTCTACTATCCCTTATTTGACACGAGAACCGGCGAGCATGTcggcga                                                                                                <  2:2302861/98‑2 (MQ=255)
       gTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggcgg                                                                                                >  1:1072952/1‑100 (MQ=255)
        tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCggcggc                                                                                              >  2:1115228/1‑100 (MQ=255)
        tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCggcggc                                                                                              <  2:2824414/100‑1 (MQ=255)
         ttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcggtt                                                                                              <  1:1534585/100‑3 (MQ=255)
          tAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGttt                                                                                           <  1:2237582/100‑1 (MQ=255)
           agagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCttt                                                                                          >  1:3151070/1‑100 (MQ=255)
           agagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGtttt                                                                                          <  1:2148358/100‑1 (MQ=255)
            gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTc                                                                                         <  1:3169085/100‑1 (MQ=255)
            gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGGCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCggcggc                                                                                              <  2:198983/96‑1 (MQ=255)
             agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCt                                                                                        <  2:401269/100‑1 (MQ=255)
             agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGtttttc                                                                                        <  1:2727489/100‑3 (MQ=38)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGtttttct                                                                                       <  1:2543545/100‑4 (MQ=255)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTg                                                                                       <  2:3068363/100‑1 (MQ=255)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGtttttct                                                                                       <  1:1253871/100‑4 (MQ=38)
               tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCt                                                                                        <  2:3555650/98‑1 (MQ=255)
                aTTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAGCCGCCGACATGCTCGGCGGtttttctgaa                                                                                    <  2:1359149/100‑7 (MQ=21)
                aTTGTGGCGATAGTCGTACACAGTTCCCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGAc                                                                                     <  2:3230978/100‑1 (MQ=255)
                 ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACt                                                                                    >  2:3578190/1‑100 (MQ=255)
                  tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACtt                                                                                   <  1:1691541/100‑1 (MQ=255)
                  tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAGCCGCCGACATGCTCGGCGGtttttctgaaaa                                                                                  <  1:1178273/100‑9 (MQ=21)
                  tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAGCCGCCGACATGCTCGGCGGtttttctgaaaa                                                                                  <  1:2071804/100‑9 (MQ=18)
                   ttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTa                                                                                  <  2:3393522/100‑1 (MQ=255)
                      gCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaa                                                                               <  2:2704286/100‑1 (MQ=255)
                         tctgTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGt                                                                            <  1:1731089/97‑1 (MQ=255)
                          tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                           <  2:1273005/100‑1 (MQ=255)
                          tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                           <  1:2518403/100‑1 (MQ=255)
                          tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                           <  2:2669169/100‑1 (MQ=255)
                          tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTATTCAGCAAACCCGCCAAGCATGTCGGCGGCTTTCTAACTTACAAAGTa                                                                           >  2:700722/1‑100 (MQ=255)
                            gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTaga                                                                         >  2:1731081/1‑98 (MQ=255)
                            gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                         >  1:122922/1‑100 (MQ=255)
                            gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                         <  2:2180652/100‑1 (MQ=255)
                            gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                         >  1:119219/1‑100 (MQ=255)
                              cGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAg                                                                       >  1:1933059/1‑100 (MQ=255)
                               gTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGa                                                                      >  2:686451/1‑100 (MQ=255)
                                tACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAt                                                                     >  1:1101994/1‑100 (MQ=255)
                                 acacAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAtt                                                                    <  1:791207/100‑1 (MQ=255)
                                 acacAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAtt                                                                    >  1:660280/1‑100 (MQ=255)
                                  cacaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTa                                                                   <  1:2083576/100‑1 (MQ=255)
                                   acacTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTAc                                                                  <  2:905452/100‑1 (MQ=255)
                                    caGTTCACTTAACGCTCTGCGTCTAGCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACa                                                                 <  2:1623518/100‑1 (MQ=255)
                                      gTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATc                                                               >  2:2472546/1‑100 (MQ=255)
                                            ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                         >  1:3606805/1‑100 (MQ=255)
                                            ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                         >  2:3276239/1‑100 (MQ=255)
                                             tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                        <  1:1587197/100‑1 (MQ=255)
                                             tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCACCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                         >  2:2985154/1‑99 (MQ=255)
                                             tAACGCTCTGCGTCTATCTACTACCCCTTCTTTTCCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                        <  2:599730/100‑1 (MQ=255)
                                                  ctctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGcc                                                   >  1:959586/1‑100 (MQ=255)
                                                    ctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaa                                                 >  2:2904520/1‑100 (MQ=255)
                                                      gCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAAc                                               >  2:2834748/1‑100 (MQ=255)
                                                        gTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGc                                             <  1:1931021/100‑1 (MQ=255)
                                                         tctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGca                                            >  1:3291100/1‑100 (MQ=255)
                                                          ctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGcac                                           <  2:1798013/100‑1 (MQ=255)
                                                              ctactaTCCCTTATTTTTCAGAAAAACCGCCGAGCATGTCGACGGCTTTTTGACTTAAAAAATAACAGATTACATCAGCGGCATTGCCAAACGCACaa                                         >  2:2343950/1‑98 (MQ=255)
                                                               tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                      >  1:583112/1‑100 (MQ=255)
                                                                actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                     >  1:3223562/1‑100 (MQ=255)
                                                                actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                     >  2:352087/1‑100 (MQ=255)
                                                                      ccTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAg                               >  2:257614/1‑100 (MQ=255)
                                                                       cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGTCATTGCTAAACGCACAATGCTAATCAg                               >  2:430868/1‑99 (MQ=255)
                                                                       cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGc                              >  1:3628456/1‑100 (MQ=255)
                                                                           tttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGACATTGCCAAACGCACAATGCTAATCAGCGGtt                          >  2:512842/1‑100 (MQ=255)
                                                                              ttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgtg                       <  2:1700498/100‑1 (MQ=255)
                                                                              agtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCCGCGGTtgt                        >  1:1510463/11‑100 (MQ=18)
                                                                              agtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  2:116329/11‑100 (MQ=21)
                                                                              agtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  1:827037/11‑100 (MQ=21)
                                                                              agtcagaaagCCGCCGACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  2:2232320/7‑100 (MQ=21)
                                                                               tCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      <  2:2802691/100‑1 (MQ=255)
                                                                               tCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGAATTACAAAGTAACAGATTACAACAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      >  1:304346/1‑100 (MQ=255)
                                                                                tcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      >  1:1281711/9‑100 (MQ=21)
                                                                                gaGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATAGCCAAACGCACAATGCTAATCAGCGGTCGTGGc                     >  2:648138/2‑100 (MQ=255)
                                                                                cggAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  1:2912965/3‑100 (MQ=255)
                                                                                cAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  1:3121419/1‑100 (MQ=255)
                                                                                cAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  1:2200413/1‑100 (MQ=255)
                                                                                cAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGCTGTGGc                     >  2:1527771/1‑100 (MQ=255)
                                                                                 cagaaagccgccgACATGCTCGGCGGTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGcc                    >  2:2074672/8‑100 (MQ=37)
                                                                                 cagaaagCCGCCGACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  1:2025525/4‑100 (MQ=25)
                                                                                 aGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     <  2:1911299/99‑1 (MQ=255)
                                                                                  gAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCa                   >  2:1345329/1‑100 (MQ=255)
                                                                                    aaaaCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                 <  1:1817061/100‑1 (MQ=255)
                                                                                      agccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATac                >  1:806171/3‑100 (MQ=37)
                                                                                      aaCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATaca               <  1:1654873/100‑1 (MQ=255)
                                                                                      aaCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATaca               >  1:2329023/1‑100 (MQ=255)
                                                                                       aCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac              >  1:1460624/1‑100 (MQ=255)
                                                                                       aCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac              >  2:2326504/1‑100 (MQ=255)
                                                                                       aCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATACTAATCAGCGGTTGTGGCCATacac              <  1:1326309/100‑1 (MQ=255)
                                                                                        ccgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAcc             <  1:304396/100‑1 (MQ=255)
                                                                                            cgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc         <  1:3566177/100‑1 (MQ=255)
                                                                                            cgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCCGCATGGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc         <  2:619369/100‑1 (MQ=255)
                                                                                              aGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAc       <  1:2274246/100‑1 (MQ=255)
                                                                                               gCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAcc      <  1:3025179/100‑1 (MQ=255)
                                                                                                 aTGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAg    <  1:921023/100‑1 (MQ=255)
                                                                                                  tGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGt   >  2:1327350/1‑100 (MQ=255)
                                                                                                   gTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTa  <  2:1336993/100‑1 (MQ=255)
                                                                                                   gTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTa  >  2:297113/1‑100 (MQ=255)
                                                                                                    |                                                                                                     
CAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTA  >  NC_000913/2389917‑2390115

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: