Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,593 A→G K317K (AAA→AAG wzzE → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,5930AG92.2% 1561.4 / 94.6 549K317K (AAA→AAGwzzEputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (26/17);  new base G (279/227);  total (305/244)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 561 total aligned reads displayed.

GGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGC  >  NC_000913/3967502‑3967683
                                                                                           |                                                                                            
gggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAAc                                                                                             >  1:2652192/1‑93 (MQ=255)
  gCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAgcgc                                                                                           >  2:6431732/1‑93 (MQ=255)
  gCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAgcgc                                                                                           <  2:1620810/93‑1 (MQ=255)
  gCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAgcgc                                                                                           >  2:2350677/1‑93 (MQ=255)
    cATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGa                                                                                         >  2:891101/1‑93 (MQ=255)
    cATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGa                                                                                         <  1:800487/93‑1 (MQ=255)
      tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATa                                                                                       >  1:3340752/1‑93 (MQ=255)
      tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATa                                                                                       >  1:2510946/1‑93 (MQ=255)
      tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATa                                                                                       >  2:6872766/1‑93 (MQ=255)
      tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATa                                                                                       >  2:1491940/1‑93 (MQ=255)
      tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATa                                                                                       >  2:6648776/1‑93 (MQ=255)
       gTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAg                                                                                      <  2:1354715/93‑1 (MQ=255)
       gTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAg                                                                                      >  2:5413634/1‑93 (MQ=255)
        ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGc                                                                                     <  2:7140580/93‑1 (MQ=255)
           aaCACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCa                                                                                  >  2:3322373/1‑93 (MQ=255)
           aaCACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCa                                                                                  <  2:1945987/93‑1 (MQ=255)
              acCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCAcgt                                                                               >  2:2769202/1‑93 (MQ=255)
              acCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCAcgt                                                                               <  1:5429121/93‑1 (MQ=255)
              acCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCAcgt                                                                               <  2:5730183/93‑1 (MQ=255)
              acCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCAcgt                                                                               <  2:4832467/93‑1 (MQ=255)
                 cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCAcgtcgt                                                                            >  1:1035539/1‑93 (MQ=255)
                  tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTg                                                                          >  1:1693521/1‑93 (MQ=255)
                   gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGc                                                                         >  2:7270656/1‑93 (MQ=255)
                   gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGc                                                                         >  1:1330050/1‑93 (MQ=255)
                   gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGc                                                                         >  1:4347375/1‑93 (MQ=255)
                     aTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCt                                                                       >  2:178667/1‑93 (MQ=255)
                     aTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCt                                                                       >  1:7288297/1‑93 (MQ=255)
                     aTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCt                                                                       <  1:2996073/93‑1 (MQ=255)
                       gTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTc                                                                     <  2:5262356/93‑1 (MQ=255)
                       gTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTc                                                                     <  2:6981207/93‑1 (MQ=255)
                       gTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTc                                                                     >  1:6614272/1‑93 (MQ=255)
                        ttGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTcc                                                                    <  1:3237963/93‑1 (MQ=255)
                        ttGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTcc                                                                    >  2:4639012/1‑93 (MQ=255)
                         tGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCt                                                                   >  2:5870379/1‑93 (MQ=255)
                               aaCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTgatgat                                                             >  1:1553974/1‑93 (MQ=255)
                                aCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGAtt                                                            >  2:6572659/1‑93 (MQ=255)
                                 cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTa                                                           >  1:2773684/1‑93 (MQ=255)
                                  ccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAt                                                          <  1:4831067/93‑1 (MQ=255)
                                  ccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAt                                                          <  1:5540730/93‑1 (MQ=255)
                                  ccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAt                                                          <  1:6960024/93‑1 (MQ=255)
                                   cTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGGCCTTCCTGATGATTAt                                                          >  1:6146536/1‑93 (MQ=255)
                                   cTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAtg                                                         >  1:3962763/1‑93 (MQ=255)
                                   cTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAtg                                                         >  1:64695/1‑93 (MQ=255)
                                    tGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAtgt                                                        >  1:3575795/1‑93 (MQ=255)
                                    tGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAtgt                                                        >  1:4996064/1‑93 (MQ=255)
                                     ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTAtgtg                                                       <  2:6371650/93‑1 (MQ=255)
                                      gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTgg                                                      >  2:2402794/1‑93 (MQ=255)
                                      gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTgg                                                      <  2:2949961/93‑1 (MQ=255)
                                      gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTgg                                                      <  1:4586541/93‑1 (MQ=255)
                                       aTCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTggg                                                     <  2:2162897/93‑1 (MQ=255)
                                        tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTgggg                                                    >  1:5783103/1‑93 (MQ=255)
                                        tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTgggg                                                    >  2:4466365/1‑93 (MQ=255)
                                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAt                                                 >  1:6839798/1‑93 (MQ=255)
                                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAt                                                 >  1:2414765/1‑93 (MQ=255)
                                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAt                                                 <  2:4299985/93‑1 (MQ=255)
                                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAt                                                 <  1:1411049/93‑1 (MQ=255)
                                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAt                                                 >  1:950049/1‑93 (MQ=255)
                                            cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAtt                                                <  2:2317821/93‑1 (MQ=255)
                                             gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTg                                               >  2:1703775/1‑93 (MQ=255)
                                             gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTg                                               >  2:6141049/1‑93 (MQ=255)
                                              ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGt                                              <  2:4904/93‑1 (MQ=255)
                                               tttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTc                                             <  1:6314742/93‑1 (MQ=255)
                                                ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCg                                            <  1:7138907‑M1/93‑2 (MQ=255)
                                                 tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCgg                                           <  1:1740811‑M1/93‑3 (MQ=255)
                                                 tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCgg                                           >  1:5703893‑M1/1‑91 (MQ=255)
                                                   aGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCgggg                                         >  2:1868130‑M1/1‑89 (MQ=255)
                                                       cTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCt                                    >  1:5190013/1‑93 (MQ=255)
                                                        tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggc                                    <  1:3515971‑M1/93‑10 (MQ=255)
                                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGa                                  <  2:7402844/93‑1 (MQ=255)
                                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGa                                  >  2:2084833/1‑93 (MQ=255)
                                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGa                                  <  2:5584834/93‑1 (MQ=255)
                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGAt                                 <  1:3176595/93‑1 (MQ=255)
                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGAt                                 <  1:5149867/93‑1 (MQ=255)
                                                            gCTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgat                                >  1:1144801‑M1/1‑80 (MQ=255)
                                                            gCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCg                               <  1:4739804/93‑1 (MQ=255)
                                                             cTATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatc                               >  2:4150619‑M1/1‑79 (MQ=255)
                                                              tATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCgg                              >  2:3827196/1‑93 (MQ=255)
                                                              tATTTGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCggg                             <  1:5991660/93‑1 (MQ=255)
                                                                tttGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggg                            >  2:5008035‑M1/1‑76 (MQ=255)
                                                                tttGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggg                            <  1:3850831‑M1/93‑18 (MQ=255)
                                                                tttGCGTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGc                           >  2:4010433/1‑93 (MQ=255)
                                                                     gTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGgtg                      >  2:4553527/1‑93 (MQ=255)
                                                                     gTACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGgtg                      >  1:3684046/1‑93 (MQ=255)
                                                                      tACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggt                      <  1:4088705‑M1/93‑24 (MQ=255)
                                                                      tACGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGgtgt                     >  2:5322317/1‑93 (MQ=255)
                                                                       aCGCCGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtg                     >  2:2540815‑M1/1‑69 (MQ=255)
                                                                           cGGAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCAt                <  2:4974227/93‑1 (MQ=255)
                                                                            ggAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCAtt               >  2:125498/1‑93 (MQ=255)
                                                                            ggAAGAACCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCAtt               <  1:2897334/93‑1 (MQ=255)
                                                                             gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTa              >  1:3078682/1‑93 (MQ=255)
                                                                              aagaaCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTaa             >  1:5340279/1‑93 (MQ=255)
                                                                               agaaCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtgtcgcatta             >  1:217615‑M1/1‑61 (MQ=255)
                                                                               agaaCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtgtcgcatta             <  1:202561‑M1/93‑33 (MQ=255)
                                                                                gaaCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAAcc           <  2:4756804/93‑1 (MQ=255)
                                                                                  aCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtgtcgcattaacc          <  2:5990890‑M1/93‑36 (MQ=255)
                                                                                   ccGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtgtcgcattaaccc         >  2:7511573‑M1/1‑57 (MQ=255)
                                                                                     ggTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg      >  1:4189952/1‑93 (MQ=255)
                                                                                      gTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCggggggggctgatcggggctggtgtcgcattaacccgcc      >  2:3534008‑M1/1‑54 (MQ=255)
                                                                                      gTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGt     >  1:4451339/1‑93 (MQ=255)
                                                                                       tAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCAGATCGGGGCTGGAGCCGCATTAACCCGCCGtt    >  1:1844637/1‑93 (MQ=255)
                                                                                           |                                                                                            
GGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGC  >  NC_000913/3967502‑3967683

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