Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,609,489 A→G intergenic (‑83/‑7) yhhS ← / → b3474 putative transport/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,609,4890AG99.4% 550.1 / ‑4.8 177intergenic (‑83/‑7)yhhS/b3474putative transport/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (89/87);  total (89/88)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 186 total aligned reads displayed.

TTTACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCA  >  NC_000913/3609399‑3609578
                                                                                          |                                                                                         
tttACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGc                                                                                         <  1:728887/93‑1 (MQ=255)
tttACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGc                                                                                         <  2:2921750/93‑1 (MQ=255)
tttACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGc                                                                                         >  1:1275614/1‑93 (MQ=255)
  tACCACGTGCTGACGTGAAAATCAGCAGCAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGccc                                                                                       <  2:132047/93‑1 (MQ=255)
    ccACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCa                                                                                     >  2:4206288/1‑93 (MQ=255)
     cACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCaa                                                                                    <  2:3046049/93‑1 (MQ=255)
      aggTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAAt                                                                                   <  1:1836301/91‑1 (MQ=255)
      aCGTGCTGACGTTAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAAt                                                                                   <  1:657271/93‑1 (MQ=255)
      aCGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAAt                                                                                   >  1:633340/1‑93 (MQ=255)
       cGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATg                                                                                  >  1:4065318/1‑93 (MQ=255)
        gTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATgg                                                                                 >  1:2972949/1‑93 (MQ=255)
        gTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATgg                                                                                 >  1:971615/1‑93 (MQ=255)
         tGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGa                                                                                >  1:487663/1‑93 (MQ=255)
           cTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCTTGACTCCGGCTACCCTTAATACAAAAAGTCACAAGGAGGCCCCAATGGaaa                                                                              <  2:3644448/93‑1 (MQ=255)
           cTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGaaa                                                                              <  1:2778743/93‑1 (MQ=255)
           cTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGaaa                                                                              >  2:1796673/1‑93 (MQ=255)
            tGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAc                                                                             <  1:4325584/93‑1 (MQ=255)
             gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcc                                                                            >  2:3909400/1‑93 (MQ=255)
             gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcc                                                                            <  1:3181593/93‑1 (MQ=255)
               cGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcccc                                                                          <  1:4532218/93‑1 (MQ=255)
                gTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCt                                                                         >  2:3541887/1‑93 (MQ=255)
                  gAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCa                                                                       <  2:2128389/93‑1 (MQ=255)
                  gAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCa                                                                       >  2:2724635/1‑93 (MQ=255)
                   aaaaTCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCaa                                                                      <  2:4473690/93‑1 (MQ=255)
                    aaaTCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAAc                                                                     <  2:3198622/93‑1 (MQ=255)
                       tCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCg                                                                  >  2:2421879/1‑93 (MQ=255)
                       tCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCg                                                                  >  2:2526416/1‑93 (MQ=255)
                       tCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCg                                                                  <  2:4622925/93‑1 (MQ=255)
                         agcagTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGAt                                                                <  2:1220680/93‑1 (MQ=255)
                             gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAATGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                            <  1:3306546/93‑1 (MQ=255)
                             gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                            <  1:806691/93‑1 (MQ=255)
                             gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                            <  1:1068551/93‑1 (MQ=255)
                             gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                            >  2:656421/1‑93 (MQ=255)
                             atAAGAATCAGAATATTGCTGGCGTAACTCCCGCTACACTTAATACAAAAAGCCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                            <  2:3974734/92‑1 (MQ=255)
                              tAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAAc                                                           <  2:2520786/93‑1 (MQ=255)
                               aaGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACg                                                          <  1:4644063/93‑1 (MQ=255)
                                aGAAGCAGAATATTGCCGGCGGCCCTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACgg                                                         <  1:1026883/93‑1 (MQ=255)
                                   aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                      <  1:1864584/93‑1 (MQ=255)
                                   aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                      >  1:2532343/1‑93 (MQ=255)
                                   aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                      <  1:3556837/93‑1 (MQ=255)
                                   aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                      >  2:1433529/1‑93 (MQ=255)
                                      aGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAATCCCCTCAACCCGATAAAACGGGCATGc                                                   <  2:1869142/93‑1 (MQ=255)
                                       gAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGca                                                  >  1:3591272/1‑93 (MQ=255)
                                         atatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGcaca                                                >  1:1566100/1‑93 (MQ=255)
                                         atatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGcaca                                                <  1:2717752/93‑1 (MQ=255)
                                         atatTCCTGGCGTGACTCCCGCTACACTTAATACAACAAGTCACAAGGAGGCCCCACTGGCAAGCCCTCCACACGATAAACCGGGCATGcaca                                                >  2:3697835/1‑93 (MQ=255)
                                          tatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACAt                                               <  1:4432144/93‑1 (MQ=255)
                                          tatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACAt                                               <  1:789311/93‑1 (MQ=255)
                                           atTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACAtt                                              <  2:4226665/93‑1 (MQ=255)
                                            ttGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTc                                             >  1:1244999/1‑93 (MQ=255)
                                            ttGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTc                                             >  1:2541669/1‑93 (MQ=255)
                                             tGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCt                                            <  1:578873/93‑1 (MQ=255)
                                               gtGGCGGGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGc                                          <  2:1220688/92‑1 (MQ=255)
                                                   cGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCaa                                      <  1:2118187/93‑1 (MQ=255)
                                                    gTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAg                                     <  2:2372813/93‑1 (MQ=255)
                                                     tGACGCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGc                                    <  2:886017/93‑1 (MQ=255)
                                                       aCTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTg                                  >  1:1713962/1‑93 (MQ=255)
                                                        cTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTgg                                 >  1:3961414/1‑93 (MQ=255)
                                                         actcGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAATCCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGc                                <  1:3648208/90‑1 (MQ=255)
                                                          cccGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGcc                               <  1:1452613/93‑1 (MQ=255)
                                                          cccGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGcc                               >  2:152868/1‑93 (MQ=255)
                                                           ccGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCt                              >  2:3131919/1‑93 (MQ=255)
                                                           ccGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCt                              >  2:2904377/1‑93 (MQ=255)
                                                             gCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCg                            >  1:4259703/1‑93 (MQ=255)
                                                             gCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCg                            >  1:3346521/1‑93 (MQ=255)
                                                               tACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCt                          >  2:1891398/1‑93 (MQ=255)
                                                               tACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCt                          >  2:2198877/1‑93 (MQ=255)
                                                               tACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCt                          >  1:4573548/1‑93 (MQ=255)
                                                                acacTTAATACAAAAAGTCACAAGGAGGCCCTAATGGAACCCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTg                         <  2:4102434/93‑1 (MQ=255)
                                                                acacTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTg                         >  1:1728310/1‑93 (MQ=255)
                                                                acacTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTg                         <  2:3155770/93‑1 (MQ=255)
                                                                 cacTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTgg                        <  1:2388838/93‑1 (MQ=255)
                                                                  acTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGt                       >  1:2211232/1‑93 (MQ=255)
                                                                   cTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTa                      <  1:4301090/93‑1 (MQ=255)
                                                                    ttAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAg                     <  1:4399556/93‑1 (MQ=255)
                                                                    ttAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAg                     >  1:218617/1‑93 (MQ=255)
                                                                    ttAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAg                     >  1:3826503/1‑93 (MQ=255)
                                                                     tAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGt                    >  2:3274768/1‑93 (MQ=255)
                                                                     tAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGt                    <  2:348680/93‑1 (MQ=255)
                                                                     tAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGt                    >  2:944125/1‑93 (MQ=255)
                                                                     tAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGt                    >  1:2266810/1‑93 (MQ=255)
                                                                        tACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGAt                 >  1:2402007/1‑93 (MQ=255)
                                                                        aaCAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGAt                 <  2:2785201/92‑1 (MQ=255)
                                                                         aCAAAAAGTCACAAGGAGGCCCCAATGGAAACCCGTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATc                <  1:110414/93‑1 (MQ=255)
                                                                          cAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAATCGGGCATGCCCATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATcc               <  1:343165/93‑1 (MQ=255)
                                                                           aaaaaGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCt              <  1:4242083/93‑1 (MQ=255)
                                                                           aaaaaGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCt              >  2:4221123/1‑93 (MQ=255)
                                                                               aGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTcgcg          >  2:55876/1‑93 (MQ=255)
                                                                                gTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCgg         >  1:1362435/1‑93 (MQ=255)
                                                                                 tCACAAGGAGGCCCCAATGGAAGCCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCggg        <  2:1358122/93‑1 (MQ=255)
                                                                                  cacaAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGc       >  1:2628175/1‑93 (MQ=255)
                                                                                    caAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAt     <  1:215264/93‑1 (MQ=255)
                                                                                    caAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAt     >  2:2586455/1‑93 (MQ=255)
                                                                                     aaGGAGGCCCCAATGGAACCCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAtt    <  2:4051052/93‑1 (MQ=255)
                                                                                     aaGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAtt    >  1:4189684/1‑93 (MQ=255)
                                                                                     aaGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAtt    >  2:2069752/1‑93 (MQ=255)
                                                                                      aggaggCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc   >  2:115237/1‑93 (MQ=255)
                                                                                       ggAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa  <  1:3909400/93‑1 (MQ=255)
                                                                                       ggAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa  <  1:1417670/93‑1 (MQ=255)
                                                                                       ggAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa  <  1:192959/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCA  >  NC_000913/3609399‑3609578

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