Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,436 G→A 5.8% intergenic (+23/+142) cmoB → / ← torZ tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent/trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,4360GA 5.8% 349.5 / 8.3 104intergenic (+23/+142)cmoB/torZtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (45/53);  new base A (3/3);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 107 total aligned reads displayed.

AGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATT  >  NC_000913/1954289‑1954571
                                                                                                                                                   |                                                                                                                                       
agcGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATgc                                                                                                                                        >  2:1580448/1‑149 (MQ=255)
  cGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATgcgg                                                                                                                                      <  4:830974/149‑1 (MQ=255)
   gACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATgcggc                                                                                                                                     <  1:263803/149‑1 (MQ=255)
        aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                <  3:1097956/149‑1 (MQ=255)
        aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                <  3:1565695/149‑1 (MQ=255)
        aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                <  3:1691981/149‑1 (MQ=255)
        aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                >  4:1691981/1‑149 (MQ=255)
        aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTAGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                >  4:1565695/1‑149 (MQ=255)
         ccGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAAc                                                                                                                               <  1:398233/149‑1 (MQ=255)
            aaTGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGcc                                                                                                                            >  2:445488/1‑149 (MQ=255)
             atggatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCt                                                                                                                           <  4:818897/149‑1 (MQ=255)
                gatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                        <  3:9017/149‑1 (MQ=255)
                    gTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                    >  2:719069/1‑149 (MQ=255)
                     ttACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                   >  2:1200504/1‑149 (MQ=255)
                     ttACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                   <  1:719069/149‑1 (MQ=255)
                           aaTCGCTGGCCGATTTCCTCGACCCACATGATCCAGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACGGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACaa                                                                                                             >  3:1259684/1‑149 (MQ=255)
                               gCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAcc                                                                                                         >  1:75315/1‑149 (MQ=255)
                                    ccGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCt                                                                                                    >  1:1306113/1‑149 (MQ=255)
                                      gATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTaa                                                                                                  >  1:385909/1‑149 (MQ=255)
                                       aTTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACg                                                                                                                              <  4:1824808/120‑1 (MQ=255)
                                       aTTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACg                                                                                                                              >  3:1824808/1‑120 (MQ=255)
                                       aTTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAt                                                                                                 <  1:1200504/149‑1 (MQ=255)
                                        tttCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                                                                >  4:192743/1‑149 (MQ=255)
                                        tttCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                                                                <  2:75315/149‑1 (MQ=255)
                                          tCCTCGACCCACATTATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                                                              <  2:1306113/149‑1 (MQ=255)
                                            cTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCt                                                                                                    >  4:602056/1‑141 (MQ=255)
                                            cTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAat                                                                                            <  4:1259684/149‑1 (MQ=255)
                                            cTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAat                                                                                            <  1:1556926/149‑1 (MQ=255)
                                            cTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                  <  1:57443/127‑2 (MQ=38)
                                            cTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                  >  2:57443/1‑126 (MQ=38)
                                               gACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatata                                                                                         <  1:1580448/149‑1 (MQ=255)
                                                aCCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatatat                                                                                        >  3:1146448/1‑149 (MQ=255)
                                                   cacaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGc                                                                                                     <  4:998273/133‑1 (MQ=255)
                                                   cacaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGc                                                                                                     >  3:998273/1‑133 (MQ=255)
                                                    acaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                <  3:1854297/105‑1 (MQ=255)
                                                    acaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGaa                                                                                                                                >  4:1854297/1‑105 (MQ=255)
                                                    acaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                                                                    >  2:1070248/1‑149 (MQ=255)
                                                      aTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                           >  3:1583401/1‑124 (MQ=255)
                                                      aTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                           <  4:1583401/124‑1 (MQ=255)
                                                      aTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                                  <  1:322492/149‑1 (MQ=255)
                                                        gATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGa                                                                                <  3:192743/149‑1 (MQ=255)
                                                        gATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGttc                                                                                                                                 <  1:172168/100‑3 (MQ=255)
                                                        gATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGttc                                                                                                                                 >  2:172168/1‑98 (MQ=255)
                                                         aTCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                                                               <  2:295913/149‑1 (MQ=255)
                                                          tCCGGGTAAAACGGTGGGAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACt                                                                              >  3:437546/1‑149 (MQ=255)
                                                           ccGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTa                                                                             <  4:1715342/149‑1 (MQ=255)
                                                            cGGGTAAAACGGTGGGAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                                                            <  4:437546/149‑1 (MQ=255)
                                                            cGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCaa                                                                                             >  1:151463/1‑132 (MQ=255)
                                                            cGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCaa                                                                                             <  2:151463/132‑1 (MQ=255)
                                                            cGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                                                            <  4:1146448/149‑1 (MQ=255)
                                                                tAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGTTCACGccgcatc                                                                                                                       >  2:79626/1‑97 (MQ=38)
                                                                tAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATAAGGCGTTCACGccgcatc                                                                                                                       <  1:79626/102‑6 (MQ=38)
                                                                   aaCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                                     >  3:1668726/1‑149 (MQ=255)
                                                                         ggAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATaa                                                               <  2:1802863/149‑1 (MQ=255)
                                                                          gAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGCCTTGCTaaa                                                                                                 >  1:1341370/1‑113 (MQ=255)
                                                                          gAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGCCTTGCTAAt                                                                                                 <  2:1341370/114‑1 (MQ=255)
                                                                          gAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCtt                                                                                                       >  4:1381082/1‑108 (MQ=255)
                                                                           aaGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                                     >  1:860411/1‑141 (MQ=255)
                                                                           aaGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                                     <  2:860411/141‑1 (MQ=255)
                                                                            aGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCTAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTg                                                                                                      <  1:239526/107‑1 (MQ=255)
                                                                            aGGTTATCCGGCACCTAAACGCGCGGTGCTGAATGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTg                                                                                                      >  2:239526/1‑107 (MQ=255)
                                                                             ggTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCTAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAACATATTGCAGGGACt                                                                              >  4:367182/1‑130 (MQ=255)
                                                                             ggTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCAt                                                           <  1:445488/149‑1 (MQ=255)
                                                                             ggTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAACATATTGCAGGGACt                                                                              <  3:367182/130‑1 (MQ=255)
                                                                              gTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCATAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATa                                                          <  4:1668726/149‑1 (MQ=255)
                                                                                     ggCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATc                                                   <  3:999317/149‑1 (MQ=255)
                                                                                         ccTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGc                                               <  2:385909/149‑1 (MQ=255)
                                                                                          cTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                                                          <  4:438317/105‑1 (MQ=255)
                                                                                          cTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                                                          >  3:438317/1‑105 (MQ=255)
                                                                                                         tGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                        >  3:1788101/1‑60 (MQ=255)
                                                                                                         tGATTGCGCGCAAGCCGTAAAGGACTGGTATTACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                        <  4:1788101/60‑1 (MQ=255)
                                                                                                            ttGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAg                                                                                   <  3:473580/94‑1 (MQ=255)
                                                                                                            ttGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAg                                                                                   >  4:473580/1‑94 (MQ=255)
                                                                                                             tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                                                              >  1:1570091/1‑82 (MQ=255)
                                                                                                             tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                                                              <  2:1570091/82‑1 (MQ=255)
                                                                                                             tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAat                                                                                            >  2:867948/1‑84 (MQ=255)
                                                                                                             tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAat                                                                                            <  1:867948/84‑1 (MQ=255)
                                                                                                                    aaGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatata                                                                                         >  3:1612358/1‑80 (MQ=255)
                                                                                                                    aaGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatata                                                                                         <  4:1612358/80‑1 (MQ=255)
                                                                                                                    aaGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                            >  1:13775/1‑141 (MQ=255)
                                                                                                                    aaGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                            <  2:13775/141‑1 (MQ=255)
                                                                                                                         gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCt                                                                                                        <  1:446450/60‑1 (MQ=255)
                                                                                                                         gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCt                                                                                                        >  2:446450/1‑60 (MQ=255)
                                                                                                                         gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                                  >  1:1379607/1‑114 (MQ=255)
                                                                                                                         gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                                  <  2:1379607/114‑1 (MQ=255)
                                                                                                                          tAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                   <  1:1453359/128‑1 (MQ=255)
                                                                                                                          tAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                   >  2:1453359/1‑128 (MQ=255)
                                                                                                                          tAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                   <  1:423960/128‑1 (MQ=255)
                                                                                                                          tAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                   >  2:423960/1‑128 (MQ=255)
                                                                                                                                tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                                                               >  4:1862561/1‑78 (MQ=255)
                                                                                                                                tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                                                               <  3:1862561/78‑1 (MQ=255)
                                                                                                                                  tGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                                                                    >  2:891692/1‑71 (MQ=255)
                                                                                                                                  tGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                                                                    <  1:891692/71‑1 (MQ=255)
                                                                                                                                    gtagtaCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTGAATATATTGCAgg                                                                                  >  4:994666/1‑71 (MQ=255)
                                                                                                                                    gTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                                  <  3:994666/71‑1 (MQ=255)
                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                                  <  2:1343609/94‑1 (MQ=255)
                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                                  >  1:1343609/1‑94 (MQ=255)
                                                                                                                                               ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGc                                            <  1:112245/98‑1 (MQ=255)
                                                                                                                                                 aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt  >  4:1090100/1‑138 (MQ=255)
                                                                                                                                                 aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt  <  3:1090100/138‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
AGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATT  >  NC_000913/1954289‑1954571

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