Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,902,432 G→A P265L (CCG→CTG)  bglF ← PTS system beta‑glucosides‑specific enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,902,4320GA99.5% 1363.1 / ‑5.2 378P265L (CCG→CTG) bglFPTS system beta‑glucosides‑specific enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (189/186);  total (190/188)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 388 total aligned reads displayed.

CCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902342‑3902521
                                                                                          |                                                                                         
  gccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGccc                                                                                       <  1:3705408/93‑1 (MQ=255)
   cccATTACCGTGCCCGCTAATGCAGGAACCGCCTGATAAAGCCAGTGATTACCGGCGGCAATCAGTTCGCTTATCTAGGTTGATAGCATcccc                                                                                      >  2:7035887/1‑93 (MQ=255)
   cccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcccc                                                                                      >  1:1601370/1‑93 (MQ=255)
   cccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcccc                                                                                      >  2:5551854/1‑93 (MQ=255)
   cccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcccc                                                                                      >  2:2336084/1‑93 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    >  1:570816/1‑93 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    <  1:3888583/93‑1 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAAGTTGATAGCAGCCccac                                                                                    >  2:6425839/1‑93 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCCGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    >  1:6713515/1‑93 (MQ=255)
       ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  2:2972357/1‑93 (MQ=255)
       ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  1:2775919/1‑93 (MQ=255)
        tACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcag                                                                                 <  1:5446388/93‑1 (MQ=255)
        tACCGCGCACGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcag                                                                                 <  1:2716931/93‑1 (MQ=255)
         aCCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcagc                                                                                <  2:4633659/93‑1 (MQ=255)
         aCCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcagc                                                                                <  1:1654442/93‑1 (MQ=255)
          ccGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcagca                                                                               <  1:2462039/93‑1 (MQ=255)
            gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             <  2:3105303/93‑1 (MQ=255)
            gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             <  2:7428401/93‑1 (MQ=255)
             cgcCCGCAAATGCAGAAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaa                                                                            <  2:3599768/93‑1 (MQ=255)
              gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaaa                                                                           <  2:7254341/93‑1 (MQ=255)
               cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAg                                                                          <  2:4082744/93‑1 (MQ=255)
                ccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAgg                                                                         >  1:5178546/1‑93 (MQ=255)
                  gCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTg                                                                       >  1:3539372/1‑93 (MQ=255)
                   cAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGa                                                                      <  1:2953145/93‑1 (MQ=255)
                    aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGAc                                                                     <  1:6582849/93‑1 (MQ=255)
                      aTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACgg                                                                   <  2:1137752/93‑1 (MQ=255)
                       tGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACggg                                                                  <  1:7151525/93‑1 (MQ=255)
                        gCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACgggg                                                                 >  2:5172914/1‑92 (MQ=255)
                        gCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACgggg                                                                 >  1:604499/1‑92 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  1:7540250/93‑1 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  1:5355494/93‑1 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  1:3946207/93‑1 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  1:7866041/93‑1 (MQ=255)
                                cGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTGCGCTTATCCAGTTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAc                                                         <  2:5059466/93‑1 (MQ=255)
                                 gCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAcc                                                        >  1:907585/1‑93 (MQ=255)
                                  ccTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCa                                                       <  1:4563923/93‑1 (MQ=255)
                                   cTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCAt                                                      >  2:3319568/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  2:2698856/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:4681587/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:4660561/1‑93 (MQ=255)
                                     gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCa                                                    >  1:2355434/1‑93 (MQ=255)
                                     gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCa                                                    <  2:6687932/93‑1 (MQ=255)
                                     gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCa                                                    <  2:1619378/93‑1 (MQ=255)
                                       tAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGa                                                  >  2:374266/1‑93 (MQ=255)
                                       tAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGa                                                  >  1:6905909/1‑93 (MQ=255)
                                       tAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGa                                                  <  1:1897166/93‑1 (MQ=255)
                                        aaaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGAc                                                 >  2:35054/1‑93 (MQ=255)
                                         aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACa                                                <  1:49540/93‑1 (MQ=255)
                                         aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACa                                                <  2:4338651/93‑1 (MQ=255)
                                         aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACa                                                <  1:3298787/93‑1 (MQ=255)
                                           gCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATa                                              <  1:7387842/93‑1 (MQ=255)
                                           gCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATa                                              <  1:1165943/93‑1 (MQ=255)
                                            caagagATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATATCCATCAGACATAg                                             <  2:208303/91‑1 (MQ=255)
                                             cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGc                                            <  2:2900451/93‑1 (MQ=255)
                                               gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          <  2:3954538/93‑1 (MQ=255)
                                               gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          >  2:850966/1‑93 (MQ=255)
                                               gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          <  2:1473200/93‑1 (MQ=255)
                                                agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAAt                                         >  1:2223205/1‑93 (MQ=255)
                                                 gaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATg                                        >  2:6058441/1‑93 (MQ=255)
                                                  aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATgg                                       <  1:3434339/93‑1 (MQ=255)
                                                   tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      >  2:7724938/1‑93 (MQ=255)
                                                    aaCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                     >  1:6439649/1‑93 (MQ=255)
                                                     aCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtgt                                    <  2:1903641/93‑1 (MQ=255)
                                                      ccGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGtgtg                                   <  1:624117/93‑1 (MQ=255)
                                                       cggcggCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                  <  1:4826075/93‑1 (MQ=255)
                                                       cggcggCAATCAGTTCGCTTATCCAGGTGGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                  <  1:375754/93‑1 (MQ=255)
                                                        ggcggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                 <  1:566561/93‑1 (MQ=255)
                                                         gcggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACTTAGCAATGGTGTgaag                                <  2:7509497/93‑1 (MQ=255)
                                                          cggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                               <  2:157914/93‑1 (MQ=255)
                                                          cggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                               <  1:4396068/93‑1 (MQ=255)
                                                             cAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAt                            <  2:6958410/93‑1 (MQ=255)
                                                             cAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAt                            <  1:1299380/93‑1 (MQ=255)
                                                             cAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAt                            >  1:7559348/1‑93 (MQ=255)
                                                              aaTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtt                           <  2:5230144/93‑1 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAACttt                          >  1:4049967/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  1:5206009/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  1:3124207/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  2:1359092/1‑93 (MQ=255)
                                                                tCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                         <  1:1047312/93‑1 (MQ=255)
                                                                tCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                         >  2:6280721/1‑93 (MQ=255)
                                                                  aGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTg                       >  2:4250944/1‑93 (MQ=255)
                                                                   gTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                      <  2:4867183/93‑1 (MQ=255)
                                                                   gTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                      >  2:2518107/1‑93 (MQ=255)
                                                                   gTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                      <  1:2091215/93‑1 (MQ=255)
                                                                    ttCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                     <  2:7858175/93‑1 (MQ=255)
                                                                    ttCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                     <  2:4383123/93‑1 (MQ=255)
                                                                     tCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAtt                    >  1:5107793/1‑93 (MQ=255)
                                                                          tATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGa               <  2:4738903/93‑1 (MQ=255)
                                                                          tATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGa               <  2:6847687/93‑1 (MQ=255)
                                                                           aTCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc              <  1:5784365/93‑1 (MQ=255)
                                                                           aTCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc              <  2:3460454/93‑1 (MQ=255)
                                                                           aTCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc              <  2:7167277/93‑1 (MQ=255)
                                                                            tCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACg             >  1:5563556/1‑93 (MQ=255)
                                                                              cAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGt           >  2:450649/1‑93 (MQ=255)
                                                                              cAGGTTGATAGCAGCCCCACCAGCAGAAAGGCGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACgga           >  2:5738429/1‑92 (MQ=255)
                                                                               aGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCTGACGGTa          >  2:3237009/1‑93 (MQ=255)
                                                                                ggTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTaa         >  2:1734263/1‑93 (MQ=255)
                                                                                    gATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc     <  1:1722957/93‑1 (MQ=255)
                                                                                    gATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc     >  2:1014791/1‑93 (MQ=255)
                                                                                     aTAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACg    >  2:7085520/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902342‑3902521

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