Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,679 G→T 16.4% intergenic (+68/‑779) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,6790GT16.4% 434.7 / 68.4 165intergenic (+68/‑779)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (62/75);  new base T (10/17);  total (72/93)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 166 total aligned reads displayed.

TCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATT  >  NC_000913/770582‑770777
                                                                                                 |                                                                                                  
  aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                <  1:2351096/100‑1 (MQ=255)
   cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                               >  1:637154/1‑100 (MQ=255)
   cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                               >  1:2144138/1‑100 (MQ=255)
   cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                               >  2:195143/1‑100 (MQ=255)
     gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                             >  1:2987209/1‑100 (MQ=255)
      aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                            <  2:414914/100‑1 (MQ=255)
      aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                            <  1:3138462/100‑1 (MQ=255)
        gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                            <  2:265619/98‑1 (MQ=255)
        gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATg                                                                                          <  2:3152119/100‑1 (MQ=255)
        gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATg                                                                                          <  1:949728/100‑1 (MQ=255)
         ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                         <  1:680229/100‑1 (MQ=255)
          tCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGaa                                                                                        <  2:1860914/100‑1 (MQ=255)
           cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                       <  1:2949424/100‑1 (MQ=255)
           cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                       <  1:687238/100‑1 (MQ=255)
           cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                       <  2:626097/100‑1 (MQ=255)
           cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAAc                                                                                       >  1:982021/1‑100 (MQ=255)
           cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAAc                                                                                       >  1:3562913/1‑100 (MQ=255)
            aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                      <  2:2441773/100‑2 (MQ=255)
            aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                      <  1:3595728/100‑2 (MQ=255)
            aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                      <  2:3053800/100‑2 (MQ=255)
            aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCACTTTGTTCATGCCGTATCCGGCATGAAca                                                                                      <  2:3206192/100‑2 (MQ=255)
             gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                     <  2:1547930/100‑1 (MQ=255)
             gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                     >  2:3348279/1‑100 (MQ=255)
              ttACCGGCTTGCCTGAATAGCAATCAAACCGAATCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaa                                                                                     <  2:2872105/99‑3 (MQ=38)
              ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                    >  2:224107/1‑100 (MQ=255)
              ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                    <  2:3422510/100‑1 (MQ=255)
              ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                    <  2:3270598/100‑1 (MQ=255)
               tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaaag                                                                                   <  1:2313198/100‑5 (MQ=39)
               tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                   <  1:2068396/100‑1 (MQ=255)
               tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                   <  2:144528/100‑1 (MQ=255)
                 ccGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTa                                                                                 >  1:1050216/1‑100 (MQ=255)
                  cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTa                                                                                 <  2:2008533/99‑1 (MQ=255)
                  cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                <  1:1279723/100‑1 (MQ=255)
                  cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                <  2:925443/100‑1 (MQ=255)
                   ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAtt                                                                               <  2:1399405/100‑1 (MQ=255)
                   ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                >  1:1138875/1‑99 (MQ=255)
                      ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                            <  1:134294/100‑1 (MQ=255)
                      ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                            <  1:262303/100‑1 (MQ=255)
                       tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                           >  1:586191/1‑100 (MQ=255)
                         ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                         <  1:2613961/100‑1 (MQ=255)
                           tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                       <  2:3294506/100‑1 (MQ=255)
                           tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                       >  2:563652/1‑100 (MQ=255)
                              aTAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCAGTCTACaaa                                                                    >  1:2873566/1‑100 (MQ=255)
                                aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                  <  1:1788652/100‑1 (MQ=255)
                                aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                  <  2:3493625/100‑1 (MQ=255)
                                aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                  <  2:2046323/100‑1 (MQ=255)
                                 gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                 <  1:1578544/100‑1 (MQ=255)
                                 gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                 >  2:2876299/1‑100 (MQ=255)
                                  cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                <  2:1759328/100‑1 (MQ=255)
                                   aaTCAAACCGCAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                               <  2:2383763/100‑1 (MQ=255)
                                   aaTCAAACCGAAGCCACATATGCGCGGCCAGAGTGTTGACAAAGGGCGCTTTATTCAGGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                               <  2:2329783/100‑1 (MQ=255)
                                     tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                             <  1:2184802/100‑1 (MQ=255)
                                      cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCTTGAACGCTTTATTCGGTCTACACAAGc                                                                 >  2:3588975/1‑95 (MQ=255)
                                      cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCa                                                            <  1:3560033/100‑1 (MQ=255)
                                       aaaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCaa                                                           >  1:1479181/1‑100 (MQ=255)
                                       aaaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCaa                                                           >  1:1687891/1‑100 (MQ=255)
                                        aaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCaaa                                                          >  2:1163532/1‑100 (MQ=255)
                                        aaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAAGGCACTTTGTTCCTGACGGATACAGCATGCAAGCTTAATTCATTCTACAAAAGCAGGCaaa                                                          >  2:696821/1‑100 (MQ=255)
                                          ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                        <  1:14003/100‑1 (MQ=255)
                                          ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                        >  1:348520/1‑100 (MQ=255)
                                           cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                       <  2:583744/100‑1 (MQ=255)
                                              aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCTGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                    <  2:2289764/100‑1 (MQ=255)
                                              aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCTGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCATGCAAATTCAat                                                    <  2:742690/100‑1 (MQ=255)
                                               gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                                   <  1:454880/100‑1 (MQ=255)
                                                  acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCCGGCAAATTAAatatat                                                >  2:3477080/1‑100 (MQ=255)
                                                   caTATGCCCGGCCAGACTGTTGACAAAGGGCGCTTTGTTCATGCAGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTAAAGATAtt                                               >  2:3522753/1‑100 (MQ=25)
                                                     tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                             <  2:2092855/100‑1 (MQ=255)
                                                     tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                             <  2:2655566/100‑1 (MQ=255)
                                                      atGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                            >  2:41086/1‑100 (MQ=255)
                                                       tGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCag                                           <  1:2434191/100‑1 (MQ=255)
                                                        gcgcGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCaga                                          >  1:546719/1‑100 (MQ=255)
                                                         cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                         <  2:1423209/100‑1 (MQ=255)
                                                         cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                         >  1:1333488/1‑100 (MQ=255)
                                                             gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                     <  1:3335712/100‑1 (MQ=255)
                                                              ccAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatgat                                    <  1:570843/100‑1 (MQ=255)
                                                              ccAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatgat                                    <  1:1712960/100‑1 (MQ=255)
                                                              ccAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatgat                                    >  1:1266179/1‑100 (MQ=255)
                                                               cAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGatgatg                                   <  2:2128313/100‑1 (MQ=255)
                                                                  aTTGTTGACAAAGGGCGTTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAg                                <  2:428098/100‑1 (MQ=255)
                                                                    tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGac                              >  1:409651/1‑100 (MQ=255)
                                                                        gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                          >  1:2805204/1‑100 (MQ=255)
                                                                        gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                          <  1:2750163/100‑1 (MQ=255)
                                                                                  gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAg                >  2:975533/1‑100 (MQ=255)
                                                                                   cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc               >  1:3218718/1‑100 (MQ=255)
                                                                                    tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg              <  2:2080030/100‑1 (MQ=255)
                                                                                    agcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg              >  2:819861/4‑100 (MQ=255)
                                                                                      tGTTCATGCCGGATGCGGCATGAACGCTGTATGCGGTCTACAAAGGCAGTGAAATTCAATATGTTGCAGAGATGATGTAGAGACTGGCAAGCGTAGCGCa            <  2:1255732/100‑1 (MQ=21)
                                                                                      tGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcgc             >  2:926151/1‑99 (MQ=255)
                                                                                      cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa            >  2:3142161/2‑100 (MQ=255)
                                                                                      cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa            >  1:1444640/2‑100 (MQ=255)
                                                                                      cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa            >  1:2545873/2‑100 (MQ=255)
                                                                                      cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa            >  2:2363689/2‑100 (MQ=255)
                                                                                       gTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt           >  1:407882/1‑100 (MQ=255)
                                                                                        ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATc          <  2:2787334/100‑1 (MQ=255)
                                                                                        ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATc          >  2:3398590/1‑100 (MQ=255)
                                                                                         tCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCa         >  2:1643635/1‑100 (MQ=255)
                                                                                            tGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt      >  1:1347725/1‑100 (MQ=255)
                                                                                            tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAACTTCAATCTATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt      >  1:1890566/1‑100 (MQ=255)
                                                                                             gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa     >  1:184746/1‑100 (MQ=255)
                                                                                                ggATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAtt  <  1:282627/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATT  >  NC_000913/770582‑770777

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