Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,018,363 G→A P305L (CCG→CTG)  ompA ← outer membrane protein 3a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,3630GA85.5% 710.4 / 85.3 249P305L (CCG→CTG) ompAouter membrane protein 3a
Reads supporting (aligned to +/- strand):  ref base G (17/19);  new base A (122/91);  total (139/110)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 255 total aligned reads displayed.

CTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCG  >  NC_000913/1018271‑1018455
                                                                                            |                                                                                            
  ttAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCagg                                                                                            <  1:6229252/93‑1 (MQ=255)
  ttAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCagg                                                                                            <  1:1844879/93‑1 (MQ=255)
   tAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGt                                                                                           <  2:3750675/93‑1 (MQ=255)
    aaCTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGtt                                                                                          >  2:2396893/1‑93 (MQ=255)
     aCTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTg                                                                                         >  2:3330897/1‑93 (MQ=255)
       ttCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGa                                                                                       <  1:6536493/93‑1 (MQ=255)
         ggATCTCTACGCGACAATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGAtt                                                                                     <  2:4637052/92‑1 (MQ=255)
         cGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGAtt                                                                                     <  1:3973036/93‑1 (MQ=255)
           aTCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCg                                                                                   >  2:3057144/1‑93 (MQ=255)
           aTCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCg                                                                                   >  1:3029163/1‑93 (MQ=255)
            tctcTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGCGTTGCCAGTAACCAGGTTGGATTCGc                                                                                  >  2:2482602/1‑93 (MQ=255)
              tctACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGccc                                                                                >  2:2329135/1‑93 (MQ=255)
               ctACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCa                                                                               <  1:1935141/93‑1 (MQ=255)
                tACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCAt                                                                              <  2:4718283/93‑1 (MQ=255)
                 aCGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATa                                                                             <  2:5489431/93‑1 (MQ=255)
                    cgacgaTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCa                                                                          >  2:731114/1‑93 (MQ=255)
                     gacgaTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCAc                                                                         >  1:4439938/1‑93 (MQ=255)
                     gacgaTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCAc                                                                         <  1:3673809/93‑1 (MQ=255)
                      acgaTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACg                                                                        <  2:2199526/93‑1 (MQ=255)
                       cgaTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGt                                                                       >  2:1681082/1‑93 (MQ=255)
                         aTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGc                                                                     >  2:4541071/1‑93 (MQ=255)
                         aTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGc                                                                     >  1:4107706/1‑93 (MQ=255)
                         aTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGc                                                                     <  1:1534863/93‑1 (MQ=255)
                             ggAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGgag                                                                 <  2:2735389/93‑1 (MQ=255)
                             ggAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGgag                                                                 <  1:5711613/93‑1 (MQ=255)
                             ggAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGgag                                                                 >  1:1495203/1‑93 (MQ=255)
                             ggAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGCTGGATTCGCCCATACCACGTGCGgag                                                                 >  1:1736133/1‑93 (MQ=255)
                              gAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGgaga                                                                >  2:3417264/1‑93 (MQ=255)
                              gAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGgaga                                                                <  1:5440040/93‑1 (MQ=255)
                               aGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGAt                                                               >  1:716570/1‑93 (MQ=255)
                                   aGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTg                                                           >  1:2234311/1‑93 (MQ=255)
                                      cAGTCGATCAGTGCATCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCt                                                        >  2:1010435/1‑93 (MQ=255)
                                       aGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTg                                                       >  2:5220171/1‑93 (MQ=255)
                                       aGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTg                                                       >  2:3035042/1‑93 (MQ=255)
                                       aGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTg                                                       >  1:2090524/1‑93 (MQ=255)
                                        gTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGc                                                      <  1:285275/93‑1 (MQ=255)
                                         tCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGcc                                                     >  2:3242497/1‑93 (MQ=255)
                                          cGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCACACCACGTGCGGAGATCTTGTCTGCCg                                                    >  1:5356205/1‑93 (MQ=255)
                                           gATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGcccg                                                   >  1:3818312/1‑91 (MQ=255)
                                           gATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCgg                                                   <  1:5219288/93‑1 (MQ=255)
                                           gATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCgg                                                   >  1:3387366/1‑93 (MQ=255)
                                            aTCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCggg                                                  >  1:708003/1‑93 (MQ=255)
                                            aTCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCggg                                                  >  1:5015766/1‑93 (MQ=255)
                                            aTCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCggg                                                  >  1:1677756/1‑93 (MQ=255)
                                             tCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGa                                                 >  2:3521918/1‑93 (MQ=255)
                                             tCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGa                                                 >  2:1207405/1‑93 (MQ=255)
                                             tCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGa                                                 <  1:3291672/93‑1 (MQ=255)
                                              cAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGAt                                                >  2:6576295/1‑93 (MQ=255)
                                              cAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGAt                                                <  2:1045591/93‑1 (MQ=255)
                                                  gcagcaCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCt                                            <  2:4132735/93‑1 (MQ=255)
                                                   cagcaCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCtt                                           <  2:1321021/93‑1 (MQ=255)
                                                   cagcaCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCtt                                           <  1:4982486/93‑1 (MQ=255)
                                                    agcaCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGGTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCttt                                          >  2:5810637/1‑93 (MQ=255)
                                                      caCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTgg                                        >  1:3238609/1‑93 (MQ=255)
                                                       aCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGga                                       >  2:5295110/1‑93 (MQ=255)
                                                       aCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGga                                       <  2:3958947/93‑1 (MQ=255)
                                                       aCGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGga                                       <  1:6578946/93‑1 (MQ=255)
                                                          cTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGAGATACCTTTGGAGAt                                    >  2:3083627/1‑93 (MQ=255)
                                                            gTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCa                                  >  1:5120112/1‑93 (MQ=255)
                                                            gTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCa                                  <  2:1747862/93‑1 (MQ=255)
                                                            gTTTCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCa                                  >  1:3348048/1‑93 (MQ=255)
                                                               tCACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGt                               >  2:1874112/1‑93 (MQ=255)
                                                                cACGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTa                              >  1:451103/1‑93 (MQ=255)
                                                                 ccGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTaa                             <  2:6497003/92‑1 (MQ=255)
                                                                 aCGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTaa                             >  2:6532059/1‑93 (MQ=255)
                                                                 aCGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTaa                             >  1:5962776/1‑93 (MQ=255)
                                                                 aCGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTaa                             >  1:4635210/1‑93 (MQ=255)
                                                                  cGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAAt                            <  2:5063270/93‑1 (MQ=255)
                                                                  cGTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAAt                            >  2:6056613/1‑93 (MQ=255)
                                                                   gTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTACTc                           >  1:2724807/1‑93 (MQ=255)
                                                                   gTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATc                           >  1:5580893/1‑93 (MQ=255)
                                                                   gTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATc                           <  2:2684255/93‑1 (MQ=255)
                                                                   gTTGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATc                           <  1:3551737/93‑1 (MQ=255)
                                                                    ttGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATca                          >  2:1528036/1‑93 (MQ=255)
                                                                    ttGTCAAAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATca                          >  1:368756/1‑93 (MQ=255)
                                                                     tGTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaa                         <  2:1956093/93‑1 (MQ=255)
                                                                      gTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaac                        <  1:231013/93‑1 (MQ=255)
                                                                      gTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaac                        >  2:6368646/1‑93 (MQ=255)
                                                                      gTCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaac                        >  1:479001/1‑93 (MQ=255)
                                                                       tCACAGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaaca                       <  2:3758791/93‑1 (MQ=255)
                                                                        cacaGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaacaa                      <  1:6083201/93‑1 (MQ=255)
                                                                        cacaGGCGTTCCCAGTAACCAGGTTAGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATcaacaa                      <  2:6187849/93‑1 (MQ=255)
                                                                          caGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAacaaca                    >  1:1050363/1‑93 (MQ=255)
                                                                          caGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAacaaca                    <  2:4314239/93‑1 (MQ=255)
                                                                           aGGTGTTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAg                   >  1:2912862/1‑93 (MQ=255)
                                                                            ggTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGa                  >  2:6611787/1‑93 (MQ=255)
                                                                             gtgtTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGAc                 <  2:6009143/93‑1 (MQ=255)
                                                                             gtgtTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGAc                 >  2:4794558/1‑93 (MQ=255)
                                                                               gtTGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTg               >  1:5556459/1‑93 (MQ=255)
                                                                                ttGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGa              <  1:5909848/93‑1 (MQ=255)
                                                                                ttGCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGa              >  1:1307325/1‑93 (MQ=255)
                                                                                  gCCAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGc            >  1:786002/1‑93 (MQ=255)
                                                                                    cAGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCAc          >  2:5629212/1‑93 (MQ=255)
                                                                                     aGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACg         >  2:2569022/1‑93 (MQ=255)
                                                                                     aGTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACg         >  2:2037419/1‑93 (MQ=255)
                                                                                      gTAACCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg        <  2:5288181/93‑1 (MQ=255)
                                                                                       tAACCAGGTTGGATTCGCCAATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTCGGAGATCAGGTAATCAACAACAGACTGAGCACGgc       <  2:4290009/93‑1 (MQ=255)
                                                                                         aCCAGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGgcgc     >  2:4213606/1‑93 (MQ=255)
                                                                                          ccGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCt    <  1:3810283/93‑1 (MQ=255)
                                                                                            gggTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCg  >  1:4355701/1‑93 (MQ=255)
                                                                                            |                                                                                            
CTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCG  >  NC_000913/1018271‑1018455

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