Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,795,982 A→G S320P (TCC→CCC)  rfaK ← putative hexose transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,795,9820AG98.9% 855.8 / ‑1.7 268S320P (TCC→CCC) rfaKputative hexose transferase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (164/101);  total (167/101)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 279 total aligned reads displayed.

TCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAA  >  NC_000913/3795892‑3796071
                                                                                          |                                                                                         
tctTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGcc                                                                                         >  1:1042704/1‑93 (MQ=255)
tctTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGcc                                                                                         >  1:488560/1‑93 (MQ=255)
  tttACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCaa                                                                                       >  1:4081969/1‑93 (MQ=255)
   ttACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCaaa                                                                                      <  1:1383827/93‑1 (MQ=255)
   ttACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCaaa                                                                                      <  1:3414124/93‑1 (MQ=255)
   ttACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCaaa                                                                                      >  1:3886823/1‑93 (MQ=255)
   ttACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAACGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCaaa                                                                                      >  2:1037652/1‑93 (MQ=255)
    tACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAAt                                                                                     <  2:4558807/93‑1 (MQ=255)
    tACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAAt                                                                                     <  2:1981818/93‑1 (MQ=255)
     aCTTAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAAtt                                                                                    >  1:3389253/1‑93 (MQ=255)
      cTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTa                                                                                   >  2:2927431/1‑93 (MQ=255)
       tGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGTGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTaa                                                                                  >  1:1926375/1‑93 (MQ=255)
       tGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTaa                                                                                  <  1:4120908/93‑1 (MQ=255)
        gAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTaaa                                                                                 >  2:3948825/1‑93 (MQ=255)
          aTTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATg                                                                               >  1:1560109/1‑93 (MQ=255)
          aTTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATg                                                                               >  1:287176/1‑93 (MQ=255)
           ttAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGc                                                                              >  1:3720512/1‑93 (MQ=255)
            tAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCg                                                                             >  1:3294042/1‑93 (MQ=255)
              gTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAg                                                                           >  2:1196336/1‑93 (MQ=255)
              gTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAg                                                                           <  1:431301/93‑1 (MQ=255)
              gTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAg                                                                           >  2:1571155/1‑93 (MQ=255)
               ttAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAgg                                                                          >  2:1832246/1‑93 (MQ=255)
                 aaaGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGa                                                                        >  1:2053670/1‑93 (MQ=255)
                   aGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGata                                                                      <  2:4710063/93‑1 (MQ=255)
                    gTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGatat                                                                     >  1:4603770/1‑93 (MQ=255)
                     tttCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATAtt                                                                    >  1:1532613/1‑93 (MQ=255)
                     tttCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATAtt                                                                    >  1:2244082/1‑93 (MQ=255)
                       tCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTcc                                                                  <  2:3563226/93‑1 (MQ=255)
                         aGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCgg                                                                >  2:3281750/1‑93 (MQ=255)
                         aGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATCAAATGCGAGGGATATTCCgg                                                                >  1:2146697/1‑93 (MQ=255)
                          gTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGa                                                               >  1:2749863/1‑93 (MQ=255)
                          gTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGa                                                               >  2:2779619/1‑93 (MQ=255)
                           tATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAg                                                              >  1:797699/1‑93 (MQ=255)
                           tATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAg                                                              <  1:1761111/93‑1 (MQ=255)
                           tATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAg                                                              >  2:4008474/1‑93 (MQ=255)
                             tCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTa                                                            >  2:2755951/1‑93 (MQ=255)
                             tCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTa                                                            >  2:131755/1‑93 (MQ=255)
                             tCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTa                                                            <  2:3684128/93‑1 (MQ=255)
                              cGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTAc                                                           <  2:4312668/93‑1 (MQ=255)
                                atatCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGt                                                         <  1:4674776/93‑1 (MQ=255)
                                  atCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGa                                                       >  2:1092752/1‑93 (MQ=255)
                                   tCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGaa                                                      <  2:3876537/93‑1 (MQ=255)
                                       ttACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAg                                                  <  2:2827066/93‑1 (MQ=255)
                                           cGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTa                                              >  1:1197364/1‑93 (MQ=255)
                                            gTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaa                                             >  1:2114697/1‑93 (MQ=255)
                                            gTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaa                                             <  2:350963/93‑1 (MQ=255)
                                            gTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGCGAATTCGGGTaa                                             >  2:3014813/1‑93 (MQ=255)
                                             tataTGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaaa                                            >  2:1909748/1‑93 (MQ=255)
                                             tataTGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaaa                                            >  2:1358653/1‑93 (MQ=255)
                                              atatGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaaaa                                           <  1:4189524/93‑1 (MQ=255)
                                              atatGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaaaa                                           <  1:598391/93‑1 (MQ=255)
                                              atatGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTaaaa                                           >  1:4037452/1‑93 (MQ=255)
                                                 tGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGc                                        >  2:3853570/1‑93 (MQ=255)
                                                 tGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGc                                        <  2:4581509/93‑1 (MQ=255)
                                                    tGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAg                                     <  2:4089695/93‑1 (MQ=255)
                                                    tGGGGATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAg                                     <  2:3790133/93‑1 (MQ=255)
                                                      gggATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGTGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCg                                   >  1:3254050/1‑93 (MQ=255)
                                                      gggATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCg                                   >  1:1070477/1‑93 (MQ=255)
                                                        gATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGtt                                 <  2:874702/93‑1 (MQ=255)
                                                        gATACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGtt                                 >  1:3057429/1‑93 (MQ=255)
                                                         aTACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTg                                <  2:4129093/93‑1 (MQ=255)
                                                          tACTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTgg                               <  1:2289478/93‑1 (MQ=255)
                                                           aCTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGt                              >  2:904756/1‑93 (MQ=255)
                                                           aCTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGt                              <  1:2161817/93‑1 (MQ=255)
                                                            cTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGtt                             >  1:2699955/1‑93 (MQ=255)
                                                            cTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGtt                             <  2:3441732/93‑1 (MQ=255)
                                                            cTATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGtt                             <  1:3869941/93‑1 (MQ=255)
                                                             tATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGTTAAAATGCAAGCTTTGGttt                            >  1:1607531/1‑93 (MQ=255)
                                                             tATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGttt                            >  1:2579650/1‑93 (MQ=255)
                                                             tATTTGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGttt                            >  1:4094889/1‑93 (MQ=255)
                                                             tATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGAAATTCCGGAGTACGTGAATTAGGGTAACATGCAAGCGTTGGttt                            >  1:3669499/1‑93 (MQ=255)
                                                               tttGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTAt                          >  1:3638451/1‑93 (MQ=255)
                                                               tttGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTAt                          >  2:3344625/1‑93 (MQ=255)
                                                               tttGTATAGATTTATGATGATTCGGGGGCAAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTag                          >  1:844284/1‑92 (MQ=255)
                                                                ttGTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATg                         <  1:3538459/93‑1 (MQ=255)
                                                                  gTATAGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGAt                       >  2:2787557/1‑93 (MQ=255)
                                                                    ataGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATaa                     >  1:4372455/1‑93 (MQ=255)
                                                                    ataGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATaa                     <  1:3164718/93‑1 (MQ=255)
                                                                     taGATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAg                    <  1:4423837/93‑1 (MQ=255)
                                                                       gATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGGTATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGcc                  <  2:2716402/93‑1 (MQ=255)
                                                                       gATTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGAACGTGCCTTAGGGAATAATGCAAGCATTGGTATATGCTAAGcc                  >  1:1409351/1‑93 (MQ=255)
                                                                        aTTTATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCg                 <  2:1460320/93‑1 (MQ=255)
                                                                           tATGATGATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAg              >  2:194329/1‑93 (MQ=255)
                                                                            atgatgATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGc             >  1:1282349/1‑93 (MQ=255)
                                                                            atgatgATTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGc             >  2:2551315/1‑93 (MQ=255)
                                                                             tgatgaTTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCg            >  2:2294710/1‑93 (MQ=255)
                                                                             tgatgaTTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCg            >  2:3084407/1‑93 (MQ=255)
                                                                             tgatgaTTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCg            >  1:716908/1‑93 (MQ=255)
                                                                             tgatgaTTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCg            >  1:3943427/1‑93 (MQ=255)
                                                                              gatgatTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGa           <  2:3227837/93‑1 (MQ=255)
                                                                               atgatTCGGGGGACAAATTAAATGCGAGGGATATTCCGAAGTACGTGAATTAAGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGat          <  1:4223177/93‑1 (MQ=255)
                                                                                tgatTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGata         <  2:1133747/93‑1 (MQ=255)
                                                                                 gatTCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGatat        >  1:3216908/1‑93 (MQ=255)
                                                                                    tCGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGaa     <  1:4166881/93‑1 (MQ=255)
                                                                                     cGGGGGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAAc    >  2:1406493/1‑93 (MQ=255)
                                                                                      ggggGCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACa   >  2:1734305/1‑93 (MQ=255)
                                                                                       ggggCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACaa  <  1:4247306/93‑1 (MQ=255)
                                                                                       ggggCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACaa  >  1:1827554/1‑93 (MQ=255)
                                                                                       ggggCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACaa  <  2:2335634/93‑1 (MQ=255)
                                                                                       ggggCCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACaa  >  2:3140796/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAA  >  NC_000913/3795892‑3796071

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