Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,245,958 C→T R271H (CGC→CAC)  cirA ← colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,245,9580CT98.2% 564.9 / NA 170R271H (CGC→CAC) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (74/92);  minor base G (1/2);  total (76/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGA  >  NC_000913/2245861‑2246054
                                                                                                 |                                                                                                
gTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTgg                                                                                                <  2:4436196/100‑1 (MQ=255)
 tttttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAggtggt                                                                                               <  2:4988756/100‑1 (MQ=255)
  ttttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGtt                                                                                              >  1:2532311/1‑100 (MQ=255)
   tttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGctt                                                                                             >  1:4336781/1‑100 (MQ=255)
    ttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGtttt                                                                                            <  1:1899477/100‑1 (MQ=255)
      gTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTg                                                                                          >  2:4114341/1‑100 (MQ=255)
       ttCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGt                                                                                         >  2:2002975/1‑100 (MQ=255)
       ttCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGt                                                                                         <  1:3013604/100‑1 (MQ=255)
        tctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTc                                                                                        >  2:3243178/1‑100 (MQ=255)
         ctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTcc                                                                                       <  1:3962159/100‑1 (MQ=255)
          tcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCa                                                                                      >  2:1754161/1‑100 (MQ=255)
          tcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCa                                                                                      >  1:801188/1‑100 (MQ=255)
           cGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAg                                                                                     >  2:3725549/1‑100 (MQ=255)
            gACTTTCTGACCGTGGTATTTCAGTTCGCTGGTGCCGTAATGCCAACGCCAATTTTGGACGACGGAGTAGTTGTTGCGTGCCAGGGGGGTTTTGTCCAGc                                                                                    >  2:3518501/1‑100 (MQ=255)
            gACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCATCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGc                                                                                    <  2:4758276/100‑1 (MQ=255)
            gACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGc                                                                                    <  1:3130662/100‑1 (MQ=255)
            gACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGc                                                                                    <  1:3442514/100‑1 (MQ=255)
            gACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGc                                                                                    >  1:2456500/1‑100 (MQ=255)
             aCTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCg                                                                                   <  1:521651/100‑1 (MQ=255)
             aCTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCg                                                                                   <  2:2893479/100‑1 (MQ=255)
                  ctcACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCg                                                                              <  2:583038/100‑1 (MQ=255)
                  ctcACCGTAGTATTTCAGTTCGCTGATGCCGTCATCACAACGCCCATTATGGCTCACGCAGTAGTTCTGACGTTCCAAGTGGTCTTTCTCCAACGAGTCg                                                                              >  2:3330085/1‑100 (MQ=255)
                   tcACCGTAGTATTTCAGTGCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCgg                                                                             <  1:1915407/100‑1 (MQ=255)
                      ccGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAAt                                                                          >  1:4656026/1‑100 (MQ=255)
                      ccGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAAt                                                                          >  2:2422878/1‑100 (MQ=255)
                      ccGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAAt                                                                          >  2:23814/1‑100 (MQ=255)
                      ccGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAAt                                                                          <  1:2241676/100‑1 (MQ=255)
                        gtagtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCa                                                                        >  2:4015522/1‑100 (MQ=255)
                          agtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACg                                                                      <  1:1080131/100‑1 (MQ=255)
                          agtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACg                                                                      >  2:3526139/1‑100 (MQ=255)
                             aTTTCAGTTCGCTGGTGTCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATc                                                                   >  1:1594779/1‑100 (MQ=255)
                               ttCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCt                                                                 >  1:850378/1‑100 (MQ=255)
                               ttCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCt                                                                 <  1:3987075/100‑1 (MQ=255)
                                tCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTg                                                                >  1:1918924/1‑100 (MQ=255)
                                  aGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGAc                                                              <  1:3565216/100‑1 (MQ=255)
                                   gTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGCGTCGGAATCACGATCCTGACg                                                             >  2:3881383/1‑100 (MQ=255)
                                   gTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACg                                                             >  2:2162976/1‑100 (MQ=255)
                                    ttCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACgg                                                            <  1:3199570/100‑1 (MQ=255)
                                     tCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGt                                                           <  2:1828601/100‑1 (MQ=255)
                                       gCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCg                                                         <  1:4437487/100‑1 (MQ=255)
                                        cTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGa                                                        <  1:195454/100‑1 (MQ=255)
                                         tGGTACCGTAATCCCAACGCCCATTATGGCTGACGGAGGAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGaa                                                       >  2:2694332/1‑100 (MQ=255)
                                          ggTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGaaa                                                      <  2:2725613/100‑1 (MQ=255)
                                           gTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTGCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAAc                                                     <  1:2473632/100‑1 (MQ=255)
                                             gCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCg                                                   <  2:3080939/100‑1 (MQ=255)
                                             gCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCg                                                   <  2:3468167/100‑1 (MQ=255)
                                              ccGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGt                                                  >  1:1461525/1‑100 (MQ=255)
                                              ccGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGt                                                  >  2:3642151/1‑100 (MQ=255)
                                               cGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTa                                                 >  2:3238610/1‑100 (MQ=255)
                                                gTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTAt                                                <  2:2889974/100‑1 (MQ=255)
                                                gTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTAt                                                >  1:1197040/1‑100 (MQ=255)
                                                 tAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATc                                               >  2:1381455/1‑100 (MQ=255)
                                                  aaTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGCTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATcc                                              <  2:715939/100‑1 (MQ=255)
                                                     cccAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                           <  1:1633081/100‑1 (MQ=255)
                                                     cccAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                           >  2:2603648/1‑100 (MQ=255)
                                                     cccAACGCCCATTATGACTGACGGAGGAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                           >  1:4519277/1‑100 (MQ=255)
                                                        aaCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGt                                        >  2:1504974/1‑100 (MQ=255)
                                                        aaCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGt                                        >  2:4670755/1‑100 (MQ=255)
                                                        aaCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGt                                        <  2:4380987/100‑1 (MQ=255)
                                                           gCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTaaa                                     <  2:1737591/100‑1 (MQ=255)
                                                            cccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTaaaa                                    >  1:3914113/1‑100 (MQ=255)
                                                              cATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATc                                  >  2:2320155/1‑100 (MQ=255)
                                                              cATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATc                                  <  1:2177587/100‑1 (MQ=255)
                                                               aTTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCg                                 <  2:4644146/100‑1 (MQ=255)
                                                                 tATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTg                               <  2:1431514/100‑1 (MQ=255)
                                                                 tATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTg                               <  1:2028154/100‑1 (MQ=255)
                                                                    ggCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtt                            >  1:626/1‑100 (MQ=255)
                                                                    ggCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtt                            >  1:2771264/1‑100 (MQ=255)
                                                                    ggCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtt                            <  1:1760496/100‑1 (MQ=255)
                                                                     gCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAttt                           >  1:2942903/1‑100 (MQ=255)
                                                                     gCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAttt                           <  1:3279606/100‑1 (MQ=255)
                                                                     gCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAttt                           >  1:3905209/1‑100 (MQ=255)
                                                                      cTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                          <  1:4516489/100‑1 (MQ=255)
                                                                        gACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGa                        >  1:2892272/1‑100 (MQ=255)
                                                                        gACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGa                        <  2:3418032/100‑1 (MQ=255)
                                                                         aCGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGATACCGTATCCGGCAGTAAAATCGTGATTTTGAt                       <  2:4577694/100‑1 (MQ=255)
                                                                         aCGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGAt                       <  2:2384206/100‑1 (MQ=255)
                                                                            gAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCg                    <  2:4463974/100‑1 (MQ=255)
                                                                             agtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCgg                   <  2:837015/100‑1 (MQ=255)
                                                                              gtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGTCAGTAAAATCGTGATTTTGATTCgct                  <  2:2151080/100‑3 (MQ=255)
                                                                              gtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGgt                  <  2:2859952/100‑1 (MQ=255)
                                                                              gtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAAACGTGATTTTGATTCGgt                  >  1:511067/1‑100 (MQ=255)
                                                                                 gtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc               >  1:3164903/1‑100 (MQ=255)
                                                                                 gtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGACGGAATCACGATCCTGACTGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc               >  1:4781913/1‑100 (MQ=255)
                                                                                  ttCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTcc              <  1:5147009/100‑1 (MQ=255)
                                                                                     tGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAgg           <  2:246211/100‑1 (MQ=255)
                                                                                     tGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAgg           <  1:2633789/100‑1 (MQ=255)
                                                                                     tGGCGTTCCAGGGGGTGTTTGTCCAGCAAGCCGGAATTACGACCCGGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAgg           <  2:3568332/100‑1 (MQ=255)
                                                                                      ggCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAGTTTGATTCGGTGTCCAGGc          >  2:3629104/1‑100 (MQ=255)
                                                                                       gCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCa         <  2:3943054/100‑1 (MQ=255)
                                                                                        cGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCaa        >  1:3403022/1‑100 (MQ=255)
                                                                                            ccAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTc    <  2:4209089/100‑1 (MQ=255)
                                                                                            ccAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTc    >  2:2404254/1‑100 (MQ=255)
                                                                                            ccAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAAGTc    >  2:1608125/1‑100 (MQ=255)
                                                                                             cAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCg   >  2:2736333/1‑100 (MQ=255)
                                                                                             cAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCg   <  2:58357/100‑1 (MQ=255)
                                                                                              aGGTGGTTTTTGTCTAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTACAATCCTGCTTTTGTTTCGGTGTCCAGGCAGATTCGa  >  2:459140/1‑100 (MQ=255)
                                                                                              aGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGa  >  2:3287964/1‑100 (MQ=255)
                                                                                              aGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGa  >  1:577853/1‑100 (MQ=255)
                                                                                              aGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGa  >  1:417050/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGA  >  NC_000913/2245861‑2246054

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: