Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,682 A→C 16.5% intergenic (+71/‑776) 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,6820AC16.5% 402.1 / 76.0 165intergenic (+71/‑776)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (62/75);  new base C (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 = 7.07e-01

Only 100 of 168 total aligned reads displayed.

CACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGC  >  NC_000913/770583‑770781
                                                                                                   |                                                                                                   
cACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATa                                                                                                     <  1:2916974/100‑1 (MQ=255)
 aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                    <  1:2351096/100‑1 (MQ=255)
  cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                                   >  1:637154/1‑100 (MQ=255)
  cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                                   >  1:2144138/1‑100 (MQ=255)
    gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                                 >  1:2987209/1‑100 (MQ=255)
     aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                <  1:3138462/100‑1 (MQ=255)
     aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                <  2:414914/100‑1 (MQ=255)
       gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                <  2:265619/98‑1 (MQ=255)
       gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATg                                                                                              <  1:949728/100‑1 (MQ=255)
       gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATg                                                                                              <  2:3152119/100‑1 (MQ=255)
        ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                             <  2:338669/100‑1 (MQ=255)
          cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  1:687238/100‑1 (MQ=255)
          cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  2:626097/100‑1 (MQ=255)
          cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  1:2836688/100‑1 (MQ=255)
          cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAAc                                                                                           >  1:3562913/1‑100 (MQ=255)
           aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                          <  1:3595728/100‑2 (MQ=255)
           aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                          <  2:3053800/100‑2 (MQ=255)
           aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAca                                                                                          <  2:2441773/100‑2 (MQ=255)
           aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCACTTTGTTCATGCCGTATCCGGCATGAAca                                                                                          <  2:3206192/100‑2 (MQ=255)
            gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                         >  2:3348279/1‑100 (MQ=255)
             ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                        >  2:224107/1‑100 (MQ=255)
             ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                        <  2:3422510/100‑1 (MQ=255)
              tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaaag                                                                                       <  1:2313198/100‑5 (MQ=39)
              tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                       <  2:144528/100‑1 (MQ=255)
              tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                       <  1:2068396/100‑1 (MQ=255)
                ccGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTa                                                                                     >  1:1050216/1‑100 (MQ=255)
                 cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTa                                                                                     <  2:2008533/99‑1 (MQ=255)
                 cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                    <  1:3534812/100‑1 (MQ=255)
                  ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                    >  1:1138875/1‑99 (MQ=255)
                     ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                <  1:134294/100‑1 (MQ=255)
                      tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                               >  1:586191/1‑100 (MQ=255)
                       gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                              <  2:2796441/100‑1 (MQ=255)
                        ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                             >  2:1457755/1‑100 (MQ=255)
                        ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                             <  1:2613961/100‑1 (MQ=255)
                          tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                           <  2:877999/100‑1 (MQ=255)
                          tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                           <  2:3294506/100‑1 (MQ=255)
                           gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCTCTTTGTTCACGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                           <  2:3615628/99‑1 (MQ=255)
                           gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                          >  1:3181125/1‑100 (MQ=255)
                             aTAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaa                                                                        >  2:1419281/1‑100 (MQ=255)
                              tAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaaa                                                                       <  1:132212/100‑1 (MQ=255)
                              tAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaaa                                                                       >  2:1852451/1‑100 (MQ=255)
                               aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                      <  1:1788652/100‑1 (MQ=255)
                                gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTTTTCATTCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                     >  2:2009052/1‑100 (MQ=255)
                                gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                     <  1:1578544/100‑1 (MQ=255)
                                gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                     >  2:2876299/1‑100 (MQ=255)
                                 cAATCAAACCGAAGCCATATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                    <  2:1586673/100‑1 (MQ=255)
                                 cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                    <  2:1759328/100‑1 (MQ=255)
                                  aaTCAAACCGCAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                   <  2:2383763/100‑1 (MQ=255)
                                  aaTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                   <  2:1202827/100‑1 (MQ=255)
                                  aaTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                   >  1:2733382/1‑100 (MQ=255)
                                  aaTCAAACCGAAGCCACATATGCGCGGCCAGAGTGTTGACAAAGGGCGCTTTATTCAGGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                   <  2:2329783/100‑1 (MQ=255)
                                     cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCTTGAACGCTTTATTCGGTCTACACAAGc                                                                     >  2:3588975/1‑95 (MQ=255)
                                     cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCa                                                                >  2:555288/1‑100 (MQ=255)
                                      aaaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCaa                                                               >  1:1687891/1‑100 (MQ=255)
                                       aaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCaaa                                                              >  1:2789812/1‑100 (MQ=255)
                                         ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                            <  1:14003/100‑1 (MQ=255)
                                         ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                            >  1:348520/1‑100 (MQ=255)
                                         ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                            <  1:2520942/100‑1 (MQ=255)
                                          cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                           <  2:583744/100‑1 (MQ=255)
                                             aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCTGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                        <  2:2289764/100‑1 (MQ=255)
                                             aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCTGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCATGCAAATTCAat                                                        <  2:742690/100‑1 (MQ=255)
                                             aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                        <  2:2067847/100‑1 (MQ=255)
                                              gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                                       <  1:454880/100‑1 (MQ=255)
                                                 acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCCGGCAAATTAAatatat                                                    >  2:3477080/1‑100 (MQ=255)
                                                  caTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCTGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                   >  2:1045831/1‑100 (MQ=255)
                                                  caTATGCCCGGCCAGACTGTTGACAAAGGGCGCTTTGTTCATGCAGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTAAAGATAtt                                                   >  2:3522753/1‑100 (MQ=25)
                                                    tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                                 <  2:2092855/100‑1 (MQ=255)
                                                     atGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                                <  2:2733364/100‑1 (MQ=255)
                                                       gcgcGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCaga                                              >  1:546719/1‑100 (MQ=255)
                                                        cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                             >  1:1333488/1‑100 (MQ=255)
                                                            gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                         <  1:3335712/100‑1 (MQ=255)
                                                            gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                         <  2:663480/100‑1 (MQ=255)
                                                             ccAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatgat                                        >  1:1266179/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGatgatg                                       <  2:2128313/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGatgatg                                       <  2:3216531/98‑1 (MQ=255)
                                                                gATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTa                                     <  1:2012188/100‑1 (MQ=255)
                                                                 aTTGTTGACAAAGGGCGTTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAg                                    <  2:428098/100‑1 (MQ=255)
                                                                   tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATCTTGCAGAGCTGATGTAGa                                   >  1:1576129/1‑99 (MQ=255)
                                                                   tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGac                                  >  1:409651/1‑100 (MQ=255)
                                                                       gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                              >  1:684800/1‑100 (MQ=255)
                                                                       gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                              >  1:2805204/1‑100 (MQ=255)
                                                                        aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTgg                             <  2:1910884/100‑1 (MQ=255)
                                                                         cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                            <  1:1718978/100‑1 (MQ=255)
                                                                                 gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAg                    >  2:975533/1‑100 (MQ=255)
                                                                                   tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg                  <  2:2080030/100‑1 (MQ=255)
                                                                                   agcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg                  >  1:2243745/4‑100 (MQ=37)
                                                                                   agcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg                  >  2:819861/4‑100 (MQ=255)
                                                                                     tGTTCATGCCGGATGCGGCATGAACGCTGTATGCGGTCTACAAAGGCAGTGAAATTCAATATGTTGCAGAGATGATGTAGAGACTGGCAAGCGTAGCGCa                <  2:1255732/100‑1 (MQ=21)
                                                                                     tGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcgc                 >  2:926151/1‑99 (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:57817/1‑100 (MQ=255)
                                                                                          aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg           <  2:348517/100‑1 (MQ=255)
                                                                                           tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt          >  2:3127761/1‑100 (MQ=255)
                                                                                            gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa         >  1:184746/1‑100 (MQ=255)
                                                                                            gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa         >  1:404325/1‑100 (MQ=255)
                                                                                             ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTaa        >  1:1366421/1‑100 (MQ=255)
                                                                                                  tACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTg   <  1:418478/100‑1 (MQ=255)
                                                                                                   aCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGc  <  1:2949077/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGC  >  NC_000913/770583‑770781

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: