Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 999,571 A→G Q378R (CAG→CGG)  ycbS → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913999,5710AG100.0% 487.7 / NA 153Q378R (CAG→CGG) ycbSputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (69/84);  total (69/84)

Only 100 of 159 total aligned reads displayed.

CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGT  >  NC_000913/999481‑999662
                                                                                          |                                                                                           
cGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCggg                                                                                           <  1:639601/93‑1 (MQ=255)
 gTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGga                                                                                          >  2:2551008/1‑93 (MQ=255)
  ttCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgag                                                                                         >  1:550928/1‑93 (MQ=255)
   tctcCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgaga                                                                                        <  2:1580904/93‑1 (MQ=255)
   tctcCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgaga                                                                                        <  2:4248947/93‑1 (MQ=255)
     tcCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGc                                                                                      <  2:2145332/93‑1 (MQ=255)
        aGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaaa                                                                                   <  2:1104556/93‑1 (MQ=255)
         gtgtGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAt                                                                                  <  1:2401162/93‑1 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaaact                                                                                 >  1:1714053/1‑91 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAtt                                                                                 >  1:4383661/1‑93 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAtt                                                                                 >  2:491328/1‑93 (MQ=255)
            tgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTg                                                                               >  1:1007925/1‑93 (MQ=255)
             gCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGc                                                                              <  1:3851571/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                           <  2:359552/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                           <  2:2913528/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                           >  1:505146/1‑93 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCCAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                           >  2:2693573/1‑93 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTAGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                           >  1:2506947/1‑93 (MQ=255)
                 ttACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAg                                                                          >  1:4631657/1‑93 (MQ=255)
                 ttACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAg                                                                          <  2:3015306/93‑1 (MQ=255)
                 ttACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAg                                                                          <  2:2674102/93‑1 (MQ=255)
                         gCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGtt                                                                  <  2:2409758/93‑1 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                                >  1:444908/1‑93 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                                <  2:1886863/93‑1 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                                >  2:2509603/1‑93 (MQ=255)
                            tCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCa                                                               <  1:1522245/93‑1 (MQ=255)
                              aGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGt                                                             <  1:2357596/93‑1 (MQ=255)
                               ggggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTg                                                            >  1:1122064/1‑93 (MQ=255)
                               ggggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTg                                                            >  2:3328247/1‑93 (MQ=255)
                               ggggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTg                                                            >  1:4243676/1‑93 (MQ=255)
                                gggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgg                                                           >  2:432562/1‑93 (MQ=255)
                                 ggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTggg                                                          >  1:320794/1‑93 (MQ=255)
                                 ggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTggg                                                          <  1:3709183/93‑1 (MQ=255)
                                 ggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTggg                                                          >  2:2585000/1‑93 (MQ=255)
                                  gCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgggg                                                         <  2:2239726/93‑1 (MQ=255)
                                    gAATCAAATACGCGGTGACACTGGCGAGATCCAGAACCAATAGTAATGAACAGCGGGCAAGCAAATTTGCCCAGGCCACGTTGCTGCGGGGCg                                                       >  2:526363/1‑93 (MQ=255)
                                       tCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGAc                                                    >  2:731074/1‑93 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCgg                                                 >  2:492512/1‑92 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGt                                                 <  1:2707945/93‑1 (MQ=255)
                                           aTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                                >  1:3896448/1‑93 (MQ=255)
                                           aTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                                <  2:3889715/93‑1 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                               >  1:1346686/1‑93 (MQ=255)
                                              cgcgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgggg                                             <  1:501870/93‑1 (MQ=255)
                                                cgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGaa                                           <  1:2230966/93‑1 (MQ=255)
                                                 ggTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAAc                                          >  2:1200952/1‑93 (MQ=255)
                                                 ggTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAAc                                          <  2:4667075/93‑1 (MQ=255)
                                                  gTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACg                                         >  1:2944383/1‑93 (MQ=255)
                                                  gTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACg                                         >  2:3523340/1‑93 (MQ=255)
                                                   tGACACTGGCGAAATACAGAACGAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGa                                        <  2:1385021/93‑1 (MQ=255)
                                                   tGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACgc                                        >  1:2011770/1‑92 (MQ=255)
                                                     acacTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGAcc                                      <  2:3044246/93‑1 (MQ=255)
                                                      cacTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCt                                     <  1:2622212/93‑1 (MQ=255)
                                                       acTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTg                                    >  1:4224802/1‑93 (MQ=255)
                                                        cTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTgg                                   >  1:286315/1‑93 (MQ=255)
                                                         tGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGt                                  >  1:2940143/1‑93 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                                 <  2:4411060/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                                 <  2:3354520/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                                <  1:1165876/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                                <  1:1413594/93‑1 (MQ=255)
                                                            cGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtg                               <  2:1348556/93‑1 (MQ=255)
                                                            cGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtg                               <  2:2145118/93‑1 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                              >  1:3354082/1‑93 (MQ=255)
                                                             gAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtgg                              <  1:1934427/93‑1 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggt                             <  2:3660652/93‑1 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggt                             <  2:968129/93‑1 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggt                             <  1:2885254/93‑1 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggg                             >  2:3227924/1‑92 (MQ=255)
                                                               aaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGTCCGTGGGGAACGACCTGGTAtggtg                            >  1:4607912/1‑93 (MQ=255)
                                                               aaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtg                            >  2:785785/1‑93 (MQ=255)
                                                               aaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtg                            <  1:759762/93‑1 (MQ=255)
                                                                 tACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggt                          <  1:2159167/93‑1 (MQ=255)
                                                                 tACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggt                          <  1:2610826/93‑1 (MQ=255)
                                                                 tACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggg                          >  2:71369/1‑92 (MQ=255)
                                                                 tACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAAAGACCTGGTAtggtggt                          >  1:3207638/1‑93 (MQ=255)
                                                                  aCAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtg                         <  2:2339985/93‑1 (MQ=255)
                                                                  aCAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtg                         >  1:1619698/1‑93 (MQ=255)
                                                                   cAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtgg                        <  1:4531071/93‑1 (MQ=255)
                                                                        ccAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAat                   >  2:3967934/1‑93 (MQ=255)
                                                                        ccAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAat                   <  2:3703062/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                  <  1:1552063/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                  <  2:779081/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata                  <  1:337231/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGCGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGGATGGTGGTGGACAatat                 >  2:3818868/1‑93 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat                 >  1:1381251/1‑93 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat                 >  1:2621753/1‑93 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat                 <  2:153957/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat                 <  2:518654/93‑1 (MQ=255)
                                                                           aTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATg                >  1:2962586/1‑93 (MQ=255)
                                                                               tAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGa            >  1:945656/1‑93 (MQ=255)
                                                                               tAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGa            <  1:1906480/93‑1 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGaa           <  2:282891/93‑1 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGaa           >  1:1628187/1‑93 (MQ=255)
                                                                                aaTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGaa           <  1:633831/93‑1 (MQ=255)
                                                                                 aTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAat          >  2:1884387/1‑93 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata         <  2:4529283/93‑1 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata         <  1:1207662/93‑1 (MQ=255)
                                                                                   gAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAatat        <  2:1036306/93‑1 (MQ=255)
                                                                                     aCAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTa      <  2:615902/93‑1 (MQ=255)
                                                                                       agcggGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTAcc    >  2:2856098/1‑93 (MQ=255)
                                                                                         cagGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGt  <  1:752695/93‑1 (MQ=255)
                                                                                          |                                                                                           
CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTACCGT  >  NC_000913/999481‑999662

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