Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,653 T→C 8.2% Q100R (CAG→CGG)  rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,194,6530TC 8.2% 340.7 / 17.0 123Q100R (CAG→CGG) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base T (58/54);  new base C (4/6);  total (62/61)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-02

Only 100 of 125 total aligned reads displayed.

ACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  NC_000913/1194507‑1194802
                                                                                                                                                   |                                                                                                                                                     
aCGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTggg                                                                                                                                                    >  2:1326385/1‑150 (MQ=255)
 cgcgGGTTCGTCAACCCGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTCAGCTTACGCCATTGCGCAAGGCTTGAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGt                                                                                                                                                   <  2:1032542/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACTCGGCTGGGt                                                                                                                                                   <  2:480782/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGt                                                                                                                                                   <  1:1412003/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGGGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGt                                                                                                                                                   <  1:834942/150‑1 (MQ=255)
      gTTCGTCAACCAGCTCCTCGCCGGCGGGCAGGGTAGGGCCATCAGTTAAGGTTACGCCATTGCGCAATGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTATCCGGTGCGTTT‑ACCCGGCTGGGTTAAAc                                                                                                                                              <  1:882026/150‑1 (MQ=255)
      gTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAc                                                                                                                                              <  2:15251/150‑1 (MQ=255)
        tCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgc                                                                                                                                            <  1:385102/150‑1 (MQ=255)
         cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgcg                                                                                                                                           >  2:203667/1‑150 (MQ=255)
         cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgcg                                                                                                                                           >  2:578251/1‑150 (MQ=255)
          gtaaaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGACATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAAcgcac                                                                                                                                          <  2:1064891/147‑3 (MQ=255)
          gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAGCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATGTTTCCGGTGCGTTT‑ACCCGGCTGGGTTACACgcgc                                                                                                                                          >  1:1661489/1‑150 (MQ=255)
           tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGAGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                         <  2:1681557/150‑1 (MQ=255)
           gcAACCAGCACCGCGCCGGCGGGGAGTGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCCTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCCGCTGGGTTAAACGCGcc                                                                                                                                         <  2:290014/149‑1 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAc                                                                                                                                              >  2:742574/1‑144 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAc                                                                                                                                              <  1:742574/144‑1 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAccg                                                                                                                                        <  1:314061/150‑2 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAccg                                                                                                                                        <  1:1435613/150‑2 (MQ=255)
             aaacAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAACATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTCCACATAATAGATATTTCCGGTGCGTTT‑GCCCGGCTGGGTTAACCGCCCCTg                                                                                                                                       >  1:425419/4‑150 (MQ=255)
              aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAACGCAccgga                                                                                                                                      <  2:508917/150‑4 (MQ=255)
              aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                                      >  1:1362462/1‑150 (MQ=255)
                  gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTggg                                                                                                                                                    <  2:876744/132‑1 (MQ=255)
                     ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGTGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACgcgc                                                                                                                               <  2:1061891/150‑1 (MQ=255)
                     ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcac                                                                                                                               >  2:1632142/1‑148 (MQ=255)
                      cgcgCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              <  2:1465720/150‑1 (MQ=255)
                       gcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCg                                                                                                                             >  1:1708738/1‑150 (MQ=255)
                            ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACTCGCCTGCAACGCGCCGTTAtt                                                                                                                        >  2:1158366/1‑150 (MQ=255)
                            ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                        <  2:671620/150‑1 (MQ=255)
                              ggggCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGTCTTCTAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                      <  2:765997/149‑1 (MQ=255)
                              cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                      <  2:1362462/150‑1 (MQ=255)
                              cggGCAGGGGAGGGCGATGATTTAAGGTTACGCCTTTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCTCATTATAGATTTTTCCGGTGCGGTT‑ACCCGGCTTGGTTAAACGCGACTGCAACGCGCCGTTATTgg                                                                                                                      <  2:711203/150‑1 (MQ=255)
                                 gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                   <  2:1047076/150‑1 (MQ=255)
                                   aGGGTAGGGCCAGCTTTTAAGGTTTCGCCATTGCGCAAGGCTGCACGTGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGTTTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGc                                                                                                                 <  1:27367/150‑1 (MQ=255)
                                    gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                                >  1:1437420/1‑150 (MQ=255)
                                     ggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGAGCGCCTTGTGTGGGAATACCTTCCACCTGCACATAATACATTTTTCCGGTGCGATT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGcgc                                                                                                               >  2:715241/1‑148 (MQ=255)
                                      gTATGGCCATCATTTAAGGTTACGCCATTGCGCAATGCTTCAAGTGCGTCTTGTGTTGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAcc                                                                                                              <  2:847029/150‑1 (MQ=255)
                                        aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              >  2:1243877/1‑132 (MQ=255)
                                        aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              <  1:1243877/132‑1 (MQ=255)
                                         gggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGc                                                                                                           <  2:730359/150‑1 (MQ=255)
                                          ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                          >  2:898097/1‑150 (MQ=255)
                                           gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                         >  1:40064/1‑150 (MQ=255)
                                               tcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGAAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCt                                                                                                     >  1:798279/1‑150 (MQ=255)
                                                        gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            >  1:1721606/1‑150 (MQ=255)
                                                        gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            <  2:1155512/150‑1 (MQ=255)
                                                         ttACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcgc                                                                                           <  1:249695/150‑1 (MQ=255)
                                                            cGCCATAGCGCAAGGCTTCAAGTGCGTCTTGTGTTGGAAAACCTTCCACCTGCACATAATAGATATTACCGGTGCTTTT‑ACCCTGATAGGTTAAACGCGCATGCAACGCTCCGTTATATGACAGAAACAGCAACCCATCGCTATCGCGGt                                                                                        >  1:98198/1‑150 (MQ=255)
                                                              ccATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATCCCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCa                                                                                      >  2:1261038/1‑150 (MQ=255)
                                                                    cgcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCgc                                                                                >  1:689797/1‑149 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           <  1:894068/150‑1 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           <  1:954672/150‑1 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           >  1:1119019/1‑150 (MQ=255)
                                                                                     gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                               >  1:1550485/1‑150 (MQ=255)
                                                                                     gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                               >  1:1165992/1‑150 (MQ=255)
                                                                                       cTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCt                                                             >  1:1432340/1‑150 (MQ=255)
                                                                                        ttGTGTGGTAATACCTTGCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTTCTGCATAAACACCCTg                                                            <  2:311882/150‑1 (MQ=255)
                                                                                         tgtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGa                                                           >  1:1203184/1‑150 (MQ=255)
                                                                                           ggtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAc                                                         <  2:595928/149‑1 (MQ=255)
                                                                                            gtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTTAAcc                                                        >  2:758421/1‑150 (MQ=255)
                                                                                            gtgGGAATACCTTCCACCTGCACATAATAGAATTTTCCGGTGCGTTT‑ACCCGGCAGGGTTAAACTCGCCTGCAACGCGAAGATATAGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAcc                                                        <  1:850026/150‑1 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     <  2:514679/150‑1 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     >  1:1141116/1‑150 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGGGCGGTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     <  2:206437/150‑1 (MQ=255)
                                                                                                    aCCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCTCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGaa                                                >  2:84717/1‑150 (MQ=255)
                                                                                                     ccTTCCACCTGCACATAATAGATTATTCCGGTGCGGGG‑ACCCGGCTGGGCTAAACGCGCCTGCAACGCGCCGTTATGGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAt                                               <  1:251945/150‑1 (MQ=255)
                                                                                                      ctaccACCTGCACATAATACATTTTTCCGGTGCGGTT‑ACCCGGCTGGGGTAAACGCGCCTGCACCGCGACGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                              <  2:338296/147‑1 (MQ=255)
                                                                                                      cTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                              >  1:1022349/1‑150 (MQ=255)
                                                                                                            ccTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  2:1662562/1‑150 (MQ=255)
                                                                                                            ccTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  2:1410848/1‑150 (MQ=255)
                                                                                                             cTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                                       <  2:1500460/150‑1 (MQ=255)
                                                                                                               gCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATg                                     >  1:1549775/1‑150 (MQ=255)
                                                                                                                cacaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                                    >  1:332884/1‑150 (MQ=255)
                                                                                                                cacaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                                    >  1:415250/1‑150 (MQ=255)
                                                                                                                 acaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCGGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTCTCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATTAATTCTTTTAATGtt                                   <  2:141952/150‑1 (MQ=255)
                                                                                                                   ataataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                                 <  2:798279/150‑1 (MQ=255)
                                                                                                                    taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTg                                >  1:249204/1‑150 (MQ=255)
                                                                                                                     aataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGc                               >  1:1416647/1‑150 (MQ=255)
                                                                                                                       taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                             >  2:1427192/1‑150 (MQ=255)
                                                                                                                        aGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCCGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            >  1:705867/1‑150 (MQ=255)
                                                                                                                        aGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            >  1:157329/1‑150 (MQ=255)
                                                                                                                          aTTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                          <  1:874758/150‑1 (MQ=255)
                                                                                                                           tttttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCACCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                         <  1:1632142/150‑1 (MQ=255)
                                                                                                                                 cGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAACCCGCCTGAACCGGGATGCATTCTTTTAATGTTTTGCGTCAGGCTTCag                   >  1:820719/1‑149 (MQ=255)
                                                                                                                                 cGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCAt                   >  1:1403728/1‑150 (MQ=255)
                                                                                                                                  ggTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc                  >  1:496960/1‑150 (MQ=255)
                                                                                                                                  ggTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc                  <  1:328800/150‑1 (MQ=255)
                                                                                                                                  ggTGCGATT‑ACCCGGCTCGGTTAAACGAGCCTGCAACGCGCCGTTATTGGTCAGCACCACCAACCCTACGCAAACGCGGTCAAGGCGACCAGCTGCATAAACACCCAGAACCCGGATGAATTCTTTTAATGTATTGCGTCCGGCTTCATc                  >  2:48865/1‑150 (MQ=255)
                                                                                                                                    tGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGGCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTGCATCgg                <  2:221721/150‑1 (MQ=255)
                                                                                                                                    ggCGGTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATACACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGGCCGGCTTCATCgg                <  2:1661489/149‑1 (MQ=255)
                                                                                                                                    cgCGTTTAACCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCCCCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACGGGGATGAATTCTTTTAATGTTTTGCGTCCGGCATCATCgg                >  2:120085/2‑150 (MQ=255)
                                                                                                                                      cGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg              >  1:1122566/1‑150 (MQ=255)
                                                                                                                                      cGTTT‑ACCCGGCTGGGTTAAACGCGCTTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg              <  2:1165992/150‑1 (MQ=255)
                                                                                                                                       gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  1:1312777/1‑150 (MQ=255)
                                                                                                                                       gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  2:1157773/1‑150 (MQ=255)
                                                                                                                                       gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  1:1571986/1‑150 (MQ=255)
                                                                                                                                         tt‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAc           >  2:1500172/1‑150 (MQ=255)
                                                                                                                                           aaCCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg         >  2:352680/2‑150 (MQ=255)
                                                                                                                                            aCCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg         <  2:40064/150‑1 (MQ=255)
                                                                                                                                             cccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCGTGAGCCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc        >  1:1422695/1‑150 (MQ=255)
                                                                                                                                                  cTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCCGTCAAGGCGACCTGCAGCATAAACACCCTGACCCGGGATGAATGCTTTTAATGTTTTGCGTCCGGATTCATCGGTCAACTGCGgaag   >  2:335555/1‑150 (MQ=255)
                                                                                                                                                   tGGGTTAAACGAGCCTGCAACGCGCCGTTATTTGTCAGCACCATCAACCCTTCGCTATCGCTGTCAAGTCGACCTGCTGCATAAACACCCTGAACCGGTATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTTAACTGCGgaaga  <  2:469449/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  NC_000913/1194507‑1194802

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