Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 839,256 G→A P500L (CCG→CTG)  b0805 ← putative outer membrane receptor for iron transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913839,2560GA78.7% 577.1 / 119.1 253P500L (CCG→CTG) b0805putative outer membrane receptor for iron transport
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (96/103);  total (127/126)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 260 total aligned reads displayed.

TATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTA  >  NC_000913/839164‑839348
                                                                                            |                                                                                            
tataGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCg                                                                                              >  1:5719720/1‑93 (MQ=255)
tataGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCg                                                                                              >  2:7341406/1‑93 (MQ=255)
tataGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCg                                                                                              >  2:6570667/1‑93 (MQ=255)
tataGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCg                                                                                              <  2:1091275/93‑1 (MQ=255)
  taGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAga                                                                                            <  1:1060846/93‑1 (MQ=255)
  taGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAga                                                                                            >  1:1731056/1‑93 (MQ=255)
    gACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAgaga                                                                                          >  2:1889120/1‑93 (MQ=255)
     aCATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGaa                                                                                         <  2:2007723/93‑1 (MQ=255)
      cATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAAc                                                                                        >  1:7528156/1‑93 (MQ=255)
       aTTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAAcc                                                                                       >  2:3968095/1‑93 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      >  1:5306274/1‑93 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      <  1:6290927/93‑1 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      <  1:6735315/93‑1 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      <  2:1376697/93‑1 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      <  2:68117/93‑1 (MQ=255)
        ttGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCt                                                                                      <  1:1233035/93‑1 (MQ=255)
         tGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCtt                                                                                     <  1:21238/93‑1 (MQ=255)
          gCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCttt                                                                                    >  2:2639064/1‑93 (MQ=255)
          gCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCttt                                                                                    <  2:6066622/93‑1 (MQ=255)
          gCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCttt                                                                                    <  2:1255230/93‑1 (MQ=255)
           ccGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCtttt                                                                                   >  1:1935140/1‑93 (MQ=255)
            cGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTg                                                                                  <  2:1526162/93‑1 (MQ=255)
               tttCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCca                                                                               <  1:3880916/93‑1 (MQ=255)
               tttCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCca                                                                               <  1:7020667/93‑1 (MQ=255)
                    gTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCa                                                                          <  1:3909705/93‑1 (MQ=255)
                     tCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAg                                                                         >  2:3563791/1‑93 (MQ=255)
                     tCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAg                                                                         >  2:498744/1‑93 (MQ=255)
                       aGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGtt                                                                       <  2:4267995/93‑1 (MQ=255)
                       aGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGtt                                                                       <  2:3194626/93‑1 (MQ=255)
                        ggTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTg                                                                      <  1:5907681/93‑1 (MQ=255)
                        ggTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTg                                                                      >  2:6017314/1‑93 (MQ=255)
                        ggTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTg                                                                      >  1:5416731/1‑93 (MQ=255)
                         gTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTgg                                                                     >  2:5456452/1‑93 (MQ=255)
                         gTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTgg                                                                     <  2:6373962/93‑1 (MQ=255)
                           gATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGc                                                                   >  1:4735542/1‑93 (MQ=255)
                            aTACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCa                                                                  >  1:1570963/1‑93 (MQ=255)
                            aTACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCa                                                                  <  2:3084559/93‑1 (MQ=255)
                            aTACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCa                                                                  >  1:5119740/1‑93 (MQ=255)
                             tACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAg                                                                 >  1:1286016/1‑93 (MQ=255)
                             tACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAg                                                                 >  1:1545037/1‑93 (MQ=255)
                               cAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGt                                                               >  2:2468750/1‑93 (MQ=255)
                                 gcgcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGa                                                             <  1:5772472/93‑1 (MQ=255)
                                 gcgcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGa                                                             <  1:5554430/93‑1 (MQ=255)
                                 gcgcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGa                                                             >  2:305416/1‑93 (MQ=255)
                                  cgcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGAt                                                            <  1:6584470/93‑1 (MQ=255)
                                  cgcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGAt                                                            <  2:2373623/93‑1 (MQ=255)
                                   gcCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATg                                                           >  1:3252440/1‑93 (MQ=255)
                                    ccccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATgg                                                          >  1:6923740/1‑93 (MQ=255)
                                     cccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGc                                                         >  2:3885143/1‑93 (MQ=255)
                                     cccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGc                                                         >  2:5391448/1‑93 (MQ=255)
                                      ccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCa                                                        >  2:6197785/1‑93 (MQ=255)
                                      ccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCa                                                        >  2:5132654/1‑93 (MQ=255)
                                      ccGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCa                                                        >  2:3432450/1‑93 (MQ=255)
                                        ggCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCAcc                                                      >  1:5613373/1‑93 (MQ=255)
                                         gCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACcg                                                     >  2:1347327/1‑93 (MQ=255)
                                         gCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACcg                                                     >  1:7460916/1‑93 (MQ=255)
                                           tttCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACcgcg                                                   >  1:3775737/1‑93 (MQ=255)
                                             tCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTc                                                 <  1:7216436/93‑1 (MQ=255)
                                              ccAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTcc                                                <  2:4552698/93‑1 (MQ=255)
                                               cAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCg                                               >  1:2513835/1‑93 (MQ=255)
                                                 gTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCt                                             <  2:2200258/93‑1 (MQ=255)
                                                 gTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCt                                             >  1:3470716/1‑93 (MQ=255)
                                                  ttCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTg                                            >  2:1625865/1‑93 (MQ=255)
                                                  ttCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTg                                            >  1:2318451/1‑93 (MQ=255)
                                                      tcaGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTgccgc                                        >  1:2578053/1‑93 (MQ=255)
                                                      tcaGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTgccgc                                        >  2:6669681/1‑93 (MQ=255)
                                                       caGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTgccgcc                                       <  2:4757339/93‑1 (MQ=255)
                                                        aGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGccgccg                                      >  1:2852424/1‑93 (MQ=255)
                                                         gATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCcgccgc                                     <  2:6270241/93‑1 (MQ=255)
                                                          aTTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCa                                    >  2:2486006/1‑93 (MQ=255)
                                                          aTTGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCa                                    <  1:4084013/93‑1 (MQ=255)
                                                           ttGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAg                                   >  1:3327797/1‑93 (MQ=255)
                                                           ttGCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAg                                   <  1:5230211/93‑1 (MQ=255)
                                                             gCCCGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGc                                 <  2:4190795/93‑1 (MQ=255)
                                                                cGACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGt                              <  1:3028182/93‑1 (MQ=255)
                                                                 gACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTg                             <  2:5602888/93‑1 (MQ=255)
                                                                 gACTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTg                             <  1:6495975/93‑1 (MQ=255)
                                                                  aCTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTgg                            <  1:1027125/93‑1 (MQ=255)
                                                                  aCTTGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTgg                            >  1:3005706/1‑93 (MQ=255)
                                                                    ttGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCa                          <  2:859354/93‑1 (MQ=255)
                                                                    ttGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCa                          >  2:3714380/1‑93 (MQ=255)
                                                                     tGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCAc                         <  2:5895406/93‑1 (MQ=255)
                                                                     tGGCGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCAc                         >  1:6159461/1‑93 (MQ=255)
                                                                        cGGTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGt                      <  1:4459052/93‑1 (MQ=255)
                                                                         ggTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTc                     >  2:2847796/1‑93 (MQ=255)
                                                                         ggTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTc                     <  1:1335579/93‑1 (MQ=255)
                                                                         ggTGTCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTc                     >  2:332156/1‑93 (MQ=255)
                                                                            gtCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCata                  <  1:2037464/93‑1 (MQ=255)
                                                                            gtCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCata                  <  2:2049518/93‑1 (MQ=255)
                                                                             tCGACGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCatat                 >  1:4316121/1‑93 (MQ=255)
                                                                                gcGGTGGTGACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCa              <  2:7088811/92‑1 (MQ=255)
                                                                                      gtgACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGa        <  2:5313042/93‑1 (MQ=255)
                                                                                       tgACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGAt       <  2:5217695/93‑1 (MQ=255)
                                                                                        gACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATa      >  2:3267539/1‑93 (MQ=255)
                                                                                        gACCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATa      <  2:1740069/93‑1 (MQ=255)
                                                                                         aCCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATaa     <  2:792736/93‑1 (MQ=255)
                                                                                         aCCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATaa     <  2:4333512/93‑1 (MQ=255)
                                                                                         aCCAGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATaa     <  1:3673474/93‑1 (MQ=255)
                                                                                            ggAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTa  <  1:4408478/93‑1 (MQ=255)
                                                                                            ggAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTa  >  1:3373114/1‑93 (MQ=255)
                                                                                            |                                                                                            
TATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTA  >  NC_000913/839164‑839348

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