Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,577,803 T→C T514A (ACC→GCC)  yddA ← putative multidrug ABC transporter permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,577,8030TC99.3% 527.0 / NA 139T514A (ACC→GCC) yddAputative multidrug ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (71/67);  total (71/68)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 139 total aligned reads displayed.

TTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCA  >  NC_000913/1577656‑1577949
                                                                                                                                                   |                                                                                                                                                  
 ttATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGtt                                                                                                                                                >  2:1426203/1‑151 (MQ=255)
   ataAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCTGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTc                                                                                                                                              <  2:1519462/151‑1 (MQ=255)
   ataAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGttt                                                                                                                                               >  3:59894/1‑150 (MQ=255)
   ataAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTc                                                                                                                                              >  2:1021163/1‑151 (MQ=255)
   ataAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTc                                                                                                                                              <  4:480360/151‑1 (MQ=255)
      aaaCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTc                                                                                                                                           >  2:1431527/1‑151 (MQ=255)
       aaCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTc                                                                                                                                           >  1:365517/1‑150 (MQ=255)
       aaCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTc                                                                                                                                           >  2:1354021/1‑150 (MQ=255)
       aaCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTc                                                                                                                                           >  2:271760/1‑150 (MQ=255)
        aCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCaa                                                                                                                                         <  4:101111/151‑1 (MQ=255)
         ccGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAg                                                                                                                                        <  1:1303439/151‑1 (MQ=255)
          cgcgCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCa                                                                                                                                          >  4:589784/1‑148 (MQ=255)
          cgcgCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCa                                                                                                                                          >  3:544676/1‑148 (MQ=255)
          cgcgCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGa                                                                                                                                       <  1:258098/151‑1 (MQ=255)
           gcgcTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAg                                                                                                                                        >  2:1162416/1‑149 (MQ=255)
           gcgcTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGaa                                                                                                                                      <  2:171360/151‑1 (MQ=255)
                aaTGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATc                                                                                                                                 >  3:169473/1‑151 (MQ=255)
                   gTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCa                                                                                                                               <  2:209295/150‑1 (MQ=255)
                      acaAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTc                                                                                                                           >  4:1072/1‑151 (MQ=255)
                        aaaTATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCgg                                                                                                                         <  2:1688652/151‑1 (MQ=255)
                        aaaTATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCgg                                                                                                                         >  3:191121/1‑151 (MQ=255)
                         aaTATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCgg                                                                                                                         >  2:1513029/1‑150 (MQ=255)
                            atcatcGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGAcg                                                                                                                      <  4:59893/150‑1 (MQ=255)
                                 cGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTaa                                                                                                                <  2:863920/151‑1 (MQ=255)
                                  ggCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTaaa                                                                                                               <  2:1456743/151‑1 (MQ=255)
                                  ggCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTaa                                                                                                                >  2:1674152/1‑150 (MQ=255)
                                    ccAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAAt                                                                                                             <  1:1279751/151‑1 (MQ=255)
                                     cAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAAtt                                                                                                            <  1:1458450/151‑1 (MQ=255)
                                       ggTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTaa                                                                                                          <  2:188172/151‑1 (MQ=255)
                                       ggTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTaa                                                                                                          <  1:1192441/151‑1 (MQ=255)
                                        gTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAAt                                                                                                         <  4:265341/151‑1 (MQ=255)
                                          tCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCg                                                                                                       >  2:1062536/1‑151 (MQ=255)
                                           ccAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCg                                                                                                       <  1:837078/150‑1 (MQ=255)
                                           ccAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGt                                                                                                      <  1:920236/151‑1 (MQ=255)
                                           ccAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGt                                                                                                      <  2:1152874/151‑1 (MQ=255)
                                            cAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCg                                                                                                       <  2:753350/149‑1 (MQ=255)
                                            cAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTg                                                                                                     >  2:1654374/1‑151 (MQ=255)
                                             aGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGc                                                                                                    >  4:141690/1‑151 (MQ=255)
                                                 acCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAgg                                                                                                <  4:369868/151‑1 (MQ=255)
                                                  ccGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAggg                                                                                               <  4:544656/151‑1 (MQ=255)
                                                   cGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGc                                                                                              >  3:358141/1‑151 (MQ=255)
                                                         gATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACg                                                                                        >  1:1388547/1‑151 (MQ=255)
                                                           tgtgTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACgtt                                                                                      >  1:1107528/1‑151 (MQ=255)
                                                           tgtgTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACgtt                                                                                      >  2:988368/1‑151 (MQ=255)
                                                               tAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGttt                                                                                  >  2:1237919/1‑151 (MQ=255)
                                                                  ccATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTc                                                                               <  1:1672376/151‑1 (MQ=255)
                                                                  ccATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTc                                                                               <  2:553861/151‑1 (MQ=255)
                                                                  ccATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTc                                                                               >  1:725527/1‑151 (MQ=255)
                                                                    aTAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTcgc                                                                             >  2:378293/1‑151 (MQ=255)
                                                                      aaTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTcgccg                                                                           >  2:679941/1‑151 (MQ=255)
                                                                      aaTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTcgccg                                                                           >  2:1035130/1‑151 (MQ=255)
                                                                       aTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTcgccgc                                                                          >  2:771376/1‑151 (MQ=255)
                                                                         gACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCt                                                                         <  3:155804/150‑1 (MQ=255)
                                                                           cgccgcTTGTTGGTACTTTTTCACGCACTAAGCGCTGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAAGTATTTTCTACCAGAACTATCCATTTCGGTCGTCGTAAAAGTAATCGTGCCAGGGCGATACGTTGTTTTTCGACGCTgg                                                                       <  1:916257/150‑1 (MQ=255)
                                                                               gcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGa                                                                  <  1:113088/151‑1 (MQ=255)
                                                                                 ttGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAat                                                                >  2:635029/1‑151 (MQ=255)
                                                                                  tgtgGGTTGTTTTTCACGCACTAAACGCAGTAAGCGGATAGACTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAATAAATATCCATTTCGGACGTCGTAAAAGTAATCGTGCCAGAGCGGTACGTTGTTTTTCGCCGCTGCAAAgt                                                                  <  1:1015478/148‑2 (MQ=255)
                                                                                  tgtgGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGaa                                                                 <  1:1648784/149‑1 (MQ=255)
                                                                                  tgtgGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAat                                                                >  2:1563905/1‑150 (MQ=255)
                                                                                   gtgGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAata                                                               <  4:407162/150‑1 (MQ=255)
                                                                                     gggTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATAtct                                                            <  2:814388/151‑1 (MQ=255)
                                                                                        tAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTccc                                                         >  3:310824/1‑151 (MQ=255)
                                                                                         aGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTcccc                                                        >  2:387802/1‑151 (MQ=255)
                                                                                           tttttCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTcccc                                                        >  1:1326474/1‑149 (MQ=255)
                                                                                             tttCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGc                                                     <  2:824004/150‑1 (MQ=255)
                                                                                              ttCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCg                                                    <  3:112061/150‑1 (MQ=255)
                                                                                              ttCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCg                                                    <  2:904851/150‑1 (MQ=255)
                                                                                               tCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGa                                                   <  1:1511652/150‑1 (MQ=255)
                                                                                                cACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATc                                                 <  2:52534/151‑1 (MQ=255)
                                                                                                 aCGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATc                                                 >  2:995935/1‑150 (MQ=255)
                                                                                                  cGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCg                                                >  1:1002913/1‑150 (MQ=255)
                                                                                                    cACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTgg                                             >  1:1470135/1‑151 (MQ=255)
                                                                                                        aaaCGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTc                                           <  1:556100/149‑1 (MQ=255)
                                                                                                        aaaCGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCa                                          >  2:1486277/1‑150 (MQ=255)
                                                                                                        aaaCGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCAt                                         >  2:511013/1‑151 (MQ=255)
                                                                                                          aCGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCAt                                         >  1:755354/1‑149 (MQ=255)
                                                                                                           cGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGaa                                      <  1:249792/151‑1 (MQ=255)
                                                                                                            gCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAAt                                     <  2:1630884/151‑1 (MQ=255)
                                                                                                              aGTAAGCGGATAGCCTCTTGTTCCTAAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATAc                                   >  1:1232810/1‑151 (MQ=255)
                                                                                                                 aaGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATAcgcg                                >  2:1256381/1‑151 (MQ=255)
                                                                                                                  aGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACgcgc                               >  2:516520/1‑151 (MQ=255)
                                                                                                                  aGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACgcgc                               >  2:555854/1‑151 (MQ=255)
                                                                                                                  aGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACgcgc                               <  2:1296612/151‑1 (MQ=255)
                                                                                                                   gCGGATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCa                              <  1:1353873/151‑1 (MQ=255)
                                                                                                                     ggATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGc                            >  3:365285/1‑151 (MQ=255)
                                                                                                                      gATAGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGcc                           <  2:879383/151‑1 (MQ=255)
                                                                                                                         aGCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAAt                        <  2:615870/151‑1 (MQ=255)
                                                                                                                          gCCTCTTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAAtt                       >  2:321116/1‑151 (MQ=255)
                                                                                                                            ctctTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTc                     >  3:451272/1‑151 (MQ=255)
                                                                                                                             tctTGTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTcc                    >  2:573840/1‑151 (MQ=255)
                                                                                                                                 gTTCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCaa                 <  1:1426047/150‑1 (MQ=255)
                                                                                                                                   tCCTCAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGAc              >  2:300129/1‑151 (MQ=255)
                                                                                                                                    ccTCAAGATGAGAGGCAGTTTCGTCAAGAAATCTCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGAcc             >  3:306873/1‑151 (MQ=255)
                                                                                                                                       cAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAAc          >  2:664156/1‑151 (MQ=255)
                                                                                                                                       cAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAAc          <  1:555772/151‑1 (MQ=255)
                                                                                                                                       cAAGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTCTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAAc          <  2:436120/151‑1 (MQ=255)
                                                                                                                                         aGATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCt        <  1:271718/151‑1 (MQ=255)
                                                                                                                                          gATGAGAGGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTg       <  1:747747/151‑1 (MQ=255)
                                                                                                                                              agagGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGc   >  1:567869/1‑151 (MQ=255)
                                                                                                                                               gagGCAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCa  <  1:103806/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCA  >  NC_000913/1577656‑1577949

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑