Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,679 G→T 15.0% 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,6790GT15.0% 422.0 / 57.2 154intergenic (+68/‑779)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (54/76);  new base T (14/9);  total (68/86)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

Only 100 of 157 total aligned reads displayed.

AGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGT  >  NC_000913/770580‑770773
                                                                                                   |                                                                                              
aGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                <  1:792068/100‑1 (MQ=255)
 gTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGa                                                                                               >  1:1887968/1‑100 (MQ=255)
 gTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGa                                                                                               <  1:806682/100‑1 (MQ=255)
     cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                           >  2:1424329/1‑100 (MQ=255)
       gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                         <  1:1563242/100‑1 (MQ=255)
         cGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCAt                                                                                       <  1:1395402/100‑1 (MQ=255)
          gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATg                                                                                      <  2:2395931/100‑1 (MQ=255)
          gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATg                                                                                      <  2:1602381/100‑1 (MQ=255)
           ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                     <  2:346943/100‑1 (MQ=255)
            tCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGACATGa                                                                                     <  1:3238276/99‑1 (MQ=255)
            tCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGaa                                                                                    <  2:410696/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGACGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                   <  1:448445/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                   <  1:2582173/100‑1 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAAc                                                                                   >  1:1010365/1‑100 (MQ=255)
             cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGTATGCGGCATGAAc                                                                                   <  1:916752/100‑1 (MQ=255)
              aGTTACCGGCTTGTCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                   <  2:1335346/99‑1 (MQ=255)
              aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATGTGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                   <  1:1013895/99‑1 (MQ=255)
              aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACg                                                                                  >  2:126257/1‑100 (MQ=255)
              aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACg                                                                                  <  1:2615335/100‑1 (MQ=255)
               gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                   <  2:3327273/98‑1 (MQ=255)
               gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                 >  2:2290133/1‑100 (MQ=255)
                 tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                               <  2:2152662/100‑1 (MQ=255)
                 tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                               <  2:303985/100‑1 (MQ=255)
                   ccGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTa                                                                             <  2:2105586/100‑1 (MQ=255)
                    cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                            <  1:2154356/100‑1 (MQ=255)
                     ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAtt                                                                           >  2:86621/1‑100 (MQ=255)
                     ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAtt                                                                           >  1:346307/1‑100 (MQ=255)
                      gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTc                                                                          <  2:1914780/100‑1 (MQ=255)
                      gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTc                                                                          >  2:2760515/1‑100 (MQ=255)
                         tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                       >  2:2281323/1‑100 (MQ=255)
                         tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                       <  1:3133543/100‑1 (MQ=255)
                          gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                      <  2:1243227/100‑1 (MQ=255)
                          gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                      <  1:2045662/100‑1 (MQ=255)
                           ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                       <  2:146907/98‑1 (MQ=255)
                           ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                     <  2:3375837/100‑1 (MQ=255)
                           ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                     >  1:735070/1‑100 (MQ=255)
                            cTGAATAGCAATCAAACCTAAGCCACATAGGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTa                                                                    <  1:176913/100‑1 (MQ=255)
                             tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTAc                                                                   >  2:3271766/1‑100 (MQ=255)
                                 tAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaaa                                                               <  2:1383673/100‑1 (MQ=255)
                                  aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                              <  1:682963/100‑1 (MQ=255)
                                  aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                              <  2:2095958/100‑1 (MQ=255)
                                   gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                             >  2:2457621/1‑100 (MQ=255)
                                       tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                         <  2:3270250/100‑1 (MQ=255)
                                        cAAACCGACGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGAATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCa                                                        <  1:3353440/100‑1 (MQ=255)
                                        cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATTCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCa                                                        <  1:1816130/100‑1 (MQ=255)
                                        cAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCa                                                        >  2:2978302/1‑100 (MQ=255)
                                             cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                   <  2:2158764/100‑1 (MQ=255)
                                             cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                   <  2:1887980/100‑1 (MQ=255)
                                              gAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCa                                                  >  2:61357/1‑100 (MQ=255)
                                                aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                <  2:48392/100‑1 (MQ=255)
                                                aGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                <  1:3159060/100‑1 (MQ=255)
                                                 gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                               <  2:1074708/100‑1 (MQ=255)
                                                 gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                               <  2:1581131/100‑1 (MQ=255)
                                                       tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                         >  1:1096774/1‑100 (MQ=255)
                                                           cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                     <  1:956158/100‑1 (MQ=255)
                                                           cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                     <  1:39477/100‑1 (MQ=255)
                                                           cgcgGCCAGATTGTTGACAAAGAGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                     >  1:1965342/1‑100 (MQ=255)
                                                             cgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGATTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGTgat                                   <  1:1262569/100‑1 (MQ=255)
                                                              ggCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAAGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                   <  1:2829959/99‑1 (MQ=255)
                                                              gaccAGATTGTTGATAAAGGGCGCTTTTTGCATGCCGGATACTGCATGATCGCTTTATTCGCGCTACAAACGCATGCAAATTCAATATATTGCagag                                     <  1:2141893/95‑1 (MQ=18)
                                                                  aGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGt                              <  1:61357/100‑1 (MQ=255)
                                                                   gATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTa                             <  1:2998580/100‑1 (MQ=255)
                                                                   gATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTa                             <  2:1791206/100‑1 (MQ=255)
                                                                    aTTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAg                            <  1:1180116/100‑1 (MQ=255)
                                                                     ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGa                           <  2:2038953/100‑1 (MQ=255)
                                                                     ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGa                           <  1:2814319/100‑1 (MQ=255)
                                                                     ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGa                           <  1:2626430/100‑1 (MQ=255)
                                                                      tgttgACTAAGGGCGCTTTGTTCATGCCGGATAGGGCATGAACTTTTTATTCGGTCTACAATAGCATTCAAATTCTATATATTGCAGAGATGATGTAGac                          <  1:830081/100‑1 (MQ=21)
                                                                       gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                         >  1:2451649/1‑100 (MQ=255)
                                                                         tgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACt                       >  2:1130251/1‑100 (MQ=255)
                                                                         tgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACt                       >  2:1253814/1‑100 (MQ=255)
                                                                           aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTgg                     >  2:1867834/1‑100 (MQ=255)
                                                                            cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                    <  2:3110925/100‑1 (MQ=255)
                                                                             aaaGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCa                   <  1:1406548/100‑1 (MQ=255)
                                                                               aGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAg                 <  2:767485/100‑1 (MQ=255)
                                                                                 ggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTCGCAAGCg               <  2:1773959/100‑1 (MQ=255)
                                                                                 ggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTAAATATATTGCAGAGATGATGTAGACACTGGCAAGCg               >  2:1119879/1‑100 (MQ=255)
                                                                                  ataaagcgTTCATGCCGTATCCGGCATGAACACTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt              >  2:2841392/8‑100 (MQ=25)
                                                                                   cgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTa             <  2:2291820/100‑1 (MQ=255)
                                                                                     gtttGTTCATGCCGGATACGCCATGAACGCTTTATTCGGACTATAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc           <  2:1419125/99‑1 (MQ=255)
                                                                                     cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc           <  2:2935049/100‑1 (MQ=255)
                                                                                     aagcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc           >  2:1086864/5‑100 (MQ=37)
                                                                                        tGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  2:3299749/1‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  2:2919036/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  1:712688/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  2:2436857/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  2:737735/2‑100 (MQ=255)
                                                                                        cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  1:2423055/2‑100 (MQ=255)
                                                                                         gTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt       >  2:2559758/1‑100 (MQ=255)
                                                                                         gTTCATGCCGTATCCGGCATGAACACTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt       >  1:1657377/1‑100 (MQ=37)
                                                                                         gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt       >  1:1617481/1‑100 (MQ=255)
                                                                                         gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATACATTGCATAGATGATGTAGACACTGGCAAGCGTAGCGCAt       >  1:2381092/1‑100 (MQ=255)
                                                                                          ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATc      <  2:945979/100‑1 (MQ=255)
                                                                                            cATGCCGTATCCGGCATGAACGCTTTATTCGGTGTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGAAGCGCATCAg    >  1:2231505/1‑100 (MQ=37)
                                                                                            cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa        >  1:2665158/1‑96 (MQ=255)
                                                                                             aTGCCGGATACGGCATGAACGCTTTCTTCGCTCTACACAATCAGGCAACTTCAATATATCGCAGAGATGCTGTAGACACTGGCCAACGTCACACATCa     >  1:2671560/1‑98 (MQ=255)
                                                                                             aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg   >  2:2370738/1‑100 (MQ=255)
                                                                                              tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt  >  2:2044396/1‑100 (MQ=255)
                                                                                              tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt  >  2:715568/1‑100 (MQ=255)
                                                                                              tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt  <  1:1064233/100‑1 (MQ=255)
                                                                                                   |                                                                                              
AGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGT  >  NC_000913/770580‑770773

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: