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,9580CT100.0% 295.4 / NA 88R271H (CGC→CAC) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/56);  total (32/56)

TTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCG  >  NC_000913/2245864‑2246053
                                                                                              |                                                                                               
tttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGttt                                                                                            <  1:1081229/100‑1 (MQ=255)
   gTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCATCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTg                                                                                         <  2:107904/100‑1 (MQ=255)
   gTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTg                                                                                         <  1:1449351/100‑1 (MQ=255)
   gTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTg                                                                                         <  1:1308423/100‑1 (MQ=255)
     tctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTc                                                                                       >  1:1554376/1‑100 (MQ=255)
     tctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTc                                                                                       <  2:1610412/100‑1 (MQ=255)
     tctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTc                                                                                       <  1:1209328/100‑1 (MQ=255)
     tctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTc                                                                                       >  1:1246737/1‑100 (MQ=255)
      ctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTcc                                                                                      >  2:1183293/1‑100 (MQ=255)
      ctcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTcc                                                                                      <  1:1720743/100‑1 (MQ=255)
       tcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCa                                                                                     >  1:1565467/1‑100 (MQ=255)
       tcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTGTTTGTCCa                                                                                     <  2:2484851/100‑1 (MQ=255)
          aCTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCg                                                                                  <  1:1470472/100‑1 (MQ=255)
             ttCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGt                                                                               >  2:532504/1‑100 (MQ=255)
              tctcACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTc                                                                              >  2:1843685/1‑100 (MQ=255)
               ctcACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCg                                                                             <  1:618794/100‑1 (MQ=255)
                 cACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGa                                                                           >  2:139266/1‑100 (MQ=255)
                 cACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGa                                                                           >  1:1270022/1‑100 (MQ=255)
                    cGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATc                                                                        <  2:2570844/100‑1 (MQ=255)
                    cGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATc                                                                        <  1:1397162/100‑1 (MQ=255)
                      tagtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCAc                                                                      <  1:1152948/100‑1 (MQ=255)
                       agtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACg                                                                     <  2:1771207/100‑1 (MQ=255)
                       aataTTTCAGTGCGCAGGTGCCGTAATCCCAACTCCCACTATGGCTGACTTGGTAATTCTGGCGGTCCAGGTAGTTTTTGTGCCGCGAGTCGGAATCACg                                                                     <  2:201121/98‑1 (MQ=255)
                          aTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATc                                                                  <  2:679983/100‑1 (MQ=255)
                          aTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATc                                                                  >  1:1310186/1‑100 (MQ=255)
                            ttCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTTTGGCTGACGGAGGAGTTCTGGCGTTCCCGGTGGTTTTTGTCCCGCGGGTCGGAATCACGATCCt                                                                >  2:1619646/1‑100 (MQ=255)
                               aGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGCTCCTGAc                                                             <  2:1007771/100‑1 (MQ=255)
                               aGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGAc                                                             >  2:1471915/1‑100 (MQ=255)
                                gTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACg                                                            >  1:62934/1‑100 (MQ=255)
                                 ttCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCTTGCGTTTCGGAATCACGATCCTGACgg                                                           >  2:1965409/1‑100 (MQ=255)
                                 ttCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACgg                                                           >  1:672982/1‑100 (MQ=255)
                                 ttCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACgg                                                           <  1:2620588/100‑1 (MQ=255)
                                   cGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTc                                                         <  1:2056404/100‑1 (MQ=255)
                                    gCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCg                                                        <  1:2544818/100‑1 (MQ=255)
                                     cTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGa                                                       <  1:29225/100‑1 (MQ=255)
                                     cTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGa                                                       <  2:854693/100‑1 (MQ=255)
                                        gTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAAc                                                    <  1:868130/100‑1 (MQ=255)
                                         tGCCGTAATCCCACCGCCCATTACAGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCAAGTCGGAATCACGAGCCTGGCGGTCGAAAcc                                                   <  2:2040347/100‑1 (MQ=255)
                                         tGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAAcc                                                   <  1:2057429/100‑1 (MQ=255)
                                         tGCCGTAATCCCAACGCCCATTATGACTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAAcc                                                   <  2:2228065/100‑1 (MQ=255)
                                            cGTAATCCCAACGCCCATTATGGCTGATGGCGTAGTTCTGGCGTTCCCGGTGGTTTTTGTCCAGTGTGTCGGAACCCGGATCCTGCGGGTCGACATCGTa                                                >  2:1126386/1‑100 (MQ=255)
                                             gTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTAt                                               >  2:702584/1‑100 (MQ=255)
                                             gTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTAt                                               >  1:125369/1‑100 (MQ=255)
                                                aTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCg                                            >  2:1843978/1‑100 (MQ=255)
                                                aTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCg                                            >  1:1126019/1‑100 (MQ=255)
                                                  cccAGCGCCCATGATGGCGGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                          <  2:2259005/100‑1 (MQ=255)
                                                  cccAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                          <  1:2237788/100‑1 (MQ=255)
                                                  cccAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                          <  2:1420365/100‑1 (MQ=255)
                                                  cccAACGCCCATTATGGCTGACGGAGTAGTTCGGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGc                                          <  2:997374/100‑1 (MQ=255)
                                                    cAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAg                                        <  1:2189907/100‑1 (MQ=255)
                                                     accgtccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTAGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGt                                       <  1:2372547/95‑1 (MQ=255)
                                                     aaCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGt                                       >  2:2107320/1‑100 (MQ=255)
                                                         cccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTaaaa                                   >  1:711690/1‑100 (MQ=255)
                                                          ccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAAt                                  <  2:863163/100‑1 (MQ=255)
                                                             ttATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGt                               <  1:816378/100‑1 (MQ=255)
                                                              tATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTg                              <  1:672991/100‑1 (MQ=255)
                                                              tATGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTg                              <  1:1308050/100‑1 (MQ=255)
                                                                tGGCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAt                            <  1:1072505/100‑1 (MQ=255)
                                                                 ggCTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtt                           >  1:2066389/1‑100 (MQ=255)
                                                                   cTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                         <  2:2114467/100‑1 (MQ=255)
                                                                   cTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                         <  1:2486473/100‑1 (MQ=255)
                                                                   cTGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                         <  1:1260553/100‑1 (MQ=255)
                                                                    tGACGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTg                        <  1:1525618/100‑1 (MQ=255)
                                                                      aCGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGAt                      <  1:467643/100‑1 (MQ=255)
                                                                       cGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGAtt                     >  1:37277/1‑100 (MQ=255)
                                                                       cGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGAtn                     <  2:1126097/100‑2 (MQ=255)
                                                                       cGGAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTAATTTTGAtt                     >  2:1518954/1‑100 (MQ=255)
                                                                        ggAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTc                    <  2:798749/100‑1 (MQ=255)
                                                                         gAGTAGTTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCg                   >  1:704215/1‑100 (MQ=255)
                                                                          agtggtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAAGCACGATCCTGACGGTCGAAACTGGATCCGGCAGTAAAATCGTGATTTTGATTCgg                  <  2:1635146/100‑1 (MQ=255)
                                                                          agtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCgg                  <  1:1375458/100‑1 (MQ=255)
                                                                          agtagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCgg                  >  1:825855/1‑100 (MQ=255)
                                                                            tagtTCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGgtg                <  2:671616/100‑1 (MQ=255)
                                                                              gtTCTGGCGTTGCAGGTGGTCTCTGTCCAGCGAGTCGGAATAACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc              <  2:1228481/100‑1 (MQ=255)
                                                                                tCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCa            <  1:2088645/100‑1 (MQ=255)
                                                                                tCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCa            <  1:2567263/100‑1 (MQ=255)
                                                                                tCTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCa            >  1:1190481/1‑100 (MQ=255)
                                                                                 cTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCCGTAAAATCGTGATTTTGATTCGGTGTCCAg           >  1:1580418/1‑100 (MQ=255)
                                                                                 cTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAg           <  2:2228286/100‑1 (MQ=255)
                                                                                 cTGGCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAg           <  1:985105/100‑1 (MQ=255)
                                                                                   ggCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAACTCGTGCTTTTCATTCTGTGTcaacgc         >  2:2055025/1‑95 (MQ=255)
                                                                                   ggCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGc         <  2:245633/100‑1 (MQ=255)
                                                                                   ggCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGc         <  1:266364/100‑1 (MQ=255)
                                                                                   ggCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGc         >  2:824472/1‑100 (MQ=255)
                                                                                   ggCGTTCCAGGTAGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGc         >  1:2271379/1‑100 (MQ=255)
                                                                                    gCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCa        <  2:71767/100‑1 (MQ=255)
                                                                                    gCGTTCCAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCa        >  1:348962/1‑100 (MQ=255)
                                                                                          cAGGTGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCg  <  1:1661497/100‑1 (MQ=255)
                                                                                              |                                                                                               
TTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCG  >  NC_000913/2245864‑2246053

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: