Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,807,262 A→G S366G (AGC→GGC)  kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,807,2620AG99.6% 778.1 / NA 226S366G (AGC→GGC) kdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (115/110);  minor base T (0/1);  total (115/111)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGC  >  NC_000913/3807172‑3807351
                                                                                          |                                                                                         
ccGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCaggcaggc                                                                                         >  2:5668549/1‑93 (MQ=255)
ccGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCaggcaggc                                                                                         <  1:5939346/93‑1 (MQ=255)
 cGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCg                                                                                        <  2:6335944/93‑1 (MQ=255)
  gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCgg                                                                                       >  2:1136342/1‑93 (MQ=255)
  gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCgg                                                                                       >  1:3702838/1‑93 (MQ=255)
     ggAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCt                                                                                    <  2:570721/93‑1 (MQ=255)
     ggAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCt                                                                                    >  2:2708224/1‑93 (MQ=255)
     ggAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCt                                                                                    <  2:2614252/93‑1 (MQ=255)
     ggAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCt                                                                                    <  2:2383425/93‑1 (MQ=255)
       aaGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGa                                                                                  <  1:68034/93‑1 (MQ=255)
        aGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAt                                                                                 <  2:19052/93‑1 (MQ=255)
        aGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAt                                                                                 <  2:1535422/93‑1 (MQ=255)
         gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAtt                                                                                >  2:4861160/1‑93 (MQ=255)
            gCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttacc                                                                             >  2:5459276/1‑93 (MQ=255)
            gCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttacc                                                                             <  2:2984036/93‑1 (MQ=255)
            gCCGCACACGCTATTCCGGAATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttacc                                                                             >  1:6306520/1‑93 (MQ=255)
              cGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgt                                                                           <  2:388182/93‑1 (MQ=255)
               gCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTT‑‑GCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttac                                                                        >  1:4951637/1‑93 (MQ=255)
                cacaCGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgtta                                                                         <  1:3654443/93‑1 (MQ=255)
                 acacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGGTGAttaccgttac                                                                        >  1:3537945/1‑93 (MQ=255)
                 acacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttac                                                                        >  1:1168884/1‑93 (MQ=255)
                 acacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttac                                                                        <  1:828658/93‑1 (MQ=255)
                  cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttacc                                                                       >  1:2970501/1‑93 (MQ=255)
                  cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttacc                                                                       >  2:5928801/1‑93 (MQ=255)
                  cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttacc                                                                       <  1:1565632/93‑1 (MQ=255)
                  cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttacccttacc                                                                       >  1:4616570/1‑93 (MQ=255)
                   acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGAttaccgttaccg                                                                      >  1:5990912/1‑93 (MQ=255)
                    cGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGa                                                                     <  1:3159981/93‑1 (MQ=255)
                     gCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGAt                                                                    >  2:314069/1‑93 (MQ=255)
                       tATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGc                                                                  <  1:2018261/93‑1 (MQ=255)
                       tATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGc                                                                  <  1:3118688/93‑1 (MQ=255)
                         ttCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCa                                                                >  1:323804/1‑93 (MQ=255)
                             ggTATTGATGGGGCCGCATACTTTTATCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGACTGATTACCGTTACCGATGCCACTAc                                                            <  1:4264687/93‑1 (MQ=255)
                               tATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGc                                                          <  1:1793323/93‑1 (MQ=255)
                                aTTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCt                                                         <  1:2827465/93‑1 (MQ=255)
                                aTTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCt                                                         <  2:2050638/93‑1 (MQ=255)
                                 ttGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCtt                                                        <  1:6053179/93‑1 (MQ=255)
                                   gATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGc                                                      >  1:1579130/1‑93 (MQ=255)
                                     tGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCaa                                                    >  1:2542082/1‑93 (MQ=255)
                                       gggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCaaaa                                                  >  2:4979507/1‑93 (MQ=255)
                                       gggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCaaaa                                                  >  2:106878/1‑93 (MQ=255)
                                       gggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCaaaa                                                  >  1:1329548/1‑93 (MQ=255)
                                        ggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGTCTGATTACCGTTACCGATGCCACTACGCTTGCAAAag                                                 <  2:3931085/93‑1 (MQ=255)
                                         gCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAaga                                                <  2:4162840/93‑1 (MQ=255)
                                          ccGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAagag                                               >  2:4518343/1‑93 (MQ=255)
                                          ccGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAagag                                               <  1:5226046/93‑1 (MQ=255)
                                            gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTGCCGATGCCACTACGCTTGCAAAAGAGGt                                             <  2:428435/93‑1 (MQ=255)
                                            gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGt                                             <  1:414297/93‑1 (MQ=255)
                                            gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGt                                             >  1:5409259/1‑93 (MQ=255)
                                              aTACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGttt                                           >  1:5716472/1‑93 (MQ=255)
                                              aTACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGttt                                           >  2:3185967/1‑93 (MQ=255)
                                               tACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                          >  1:1479179/1‑93 (MQ=255)
                                               tACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                          <  1:5528690/93‑1 (MQ=255)
                                                  ttttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                       >  2:2533052/1‑93 (MQ=255)
                                                  ttttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                       >  1:381111/1‑93 (MQ=255)
                                                  ttttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                       <  1:2255664/93‑1 (MQ=255)
                                                   tttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctct                                      >  2:2667147/1‑93 (MQ=255)
                                                   tttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctct                                      <  1:1711485/93‑1 (MQ=255)
                                                    ttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCtt                                     <  2:3366950/93‑1 (MQ=255)
                                                       aCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAc                                  >  1:6016561/1‑93 (MQ=255)
                                                        cTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACt                                 >  2:3990929/1‑93 (MQ=255)
                                                          ttAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTNTTTACTCa                               >  1:3752164/1‑93 (MQ=255)
                                                          ttAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCa                               >  1:5817093/1‑93 (MQ=255)
                                                           tAAAGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAc                              >  1:2411122/1‑93 (MQ=255)
                                                            aaaGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAcc                             <  2:2280411/93‑1 (MQ=255)
                                                             aaGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCg                            >  1:648677/1‑93 (MQ=255)
                                                             aaGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCg                            <  2:549113/93‑1 (MQ=255)
                                                             aaGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCg                            >  2:3868717/1‑93 (MQ=255)
                                                              aGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGa                           <  2:4058435/93‑1 (MQ=255)
                                                              aGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGa                           >  2:142813/1‑93 (MQ=255)
                                                              aGACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGa                           >  2:4714005/1‑93 (MQ=255)
                                                               gACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGAc                          <  2:6176623/93‑1 (MQ=255)
                                                               gACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGAc                          >  2:2893272/1‑93 (MQ=255)
                                                               gACATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGAc                          >  1:2350098/1‑93 (MQ=255)
                                                               gACATGTGCGCGCGGCTGGAGCAGGCAGGCGGGCTAATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGATTTCCTCTTTACTCACCGAc                          <  1:5036396/93‑1 (MQ=255)
                                                                 cATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGc                        <  1:3294512/93‑1 (MQ=255)
                                                                 cATTTGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGc                        >  1:4069705/1‑93 (MQ=255)
                                                                   tttGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCg                      >  2:1705871/1‑93 (MQ=255)
                                                                    ttGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGa                     <  1:4058953/93‑1 (MQ=255)
                                                                     tGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGAt                    <  2:4840933/93‑1 (MQ=255)
                                                                     tGCGCGCGGCTGGAGCAGGCAGGCGGGCTGATTACCGGTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGAt                    >  1:4984039/1‑93 (MQ=255)
                                                                       cgcgcgGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                  <  2:1355492/93‑1 (MQ=255)
                                                                        gcgcgGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAc                 <  2:5943715/93‑1 (MQ=255)
                                                                         cgcgGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAcc                <  1:5139442/93‑1 (MQ=255)
                                                                         cgcgGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAcc                <  1:1335895/93‑1 (MQ=255)
                                                                         cgcgGCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGGTTCCTCTTTACTCACCGACGCCGATTAcc                >  1:1913817/1‑93 (MQ=255)
                                                                          gcgGCTGGAGCAGGCAGGCGGTCTTATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCg               <  1:4504018/93‑1 (MQ=255)
                                                                            ggCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTa             >  2:3693822/1‑93 (MQ=255)
                                                                             gCTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAg            >  1:5840402/1‑93 (MQ=255)
                                                                              cTGGAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGt           <  1:3353600/93‑1 (MQ=255)
                                                                              cTGGAGCAGGCAGGCGGCCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACCCCGATTACCGTAGt           <  1:5587341/93‑1 (MQ=255)
                                                                                 gAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc        >  2:3639892/1‑93 (MQ=255)
                                                                                 gAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc        <  2:2807561/93‑1 (MQ=255)
                                                                                 gAGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc        <  1:145709/93‑1 (MQ=255)
                                                                                  aGCAGGCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCGTCTTTACTCACCGACGCCGCTTACAGTAGTTTCt       >  2:538700/1‑93 (MQ=255)
                                                                                   gcaggcagGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa      <  1:771637/93‑1 (MQ=255)
                                                                                   gcaggcagGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa      <  1:2133023/93‑1 (MQ=255)
                                                                                     aggcaggcGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATg    <  2:2143885/93‑1 (MQ=255)
                                                                                      ggCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATgg   <  1:4688473/93‑1 (MQ=255)
                                                                                      ggCAGGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATgg   >  2:5101688/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGC  >  NC_000913/3807172‑3807351

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