Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,340,409 Δ1 bp coding (113/2268 nt) adiA ← arginine decarboxylase, degradative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,340,4080G.99.1% 1149.0 / 2.5 225coding (114/2268 nt)adiAarginine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (119/103);  minor base C (2/0);  total (122/103)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

CGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTCA  >  NC_000913/4340308‑4340504
                                                                                                    |                                                                                                
cGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTg                                                                                                  <  2:5733190/101‑1 (MQ=255)
 gTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTga                                                                                                 >  1:5733571/1‑100 (MQ=255)
  ttttGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTgat                                                                                                >  1:5733573/1‑99 (MQ=255)
   tttGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTgatt                                                                                               <  1:5733154/101‑4 (MQ=255)
    ttGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTgattt                                                                                              <  1:5733136/101‑5 (MQ=255)
    ttGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTgattt                                                                                              >  2:5733576/1‑97 (MQ=255)
     tGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTa                                                                                            <  1:5732853/101‑1 (MQ=255)
     tGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTa                                                                                            >  1:5733578/1‑101 (MQ=255)
      gATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTaa                                                                                           >  1:5733581/1‑101 (MQ=255)
          ttCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCAc                                                                                       <  2:5732940/101‑1 (MQ=255)
           tCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACg                                                                                      >  1:5733584/1‑101 (MQ=255)
             gTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGt                                                                                    >  1:5733586/1‑101 (MQ=255)
             gTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGt                                                                                    <  2:5733290/101‑1 (MQ=255)
              tCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTa                                                                                   >  2:5733587/1‑101 (MQ=255)
               ccGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTaa                                                                                  >  2:5733590/1‑101 (MQ=255)
                cGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAAc                                                                                 >  2:5733593/1‑101 (MQ=255)
                 ggATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACa                                                                                >  1:5733596/1‑101 (MQ=255)
                 ggATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACa                                                                                <  1:5733271/101‑1 (MQ=255)
                 ggATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACa                                                                                >  2:5733595/1‑101 (MQ=255)
                  gATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACAt                                                                               <  2:5733260/101‑1 (MQ=255)
                    tGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACAttt                                                                             <  1:5733237/101‑1 (MQ=255)
                     gTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACAtttt                                                                            <  1:5732974/101‑1 (MQ=255)
                     gTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACAtttt                                                                            <  1:5733144/101‑1 (MQ=255)
                      ttCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTg                                                                           <  1:5733250/101‑1 (MQ=255)
                      ttCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTg                                                                           >  2:5733601/1‑101 (MQ=255)
                      ttCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTg                                                                           <  2:5733234/101‑1 (MQ=255)
                         cATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTg                                                                        <  2:5733323/101‑1 (MQ=255)
                         cATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTg                                                                        >  1:5733603/1‑101 (MQ=255)
                          aTTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTgg                                                                       <  2:5733276/101‑1 (MQ=255)
                           tttGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGc                                                                      <  2:5733314/101‑1 (MQ=255)
                           tttGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGc                                                                      <  1:5733279/101‑1 (MQ=255)
                             tGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCtt                                                                    <  2:5733164/101‑1 (MQ=255)
                             tGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCtt                                                                    >  2:5733605/1‑101 (MQ=255)
                             tGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCtt                                                                    <  2:5733079/101‑1 (MQ=255)
                              gATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTa                                                                   <  1:5733232/101‑1 (MQ=255)
                              gATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTa                                                                   <  1:5733170/101‑1 (MQ=255)
                              gATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTa                                                                   <  2:5733167/101‑1 (MQ=255)
                              gATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTa                                                                   >  2:5733607/1‑101 (MQ=255)
                               aTAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCAGGCTTaa                                                                  >  1:5733608/1‑101 (MQ=255)
                                  gCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGc                                                               <  2:5733118/101‑1 (MQ=255)
                                   cTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTg                                                                                                  >  1:5733612/1‑66 (MQ=255)
                                     gAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATc                                                            >  1:5733615/1‑101 (MQ=255)
                                      aaCATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCt                                                           >  2:5733616/1‑101 (MQ=255)
                                        cATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGc                                                         >  2:5733618/1‑101 (MQ=255)
                                          tCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCa                                                       >  1:5733620/1‑101 (MQ=255)
                                           cAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAg                                                      >  2:5733622/1‑101 (MQ=255)
                                               cAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGc                                                  >  1:5733627/1‑101 (MQ=255)
                                                aGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCt                                                 >  2:5733628/1‑101 (MQ=255)
                                                  tCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCa                                               >  2:5733630/1‑101 (MQ=255)
                                                    aaTGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAAc                                             >  2:5733633/1‑101 (MQ=255)
                                                     aTGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACg                                            <  2:5733306/101‑1 (MQ=255)
                                                     aTGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACg                                            <  2:5733134/101‑1 (MQ=255)
                                                     aTGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACg                                            <  2:5733155/101‑1 (MQ=255)
                                                      tGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACgg                                           >  1:5733636/1‑101 (MQ=255)
                                                        gCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCg                                         >  2:5733638/1‑101 (MQ=255)
                                                         cTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACg                                                                                      <  2:5733045/55‑1 (MQ=255)
                                                          ttCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGtt                                       >  2:5733640/1‑101 (MQ=255)
                                                           tCGTTTGAAGAGAGAATGGCAAAACCATCAGCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTa                                      <  1:5733390/101‑1 (MQ=255)
                                                            cGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTAc                                     <  2:5733215/101‑1 (MQ=255)
                                                             gTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTAcc                                    >  1:5733643/1‑101 (MQ=255)
                                                              tttGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCg                                   <  2:5733273/101‑1 (MQ=255)
                                                              tttGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCg                                   >  1:5733646/1‑101 (MQ=255)
                                                              tttGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCg                                   >  2:5733647/1‑101 (MQ=255)
                                                               ttGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGa                                  >  1:5733649/1‑101 (MQ=255)
                                                                 gAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGAcc                                >  2:5733654/1‑101 (MQ=255)
                                                                 gAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGAcc                                <  1:5733398/101‑1 (MQ=255)
                                                                 gAAGAGAGAATGGCAAAACCATCATCAAAGGAGGAGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGAcc                                >  1:5733653/1‑101 (MQ=255)
                                                                  aaGAGAGAATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGAccc                               <  1:5733445/101‑1 (MQ=255)
                                                                     agagaATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAgg                            >  1:5733656/1‑101 (MQ=255)
                                                                      gagaATGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGgt                           <  2:5733264/101‑1 (MQ=255)
                                                                         aaTGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTc                        >  2:5733659/1‑101 (MQ=255)
                                                                         aaTGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTc                        <  1:5733336/101‑1 (MQ=255)
                                                                          aTGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCt                       >  1:5733660/1‑101 (MQ=255)
                                                                           tGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCtt                      <  2:5733299/101‑1 (MQ=255)
                                                                           tGGCAAAACCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCtt                      <  2:5733240/101‑1 (MQ=255)
                                                                               aaaaCCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATg                  >  2:5733665/1‑101 (MQ=255)
                                                                               aaaaCCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATg                  >  2:5733667/1‑101 (MQ=255)
                                                                                 aaCCATCATCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGc                                          >  1:5733669/1‑75 (MQ=255)
                                                                                    catcatCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGaa             >  1:5733672/1‑101 (MQ=255)
                                                                                    catcatCAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGaa             >  2:5733673/1‑101 (MQ=255)
                                                                                     atcatcAAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGaaa            <  1:5733498/101‑1 (MQ=255)
                                                                                      tcatcaAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAAc           >  1:5733676/1‑101 (MQ=255)
                                                                                       catcaAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACt          <  1:5733294/101‑1 (MQ=255)
                                                                                        atcaAAGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTc         >  2:5733679/1‑101 (MQ=255)
                                                                                           aaaGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCt      >  2:5733681/1‑101 (MQ=255)
                                                                                           aaaGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCt      >  1:5733682/1‑101 (MQ=255)
                                                                                            aaGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCtt     <  1:5733090/101‑1 (MQ=255)
                                                                                            aaGGAGGTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCtt     >  1:5733685/1‑101 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCGGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    <  1:5733331/101‑1 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    >  2:5733688/1‑101 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    >  2:5733693/1‑101 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    >  1:5733691/1‑101 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    <  2:5733321/101‑1 (MQ=255)
                                                                                             aggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCttt    <  2:5733355/101‑1 (MQ=255)
                                                                                              ggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTc   >  1:5733695/1‑101 (MQ=255)
                                                                                              ggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTc   <  2:5733177/101‑1 (MQ=255)
                                                                                              ggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTc   <  1:5733184/101‑1 (MQ=255)
                                                                                              ggaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTc   <  2:5733465/101‑1 (MQ=255)
                                                                                               gaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTca  >  1:5733698/1‑101 (MQ=255)
                                                                                               gaggTGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTca  <  2:5733255/101‑1 (MQ=255)
                                                                                                    |                                                                                                
CGTTTTGATGTTCGTCCGGATGTTCCATTTGATAGCTGAACATCAGGCAGTCAATGGCTTCGTTTGAAGAGAGAATGGCAAAACCATCATCAAAGGAGGTGGATTTAATCACGGTAACATTTTGCTGGCTTAAAGCATCTGCCAGACGCTCAACGGCGTTACCGACCCAGGTGTCTTGATGGAGAAACTCGCTTTCA  >  NC_000913/4340308‑4340504

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