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.

ACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  NC_000913/1194507‑1194802
                                                                                                                                                   |                                                                                                                                                     
aCGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                                    >  2:1326385/1‑150 (MQ=255)
 cgcgGGTTCGTCAACCCGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTCAGCTTACGCCATTGCGCAAGGCTTGAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGt                                                                                                                                                   <  2:1032542/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACTCGGCTGGGt                                                                                                                                                   <  2:480782/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGt                                                                                                                                                   <  1:1412003/150‑1 (MQ=255)
 cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGGGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGt                                                                                                                                                   <  1:834942/150‑1 (MQ=255)
      gTTCGTCAACCAGCTCCTCGCCGGCGGGCAGGGTAGGGCCATCAGTTAAGGTTACGCCATTGCGCAATGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTATCCGGTGCGTTTACCCGGCTGGGTTAAAc                                                                                                                                              <  1:882026/150‑1 (MQ=255)
      gTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAc                                                                                                                                              <  2:15251/150‑1 (MQ=255)
        tCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAcgc                                                                                                                                            <  1:385102/150‑1 (MQ=255)
         cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAcgcg                                                                                                                                           >  2:203667/1‑150 (MQ=255)
         cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAcgcg                                                                                                                                           >  2:578251/1‑150 (MQ=255)
          gtaaaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGACATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAAcgcac                                                                                                                                          <  2:1064891/147‑3 (MQ=255)
          gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAGCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATGTTTCCGGTGCGTTTACCCGGCTGGGTTACACgcgc                                                                                                                                          >  1:1661489/1‑150 (MQ=255)
           tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGAGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGcc                                                                                                                                         <  2:1681557/150‑1 (MQ=255)
           gcAACCAGCACCGCGCCGGCGGGGAGTGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCCTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCCGCTGGGTTAAACGCGcc                                                                                                                                         <  2:290014/149‑1 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAc                                                                                                                                              >  2:742574/1‑144 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAc                                                                                                                                              <  1:742574/144‑1 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                        <  1:314061/150‑2 (MQ=255)
            cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                        <  1:1435613/150‑2 (MQ=255)
             aaacAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAACATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTCCACATAATAGATATTTCCGGTGCGTTTGCCCGGCTGGGTTAACCGCCCCTg                                                                                                                                       >  1:425419/4‑150 (MQ=255)
              aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAccgga                                                                                                                                      <  2:508917/150‑4 (MQ=255)
              aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                                      >  1:1362462/1‑150 (MQ=255)
                  gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                                    <  2:876744/132‑1 (MQ=255)
                     ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGTGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACgcgc                                                                                                                               <  2:1061891/150‑1 (MQ=255)
                     ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAAcgcac                                                                                                                               >  2:1632142/1‑148 (MQ=255)
                      cgcgCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              <  2:1465720/150‑1 (MQ=255)
                       gcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCg                                                                                                                             >  1:1708738/1‑150 (MQ=255)
                            ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACTCGCCTGCAACGCGCCGTTAtt                                                                                                                        >  2:1158366/1‑150 (MQ=255)
                            ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                        <  2:671620/150‑1 (MQ=255)
                              ggggCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGTCTTCTAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                      <  2:765997/149‑1 (MQ=255)
                              cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                      <  2:1362462/150‑1 (MQ=255)
                              cggGCAGGGGAGGGCGATGATTTAAGGTTACGCCTTTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCTCATTATAGATTTTTCCGGTGCGGTTACCCGGCTTGGTTAAACGCGACTGCAACGCGCCGTTATTgg                                                                                                                      <  2:711203/150‑1 (MQ=255)
                                 gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                   <  2:1047076/150‑1 (MQ=255)
                                   aGGGTAGGGCCAGCTTTTAAGGTTTCGCCATTGCGCAAGGCTGCACGTGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGTTTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGc                                                                                                                 <  1:27367/150‑1 (MQ=255)
                                    gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                                >  1:1437420/1‑150 (MQ=255)
                                     ggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGAGCGCCTTGTGTGGGAATACCTTCCACCTGCACATAATACATTTTTCCGGTGCGATTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGcgc                                                                                                               >  2:715241/1‑148 (MQ=255)
                                      gTATGGCCATCATTTAAGGTTACGCCATTGCGCAATGCTTCAAGTGCGTCTTGTGTTGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAcc                                                                                                              <  2:847029/150‑1 (MQ=255)
                                        aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              >  2:1243877/1‑132 (MQ=255)
                                        aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                              <  1:1243877/132‑1 (MQ=255)
                                         gggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGc                                                                                                           <  2:730359/150‑1 (MQ=255)
                                          ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                          >  2:898097/1‑150 (MQ=255)
                                           gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                         >  1:40064/1‑150 (MQ=255)
                                               tcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGAAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCt                                                                                                     >  1:798279/1‑150 (MQ=255)
                                                        gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            >  1:1721606/1‑150 (MQ=255)
                                                        gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            <  2:1155512/150‑1 (MQ=255)
                                                         ttACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcgc                                                                                           <  1:249695/150‑1 (MQ=255)
                                                            cGCCATAGCGCAAGGCTTCAAGTGCGTCTTGTGTTGGAAAACCTTCCACCTGCACATAATAGATATTACCGGTGCTTTTACCCTGATAGGTTAAACGCGCATGCAACGCTCCGTTATATGACAGAAACAGCAACCCATCGCTATCGCGGt                                                                                        >  1:98198/1‑150 (MQ=255)
                                                              ccATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATCCCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCa                                                                                      >  2:1261038/1‑150 (MQ=255)
                                                                    cgcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCgc                                                                                >  1:689797/1‑149 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           <  1:894068/150‑1 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           <  1:954672/150‑1 (MQ=255)
                                                                         ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                           >  1:1119019/1‑150 (MQ=255)
                                                                                     gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                               >  1:1550485/1‑150 (MQ=255)
                                                                                     gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                               >  1:1165992/1‑150 (MQ=255)
                                                                                       cTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCt                                                             >  1:1432340/1‑150 (MQ=255)
                                                                                        ttGTGTGGTAATACCTTGCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTTCTGCATAAACACCCTg                                                            <  2:311882/150‑1 (MQ=255)
                                                                                         tgtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGa                                                           >  1:1203184/1‑150 (MQ=255)
                                                                                           ggtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAc                                                         <  2:595928/149‑1 (MQ=255)
                                                                                            gtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTTAAcc                                                        >  2:758421/1‑150 (MQ=255)
                                                                                            gtgGGAATACCTTCCACCTGCACATAATAGAATTTTCCGGTGCGTTTACCCGGCAGGGTTAAACTCGCCTGCAACGCGAAGATATAGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAcc                                                        <  1:850026/150‑1 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     <  2:514679/150‑1 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     >  1:1141116/1‑150 (MQ=255)
                                                                                               ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGGGCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     <  2:206437/150‑1 (MQ=255)
                                                                                                    aCCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCTCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGaa                                                >  2:84717/1‑150 (MQ=255)
                                                                                                     ccTTCCACCTGCACATAATAGATTATTCCGGTGCGGGGACCCGGCTGGGCTAAACGCGCCTGCAACGCGCCGTTATGGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAt                                               <  1:251945/150‑1 (MQ=255)
                                                                                                      ctaccACCTGCACATAATACATTTTTCCGGTGCGGTTACCCGGCTGGGGTAAACGCGCCTGCACCGCGACGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                              <  2:338296/147‑1 (MQ=255)
                                                                                                      cTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                              >  1:1022349/1‑150 (MQ=255)
                                                                                                            ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  2:1662562/1‑150 (MQ=255)
                                                                                                            ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  2:1410848/1‑150 (MQ=255)
                                                                                                             cTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                                       <  2:1500460/150‑1 (MQ=255)
                                                                                                               gCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATg                                     >  1:1549775/1‑150 (MQ=255)
                                                                                                                cacaTAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                                    >  1:332884/1‑150 (MQ=255)
                                                                                                                cacaTAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                                    >  1:415250/1‑150 (MQ=255)
                                                                                                                 acaTAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCGGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTCTCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATTAATTCTTTTAATGtt                                   <  2:141952/150‑1 (MQ=255)
                                                                                                                   ataataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                                 <  2:798279/150‑1 (MQ=255)
                                                                                                                    taataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTg                                >  1:249204/1‑150 (MQ=255)
                                                                                                                     aataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGc                               >  1:1416647/1‑150 (MQ=255)
                                                                                                                       taGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                             >  2:1427192/1‑150 (MQ=255)
                                                                                                                        aGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCCGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            >  1:705867/1‑150 (MQ=255)
                                                                                                                        aGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            >  1:157329/1‑150 (MQ=255)
                                                                                                                          aTTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                          <  1:874758/150‑1 (MQ=255)
                                                                                                                           tttttCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCACCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                         <  1:1632142/150‑1 (MQ=255)
                                                                                                                                 cGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAACCCGCCTGAACCGGGATGCATTCTTTTAATGTTTTGCGTCAGGCTTCag                   >  1:820719/1‑149 (MQ=255)
                                                                                                                                 cGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCAt                   >  1:1403728/1‑150 (MQ=255)
                                                                                                                                  ggTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc                  >  1:496960/1‑150 (MQ=255)
                                                                                                                                  ggTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc                  <  1:328800/150‑1 (MQ=255)
                                                                                                                                  ggTGCGATTACCCGGCTCGGTTAAACGAGCCTGCAACGCGCCGTTATTGGTCAGCACCACCAACCCTACGCAAACGCGGTCAAGGCGACCAGCTGCATAAACACCCAGAACCCGGATGAATTCTTTTAATGTATTGCGTCCGGCTTCATc                  >  2:48865/1‑150 (MQ=255)
                                                                                                                                    tGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGGCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTGCATCgg                <  2:221721/150‑1 (MQ=255)
                                                                                                                                    ggCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATACACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGGCCGGCTTCATCgg                <  2:1661489/149‑1 (MQ=255)
                                                                                                                                    cgCGTTTAACCCAGCCGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCCCCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACGGGGATGAATTCTTTTAATGTTTTGCGTCCGGCATCATCgg                >  2:120085/2‑150 (MQ=255)
                                                                                                                                      cGTTTAACCCAGCCGGGTAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg              >  1:1122566/1‑150 (MQ=255)
                                                                                                                                      cGTTTACCCGGCTGGGTTAAACGCGCTTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg              <  2:1165992/150‑1 (MQ=255)
                                                                                                                                       gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  1:1312777/1‑150 (MQ=255)
                                                                                                                                       gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  2:1157773/1‑150 (MQ=255)
                                                                                                                                       gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa             >  1:1571986/1‑150 (MQ=255)
                                                                                                                                         ttACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAc           >  2:1500172/1‑150 (MQ=255)
                                                                                                                                           aaCCCAGCCGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg         >  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)
                                                                                                                                                   |                                                                                                                                                     
ACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  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: