Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,679 G→T 16.1% 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.1% 382.0 / 55.6 140intergenic (+68/‑779)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (49/68);  new base T (13/10);  total (62/78)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 144 total aligned reads displayed.

AGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTT  >  NC_000913/770580‑770778
                                                                                                   |                                                                                                   
aGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                     <  2:126714/100‑1 (MQ=255)
aGTCACGGACGTCCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                     <  2:706556/100‑1 (MQ=255)
    aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                 >  2:430769/1‑100 (MQ=255)
      ggACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACgg                                                                                               <  2:2535590/100‑1 (MQ=255)
       gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                              <  1:566672/100‑1 (MQ=255)
       gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                              <  2:2731571/100‑1 (MQ=255)
        aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                             <  2:1553310/100‑1 (MQ=255)
        aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                             <  2:2634361/100‑1 (MQ=255)
         cGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                             <  2:2657338/99‑1 (MQ=255)
         cGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCAt                                                                                            <  2:3276876/100‑1 (MQ=255)
           ttCAGTTACCGGCTTGGCTGAATATCAATCAAACCGATGCCACATATGCGCGGCCAGATTGTTGACAGAGGGCGCTTTGTTCATGCCGGATACTGCATGa                                                                                          <  1:2780845/100‑1 (MQ=38)
           ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                          >  2:3333222/1‑100 (MQ=255)
           cttccgatctggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGa                                                                                          <  1:1431647/90‑1 (MQ=38)
             cAGTTGCCGGCTTGCCGGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                        <  2:465385/100‑1 (MQ=38)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                        <  2:1874300/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                        <  2:74945/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                        <  2:833601/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGGTCATGCCGGATACGGCATGAAc                                                                                        <  1:2065339/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATACGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                        <  1:770224/100‑1 (MQ=255)
               gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                      >  1:3066141/1‑100 (MQ=255)
                ttACCTGCTTGCCTGAATAGCAATCAATCCTAAGCAACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGTATGAACg                                                                                       <  1:1434754/98‑1 (MQ=38)
                 tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTCGTGCATGCCGTATCCGGCATGAAcaaag                                                                                    <  1:1568113/100‑5 (MQ=38)
                  aCCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATACCGGATACGGCATGAACGCttt                                                                                   <  2:3099891/100‑1 (MQ=255)
                     ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGaagatcggaag                                                                                >  2:1431654/1‑91 (MQ=38)
                      gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAttn                                                                               <  2:2944020/100‑2 (MQ=255)
                      gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTc                                                                               <  1:168408/100‑1 (MQ=255)
                       cTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                              <  2:1281069/100‑1 (MQ=255)
                        ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                             <  1:2251297/100‑1 (MQ=255)
                        ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                             <  2:2997342/100‑1 (MQ=255)
                        ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                             <  1:1831020/100‑1 (MQ=255)
                         tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                              >  1:2712575/1‑98 (MQ=255)
                          gCCTGAATAGCAATCAAATCGAAGCCGCATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTTTTTGGt                                                                            >  1:188648/1‑99 (MQ=38)
                          gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                           <  2:3048107/100‑1 (MQ=255)
                          gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                           <  2:2937134/100‑1 (MQ=255)
                           ccTGAATAGCAAGCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGCTACGGCATGAACGCTTTATTCGGTc                                                                           >  1:1482768/1‑99 (MQ=255)
                            cTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTa                                                                         <  1:2872294/100‑1 (MQ=255)
                              gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                       <  2:233843/100‑1 (MQ=255)
                              gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                       >  2:386974/1‑100 (MQ=255)
                              gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                       <  1:109492/100‑1 (MQ=255)
                                  aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                   >  1:700187/1‑100 (MQ=255)
                                   gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                  >  1:2094442/1‑100 (MQ=255)
                                    cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 <  1:430767/100‑1 (MQ=255)
                                    cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 >  1:3039784/1‑100 (MQ=255)
                                    cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 <  2:2038673/100‑1 (MQ=255)
                                    cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 <  2:1484623/100‑1 (MQ=255)
                                    cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 >  2:1434592/1‑100 (MQ=255)
                                      aTCAAACCGAAGCCACATATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                 <  1:569455/98‑1 (MQ=255)
                                       tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                              <  1:1202085/100‑1 (MQ=255)
                                       tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                              <  2:3356877/100‑1 (MQ=255)
                                       tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                              <  1:2749931/100‑1 (MQ=255)
                                            ccGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCCGGCaaact                                                         <  1:2483681/100‑3 (MQ=255)
                                             cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                        <  2:2244521/100‑1 (MQ=255)
                                                aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                     <  1:2019710/100‑1 (MQ=255)
                                                 gCCACATATGCGCGGCCAGATTGTTGCCAAAGGGCGCTTTGTTCATGCCGGATACGGCATAAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                     <  1:1213318/99‑1 (MQ=255)
                                                 gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                                    >  1:923098/1‑100 (MQ=255)
                                                  ccACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAatat                                                   <  2:484329/100‑1 (MQ=255)
                                                      atatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                >  1:1451117/1‑99 (MQ=255)
                                                       tatGCGCGGCCAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                <  1:2231433/98‑1 (MQ=255)
                                                       tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAATGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                              >  1:2494082/1‑100 (MQ=255)
                                                        atGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                             <  2:80896/100‑1 (MQ=255)
                                                             cgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                        >  1:2118252/1‑100 (MQ=255)
                                                              ggCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatg                                       <  2:1056794/100‑1 (MQ=255)
                                                               gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                        >  1:1121803/1‑98 (MQ=255)
                                                                     ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGTCAAATTCAATATATTGCCGAGATGATGTAGa                                <  2:2389546/100‑1 (MQ=255)
                                                                       gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                              >  2:2797675/1‑100 (MQ=255)
                                                                       gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                              >  2:2624686/1‑100 (MQ=255)
                                                                       gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGac                               >  1:2713471/1‑99 (MQ=255)
                                                                       gttgACAAAGGGCGCTTTGTTCATACCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                              >  1:1237102/1‑100 (MQ=255)
                                                                          gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                           <  1:2933034/100‑1 (MQ=255)
                                                                            cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                         >  2:1309336/1‑100 (MQ=255)
                                                                            cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                         >  2:3433238/1‑100 (MQ=255)
                                                                                gaataaagcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGc                     >  2:1950705/10‑100 (MQ=25)
                                                                                 ggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTCTTCGGTCTACACAAGCAGGTAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGc                     <  1:1935182/99‑1 (MQ=255)
                                                                                  gcgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                   >  2:3231920/1‑100 (MQ=255)
                                                                                  gcgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                   >  2:3041178/1‑100 (MQ=255)
                                                                                    gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAg                 <  1:2948606/100‑1 (MQ=255)
                                                                                    gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAg                 <  1:3469118/100‑1 (MQ=255)
                                                                                     cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                >  1:1797944/1‑100 (MQ=255)
                                                                                     cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                >  2:1648085/1‑100 (MQ=255)
                                                                                     cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                <  1:582380/100‑1 (MQ=255)
                                                                                      tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg               >  2:2683230/1‑100 (MQ=255)
                                                                                      tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg               <  2:2854959/100‑1 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAACAGCAGGCAACTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa             >  1:3451008/2‑100 (MQ=25)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa             >  2:2727319/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa             >  2:625994/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa             >  2:197779/2‑100 (MQ=255)
                                                                                         gTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt            >  1:2171622/1‑100 (MQ=255)
                                                                                         gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt            <  1:2617753/100‑1 (MQ=255)
                                                                                         gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt            >  2:463384/1‑100 (MQ=255)
                                                                                           tCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCa          >  2:2793332/1‑100 (MQ=255)
                                                                                            cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg         <  1:739711/100‑1 (MQ=255)
                                                                                             aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCa          <  1:1131292/98‑1 (MQ=255)
                                                                                             aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg        >  2:2683893/1‑100 (MQ=255)
                                                                                             aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg        >  2:1768396/1‑100 (MQ=255)
                                                                                              tGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt       <  2:199341/100‑1 (MQ=255)
                                                                                              tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt       >  2:3011725/1‑100 (MQ=255)
                                                                                               gCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa      >  1:967572/1‑100 (MQ=255)
                                                                                               gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa      <  1:3333198/100‑1 (MQ=255)
                                                                                                 cGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAt    <  1:3544576/100‑1 (MQ=255)
                                                                                                   gATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAttt  <  1:3268286/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTT  >  NC_000913/770580‑770778

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: