Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→6 coding (1279/2526 nt) ecpC ← ECP production outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,0480C.97.5% 355.0 / 1.9 81coding (1285/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (36/43);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307950‑308143
                                                                                                  |                                                                                               
caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcc                                                                                                >  2:1016503/1‑100 (MQ=255)
 aGCGAACTAAAGCCTCCCGGTAGAGTGTCGCTGATGCTGCCGAGGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccc                                                                                               <  1:2721844/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  2:2079427/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  2:3199199/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  1:1387154/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  1:2872564/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  2:3136837/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  1:307065/1‑100 (MQ=255)
    gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccc                                                                                            >  2:722563/1‑100 (MQ=255)
     aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccca                                                                                           <  1:45063/100‑2 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaag                                                                                         >  1:1152497/1‑97 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaag                                                                                         >  2:223715/1‑97 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccca                                                                                           >  2:1594392/1‑97 (MQ=255)
        tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        <  1:2624614/100‑5 (MQ=255)
        tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        >  1:727062/1‑96 (MQ=255)
         aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCg                                                                                      >  1:94379/1‑100 (MQ=255)
           aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGCGCTGTCACTGGCCAGCATATTTTGCAGGATAACGTTGATCGCCCCCCAAGCGGc                                                                                    >  2:849412/1‑100 (MQ=255)
             ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAg                                                                                  <  2:2433854/100‑1 (MQ=255)
             ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATACTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAg                                                                                  >  1:361087/1‑100 (MQ=255)
              cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGc                                                                                 >  1:3158292/1‑100 (MQ=255)
              cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGc                                                                                 >  1:2256440/1‑100 (MQ=255)
                nccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGt                                                                               >  1:2081523/2‑100 (MQ=255)
                 ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTc                                                                              >  1:2776022/1‑100 (MQ=255)
                 ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTc                                                                              <  1:2389689/100‑1 (MQ=255)
                     tAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGAc                                                                          <  2:470028/100‑1 (MQ=255)
                      agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGAc                                                                          <  1:3022842/100‑1 (MQ=255)
                        ngTGGCGCTGATTCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCTGCAGCGTCAGACggg                                                                       >  2:936893/2‑100 (MQ=255)
                              gcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcac                                                                 <  2:1136127/100‑1 (MQ=255)
                               cTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacc                                                                <  1:2023492/100‑1 (MQ=255)
                                tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccc                                                               <  2:557467/100‑1 (MQ=255)
                                 gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccca                                                              >  1:1846689/1‑100 (MQ=255)
                                     ctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                          <  2:2084984/100‑1 (MQ=255)
                                      tgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                         <  1:129905/100‑1 (MQ=255)
                                      tgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                         >  1:3119822/1‑100 (MQ=255)
                                       gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCt                                                        <  1:1622156/100‑1 (MQ=255)
                                        ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  2:651435/100‑1 (MQ=255)
                                        ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  2:1846667/100‑1 (MQ=255)
                                         cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGa                                                      >  1:693388/1‑100 (MQ=255)
                                            tgctgcTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTa                                                   >  2:3266835/1‑100 (MQ=255)
                                             gctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                  >  1:1337692/1‑100 (MQ=255)
                                              ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 <  1:592235/100‑1 (MQ=255)
                                              ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 >  1:128525/1‑100 (MQ=255)
                                               tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                <  2:2642646/100‑1 (MQ=255)
                                               tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                >  1:1861611/1‑100 (MQ=255)
                                                gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                               >  1:3211936/1‑100 (MQ=255)
                                                 ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                              <  1:1009540/100‑1 (MQ=255)
                                                 ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                              <  1:2454513/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  1:491829/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  2:2535312/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  1:3045137/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:345802/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:1635938/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:217941/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        <  2:2372731/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        <  2:2497221/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        >  1:2318515/1‑100 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  2:306372/100‑1 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  1:1016515/100‑1 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  2:3251078/100‑1 (MQ=255)
                                                           tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                    <  2:2074855/100‑1 (MQ=255)
                                                           tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                    >  2:1074206/1‑100 (MQ=255)
                                                            gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                   >  1:733817/1‑100 (MQ=255)
                                                             tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGt                                  >  2:269496/1‑100 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCCTATCCCTATCCCg                                   >  2:1626623/1‑99 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                 <  2:1637243/100‑1 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                 >  1:3219029/1‑100 (MQ=255)
                                                               aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                <  1:276624/100‑1 (MQ=255)
                                                               aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                >  1:1369782/1‑100 (MQ=255)
                                                                cTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATCATCGTATCCATATCCCGTTgc                               <  1:2407540/100‑1 (MQ=255)
                                                                atcggcAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgc                               <  2:1220548/95‑1 (MQ=255)
                                                                   gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgc                            >  1:1675421/1‑100 (MQ=255)
                                                                    ccAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgcc                           <  1:1594409/100‑1 (MQ=255)
                                                                    ccAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                             >  2:233771/1‑98 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          <  2:2256414/100‑1 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          >  2:407979/1‑100 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          <  1:134316/100‑1 (MQ=255)
                                                                       gCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                        >  2:1491277/1‑100 (MQ=255)
                                                                           atTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                    >  1:1740457/1‑100 (MQ=255)
                                                                            ttttGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                   <  2:2882011/100‑1 (MQ=255)
                                                                             cttccgatctAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                  <  1:943751/91‑1 (MQ=255)
                                                                                 cAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              >  1:86954/1‑100 (MQ=255)
                                                                                 cAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              >  2:1888520/1‑100 (MQ=255)
                                                                                  aGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACt             <  1:177441/100‑1 (MQ=255)
                                                                                   ggTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc            >  2:1034942/1‑100 (MQ=255)
                                                                                    gTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa           <  2:1391734/100‑1 (MQ=255)
                                                                                      tAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt         <  2:588745/100‑1 (MQ=255)
                                                                                       aaCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGTTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTTAGCGTACTCa           >  2:1385957/1‑97 (MQ=255)
                                                                                       aaCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGcagatcggaag        >  2:943739/1‑90 (MQ=255)
                                                                                       aaCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg        <  1:3136874/100‑1 (MQ=255)
                                                                                        aCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       >  2:1129614/1‑100 (MQ=255)
                                                                                         cGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       <  2:1547487/99‑1 (MQ=255)
                                                                                         cGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg      >  1:3170514/1‑100 (MQ=255)
                                                                                          gTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc     >  1:2924300/1‑100 (MQ=255)
                                                                                          gTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc     <  1:425718/100‑1 (MQ=255)
                                                                                            tGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg   >  1:2102095/1‑100 (MQ=255)
                                                                                             gATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg  <  2:727054/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307950‑308143

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