Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 615,736 A→C E527A (GAA→GCA)  entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913615,7360AC98.2% 562.9 / ‑1.5 162E527A (GAA→GCA) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (87/72);  total (87/75)
Fisher's exact test for biased strand distribution p-value = 9.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 163 total aligned reads displayed.

TCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCAC  >  NC_000913/615641‑615832
                                                                                               |                                                                                                
tCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCt                                                                                              >  2:3753598/1‑100 (MQ=255)
tCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCt                                                                                              <  1:1665094/100‑1 (MQ=255)
 ctgctgCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGagatcg                                                                                             >  1:733866/1‑95 (MQ=255)
 ctgctgCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTg                                                                                             >  1:3148592/1‑100 (MQ=255)
  tgctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTgg                                                                                            <  1:1194578/100‑1 (MQ=255)
  tgctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTgg                                                                                            >  1:288402/1‑100 (MQ=255)
  tgctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTgg                                                                                            >  1:2957081/1‑100 (MQ=255)
   gctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGa                                                                                           <  2:74301/100‑1 (MQ=255)
   gctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGa                                                                                           >  2:2079784/1‑100 (MQ=255)
    ctgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAg                                                                                          <  1:2797770/100‑1 (MQ=255)
     tgcGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGc                                                                                         <  1:3057807/100‑1 (MQ=255)
        gTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGc                                                                                      >  1:3500888/1‑100 (MQ=255)
        gTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGc                                                                                      >  2:3091965/1‑100 (MQ=255)
         tGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGcc                                                                                     <  1:85682/100‑1 (MQ=255)
         tGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGcc                                                                                     >  1:1350615/1‑100 (MQ=255)
         tGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGcc                                                                                     >  1:2241898/1‑100 (MQ=255)
          gAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCt                                                                                    <  1:3050379/100‑1 (MQ=255)
           aGCGCGGCGTTAAACCAGGGGGCAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTg                                                                                   <  1:360216/100‑1 (MQ=255)
            gcgcGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTgg                                                                                  <  1:78338/100‑1 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 <  2:2361058/100‑1 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 >  1:374213/1‑100 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 >  2:2898136/1‑100 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 <  2:929379/100‑1 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 <  1:2349577/100‑1 (MQ=255)
             cgcgGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACGACCGCGTTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                 <  1:1582833/100‑1 (MQ=255)
               cggcgTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTa                                                                               <  2:1616965/100‑1 (MQ=255)
                ggcgTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTAc                                                                              <  2:2439545/100‑1 (MQ=255)
                 gcgTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTAcc                                                                             >  2:3159849/1‑100 (MQ=255)
                   gTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGc                                                                           >  1:812411/1‑100 (MQ=255)
                   gTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGc                                                                           <  1:434152/100‑1 (MQ=255)
                    ttAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCt                                                                          <  1:3942139/100‑1 (MQ=255)
                     tAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTg                                                                         <  2:73186/100‑1 (MQ=255)
                     tAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTg                                                                         <  2:2680609/100‑1 (MQ=255)
                     caaaCCCGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTg                                                                         <  1:1204183/99‑1 (MQ=255)
                        aCCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGAt                                                                      <  2:3300449/100‑1 (MQ=255)
                        aCCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGAt                                                                      <  2:368270/100‑1 (MQ=255)
                         ccAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATa                                                                     >  1:1926401/1‑100 (MQ=255)
                         ccAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATa                                                                     >  1:1755030/1‑100 (MQ=255)
                         ccAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATa                                                                     >  1:563434/1‑100 (MQ=255)
                         ccAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATa                                                                     >  1:2006368/1‑100 (MQ=255)
                            ggggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCg                                                                  >  2:2656714/1‑100 (MQ=255)
                            ggggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCg                                                                  >  2:3549406/1‑100 (MQ=255)
                            ggggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCg                                                                  <  1:2582426/100‑1 (MQ=255)
                             gggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCgg                                                                 <  2:3016179/100‑1 (MQ=255)
                              ggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGc                                                                <  1:3191637/100‑1 (MQ=255)
                                aCAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTa                                                              <  1:408531/100‑1 (MQ=255)
                                    cGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCg                                                          <  1:3427798/100‑1 (MQ=255)
                                      tGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGa                                                        <  2:3801022/100‑1 (MQ=255)
                                          ggTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGAt                                                    >  2:3646884/1‑100 (MQ=255)
                                            tGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCg                                                  <  2:3589711/100‑1 (MQ=255)
                                               cACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCt                                               <  2:1510386/100‑1 (MQ=255)
                                                aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTg                                              >  2:917444/1‑100 (MQ=255)
                                                aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTg                                              <  1:2214441/100‑1 (MQ=255)
                                                 cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGa                                             <  2:3543113/100‑1 (MQ=255)
                                                 cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGa                                             >  2:1835775/1‑100 (MQ=255)
                                                 cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGa                                             >  1:2444545/1‑100 (MQ=255)
                                                 cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGa                                             >  2:2732053/1‑100 (MQ=255)
                                                  tACCGCTCTCGGTCTTTTAGACCCTGGCACTCCATACGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTTGATACCTGCTATCCGGACGATCGcataaa                                            >  1:3398678/1‑95 (MQ=255)
                                                  naCCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaa                                            >  1:2015421/2‑100 (MQ=255)
                                                   aCCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaaa                                           >  1:1863607/1‑100 (MQ=255)
                                                    ccGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaaaa                                          <  1:2406174/100‑1 (MQ=255)
                                                    ccGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaaaa                                          <  2:404258/100‑1 (MQ=255)
                                                     cgcgGTCGGTCTTTTTGACCCAGGCATTACAAGCGATATTTTCAGCTGAAGCAGAATGTGTACCGCTATATACCCGCTATCTGGACGACCGCCTGAAAat                                         >  2:1322295/1‑100 (MQ=255)
                                                     cgcgCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAat                                         <  2:1929667/100‑1 (MQ=255)
                                                       cgcTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAatga                                       <  2:38991/100‑1 (MQ=255)
                                                        gcTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAatgat                                      <  1:361525/100‑1 (MQ=255)
                                                         cTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAatgatg                                     >  2:364450/1‑100 (MQ=255)
                                                           ctgTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCt                                   <  1:1007567/98‑1 (MQ=255)
                                                           cGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCt                                   <  1:1202529/100‑1 (MQ=255)
                                                            ggTCTTTTTGACCGTGGCACTACATGAGATAGTTGCAGTTGGTGTGGACTGGCTACCGCGGCGTACCGGCTATCCGGAGGATCGGCTGAAAATTAAGCTg                                  >  2:161859/1‑100 (MQ=255)
                                                             gTCTTTTTTACCCTTTTACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGACAATGATGCTgg                                 >  2:1598090/1‑100 (MQ=255)
                                                             gTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTgg                                 >  1:1663976/1‑100 (MQ=255)
                                                             gTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTgg                                 >  1:2097847/1‑100 (MQ=255)
                                                             gTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTgg                                 >  2:1152280/1‑100 (MQ=255)
                                                             gTCTATTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTgg                                 >  2:3399606/1‑100 (MQ=255)
                                                               cTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGCCGATCGCCTGCACATGATGCtcgca                               >  2:3339164/1‑96 (MQ=255)
                                                               cTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGaa                               >  2:918765/1‑100 (MQ=255)
                                                               cTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGaa                               >  1:692740/1‑100 (MQ=255)
                                                                ttttttACTCTTGCACTCCATGCGATAGTTGCAGCTGGAGCGGCTTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAg                              >  2:3344717/1‑100 (MQ=255)
                                                                     gACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATgcg                         >  2:1782676/1‑100 (MQ=255)
                                                                     gACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATgcg                         <  2:1062827/100‑1 (MQ=255)
                                                                       cccTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGcgcg                       >  1:1301834/1‑100 (MQ=255)
                                                                        ccTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCGGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgt                      <  2:3274427/100‑1 (MQ=255)
                                                                          tGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCATGAAAATGATTCTGTAAGATGCAcgtcc                    >  2:3130145/1‑100 (MQ=255)
                                                                            gCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgt                  <  2:1420140/100‑1 (MQ=255)
                                                                            gCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgt                  <  1:252319/100‑1 (MQ=255)
                                                                             cACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgtc                 <  2:1234436/100‑1 (MQ=255)
                                                                             cACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgtc                 >  2:2753824/1‑100 (MQ=255)
                                                                             cACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgtc                 >  2:2145968/1‑100 (MQ=255)
                                                                              aCTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCg                >  1:906399/1‑100 (MQ=255)
                                                                               cTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGc               >  1:1255821/1‑100 (MQ=255)
                                                                                 ccATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTg             <  2:3032267/100‑1 (MQ=255)
                                                                                 ccATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTg             >  1:2827178/1‑100 (MQ=255)
                                                                                  cATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCATCCGTCGCTGt            <  2:1833654/100‑1 (MQ=255)
                                                                                      cGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAAt        <  1:2688933/100‑1 (MQ=255)
                                                                                       gATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAAtt       <  2:1305715/100‑1 (MQ=255)
                                                                                        aTAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTa      <  2:3362808/100‑1 (MQ=255)
                                                                                        aTAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTa      <  2:2624819/100‑1 (MQ=255)
                                                                                         tAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTac     >  2:501971/1‑100 (MQ=255)
                                                                                           gTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTacca   >  2:66278/1‑100 (MQ=255)
                                                                                               |                                                                                                
TCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCAC  >  NC_000913/615641‑615832

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