Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,373,180 T→C 5.1% W97R (TGG→CGG)  arnF → undecaprenyl phosphate‑alpha‑L‑ara4N exporter; flippase ArnEF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,373,1800TC 5.1% 301.0 / 4.9 99W97R (TGG→CGG) arnFundecaprenyl phosphate‑alpha‑L‑ara4N exporter; flippase ArnEF subunit
Reads supporting (aligned to +/- strand):  ref base T (54/38);  new base C (2/3);  total (57/42)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CGGATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCTATAACT  >  NC_000913/2373033‑2373328
                                                                                                                                                   |                                                                                                                                                    
cGGATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTg                                                                                                                                                     <  2:1117517/149‑1 (MQ=255)
  gATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTggg                                                                                                                                                   <  1:208064/149‑1 (MQ=255)
      cTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAgg                                                                                                                                               <  2:568832/149‑1 (MQ=255)
      cTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAgg                                                                                                                                               <  1:1533372/149‑1 (MQ=255)
      cTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCAGCCGGGAAgt                                                                                                                                               <  2:1560715/149‑2 (MQ=255)
       tGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAggg                                                                                                                                              >  1:1500946/1‑149 (MQ=255)
       tGTTACCCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAggg                                                                                                                                              >  2:1114098/1‑149 (MQ=255)
        gTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGa                                                                                                                                             <  2:198587/149‑1 (MQ=255)
          tACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCAGCCGGGAAgtacca                                                                                                                                           <  2:1023530/149‑6 (MQ=255)
            cTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACttt                                                                                                                                         >  1:345495/1‑149 (MQ=255)
             tCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACtttt                                                                                                                                        >  1:700657/1‑149 (MQ=255)
              cGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACttttt                                                                                                                                       <  1:1333249/149‑1 (MQ=255)
               gggCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTc                                                                                                                                      >  2:107264/1‑149 (MQ=255)
                ggCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCg                                                                                                                                     <  1:841916/149‑1 (MQ=255)
                ggCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCg                                                                                                                                     <  2:1047234/149‑1 (MQ=255)
                 gctgctGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGc                                                                                                                                    >  1:423520/1‑149 (MQ=255)
                 gctgctGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGc                                                                                                                                    >  1:532731/1‑149 (MQ=255)
                     ctgGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaa                                                                                                                                >  1:98822/1‑149 (MQ=255)
                        ggTTATCTGCTTTCCGTTTTTTGTTGGTATACAACGCTTCCTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGc                                                                                                                             >  2:1180633/1‑149 (MQ=255)
                        ggTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAGGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGc                                                                                                                             >  2:931878/1‑149 (MQ=255)
                        ggTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCAGGTATGGATTGCTTCAATGGTACTTCCCGCCTGGGAAGGGACTTTTTCGCTGAAAGc                                                                                                                             >  2:1180627/1‑149 (MQ=255)
                         gTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACtttt                                                                                                                                        >  2:56646/1‑137 (MQ=255)
                         gTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACtttt                                                                                                                                        <  1:56646/137‑1 (MQ=255)
                          ttATCTGCTTTCAGTTATTTGTGGGTATTAAACGCTTCATAAAGTTGCGCGAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTTGTATGGATTGCTTCATTGGTACTTCACGGCTGGGAAGGGACTTTTTCGTTGAAAGCac                                                                                                                           <  1:252740/149‑1 (MQ=255)
                                gCTTTCCGTTTTTTGTTTGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTgg                                                                                                                     <  2:655234/149‑1 (MQ=255)
                                gCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTgg                                                                                                                     <  2:467024/149‑1 (MQ=255)
                                gCTTTCCGTTTTTGGTTGGTATAAAACGCTTCATAAGCTTGCGCTAAGCAAAGCCTATGCCTTTCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGTACTTTTTCGCTGAAATCACTACTgg                                                                                                                     >  1:800392/1‑149 (MQ=255)
                                  tttCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAGGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGa                                                                                                                   <  1:931878/149‑1 (MQ=255)
                                   ttCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAg                                                                                                                  >  2:230378/1‑149 (MQ=255)
                                        ttttttGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGcc                                                                                                             >  1:717546/1‑149 (MQ=255)
                                        ttttttGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGcc                                                                                                             >  1:717542/1‑149 (MQ=255)
                                        ttttttGTTGGTATAAAACGCGTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGcc                                                                                                             >  1:811298/1‑149 (MQ=255)
                                              gttgGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGtatt                                                                                                       <  2:98822/149‑1 (MQ=255)
                                                 ggTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATg                                                                                                    >  2:359917/1‑149 (MQ=255)
                                                 ggTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATg                                                                                                    >  2:1145364/1‑149 (MQ=255)
                                                     taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCactact                                                                                                                       >  2:1254021/1‑126 (MQ=255)
                                                     taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCactact                                                                                                                       <  1:1254021/126‑1 (MQ=255)
                                                     taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCg                                                                                                >  2:1463782/1‑149 (MQ=255)
                                                     taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCg                                                                                                >  1:417878/1‑149 (MQ=255)
                                                     taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCg                                                                                                >  2:200035/1‑149 (MQ=255)
                                                         aCGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGtt                                                                                            <  2:920778/149‑1 (MQ=255)
                                                            cTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGAt                                                                                         <  2:188621/149‑1 (MQ=255)
                                                             ttCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATg                                                                                        >  1:1136514/1‑149 (MQ=255)
                                                               cATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCt                                                                                      >  2:1150310/1‑149 (MQ=255)
                                                                aTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTTGGAGTAGCCTGTATTATGAGCGGGTTGATGCTg                                                                                     >  1:728912/1‑149 (MQ=255)
                                                                aTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTg                                                                                     >  1:1388005/1‑149 (MQ=255)
                                                                aTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTg                                                                                     >  2:97056/1‑149 (MQ=255)
                                                                 tAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGa                                                                                    <  2:532731/149‑1 (MQ=255)
                                                                    aCTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGAttt                                                                                 <  1:175913/149‑1 (MQ=255)
                                                                     cTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTgg                                                                                                                     >  2:1439017/1‑112 (MQ=255)
                                                                     cTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTgg                                                                                                                     <  1:1439017/112‑1 (MQ=255)
                                                                     cTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGAtttt                                                                                >  1:778280/1‑149 (MQ=255)
                                                                      ttGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTc                                                                               <  2:778280/149‑1 (MQ=255)
                                                                      ttGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTc                                                                               <  2:508642/149‑1 (MQ=255)
                                                                        gcgcTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCt                                                                             <  1:106637/149‑1 (MQ=255)
                                                                         tactgaGCAAGGCATAGGCTTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTg                                                                            >  2:302239/6‑149 (MQ=255)
                                                                         ctcTAAGCAAAGCCTATGCCTGGCTCAGTATGAGTTATGTGCGGGTATGGAGTGCTTCAATGGTCCTTCCCGGCGGGGAAGGGACTTGTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTTATGCTGATTTTCCTg                                                                            <  1:1180633/147‑1 (MQ=255)
                                                                         cgcTAAGCAAAGCCTATGCCTTGCTCAGTTTGAGTTTTGTGCTGGTATGGATTGGTTCAATGGTAATTCCCGGTTTGGAAGGGACTTTTTCGCTGAAATCACTACTGGGAGTATCCTGTATTATGAGCGGTTTGATGCTGATTTTCCTg                                                                            <  1:1180627/149‑1 (MQ=255)
                                                                          gcTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGCGAGTAGCCTGTATAATGAGCGGGTTGATGCTGATTTTCCTGc                                                                           >  1:14541/1‑149 (MQ=255)
                                                                             aaGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCa                                                                        <  2:1326263/149‑1 (MQ=255)
                                                                              aGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCAc                                                                       >  2:72400/1‑149 (MQ=255)
                                                                               gCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACg                                                                      >  1:1478933/1‑149 (MQ=255)
                                                                               gCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACg                                                                      <  1:107264/149‑1 (MQ=255)
                                                                                cAAAGCCTATGCCTTGCTCAGTAGGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGa                                                                     <  2:116827/149‑1 (MQ=255)
                                                                                   aGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACATCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACaa                                                                  <  2:150091/149‑1 (MQ=255)
                                                                                           ccTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGt                                                          >  1:706612/1‑149 (MQ=255)
                                                                                            cTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGtt                                                         >  2:686467/1‑149 (MQ=255)
                                                                                                 tCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGa                                                    >  1:472773/1‑149 (MQ=255)
                                                                                                  cAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAg                                                   <  2:37239/149‑1 (MQ=255)
                                                                                                   aGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGt                                                  >  1:442472/1‑149 (MQ=255)
                                                                                                    gTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGtt                                                 >  2:269966/1‑149 (MQ=255)
                                                                                                       tGAGTTATGTGCTGGTCTGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTc                                              <  2:153861/149‑1 (MQ=255)
                                                                                                       tGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCGGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTc                                              >  2:1106713/1‑149 (MQ=255)
                                                                                                            tATGTGCGGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTTTGAGCGGGTTTATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGc                                         <  2:887190/149‑1 (MQ=255)
                                                                                                               gtgCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCAc                                      >  1:471744/1‑149 (MQ=255)
                                                                                                               gtgCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTACTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCAc                                      <  1:200035/149‑1 (MQ=255)
                                                                                                               gtgCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTAATATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCAc                                      >  1:1069121/1‑149 (MQ=255)
                                                                                                                 gCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACtt                                    >  2:1231475/1‑149 (MQ=255)
                                                                                                                    ggTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCCCTTTac                                 >  2:1443363/1‑149 (MQ=255)
                                                                                                                     gTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTaca                                >  1:1331313/1‑149 (MQ=255)
                                                                                                                     gTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTaca                                >  1:311876/1‑149 (MQ=255)
                                                                                                                     gTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTaca                                <  1:1231475/149‑1 (MQ=255)
                                                                                                                      tATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCCCGACAAAACAACGTTACTGAGTTTTCCCTGCCACGTTacac                               >  2:102407/1‑149 (MQ=255)
                                                                                                                            ttGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGt                         >  2:866311/1‑149 (MQ=255)
                                                                                                                             tGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGtt                        <  2:472773/149‑1 (MQ=255)
                                                                                                                               cTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTcc                      <  2:345495/149‑1 (MQ=255)
                                                                                                                               cTTCAATGGTACTTCCCGGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTcc                      >  1:1052226/1‑149 (MQ=255)
                                                                                                                                ttCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACACGGTTACTGAGTTTTCCCTGCCACTTTACTCTGCCGTTccc                     >  1:239829/1‑149 (MQ=255)
                                                                                                                                ttCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTccc                     >  1:1109982/1‑149 (MQ=255)
                                                                                                                                       gTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGa                                                    <  2:524719/111‑1 (MQ=255)
                                                                                                                                       gTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGa                                                    >  1:524719/1‑111 (MQ=255)
                                                                                                                                        tACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCg             >  2:524116/1‑149 (MQ=255)
                                                                                                                                        cccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTc                 <  2:29124/143‑1 (MQ=255)
                                                                                                                                        cccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTc                 >  1:29124/3‑145 (MQ=255)
                                                                                                                                         aCTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCgg            <  1:230378/149‑1 (MQ=255)
                                                                                                                                           ttCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCCTCgggg          >  1:130414/1‑149 (MQ=255)
                                                                                                                                           ttCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCATGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCgggg          >  1:125365/1‑149 (MQ=255)
                                                                                                                                            tCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGc         <  2:1500946/149‑1 (MQ=255)
                                                                                                                                             cccGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCt        >  2:414900/1‑149 (MQ=255)
                                                                                                                                                   tGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCTATAACt  <  2:423520/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCTATAACT  >  NC_000913/2373033‑2373328

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