Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,645:1 +A 6.7% coding (307/654 nt) 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,6451.A 6.7% 461.0 / 26.0 120coding (307/654 nt)rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base . (57/55);  new base A (2/6);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

CCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  NC_000913/1194499‑1194794
                                                                                                                                                   |                                                                                                                                                     
 cATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑Accc                                                                                                                                                   <  2:1084210/150‑1 (MQ=255)
  aTAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCg                                                                                                                                                  <  2:918763/150‑1 (MQ=255)
    aaCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGc                                                                                                                                                >  2:606717/1‑150 (MQ=255)
    aaCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGc                                                                                                                                                >  2:1308310/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)
         cgcgGGTTCGGCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACTCGGCTGGGt                                                                                                                                           <  2:375696/150‑1 (MQ=255)
          gcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGGAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGtt                                                                                                                                          <  1:616114/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:578251/1‑150 (MQ=255)
                 cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgcg                                                                                                                                   >  2:203667/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)
                     aaacAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAACATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTCCACATAATAGATATTTCCGGTGCGTTT‑GCCCGGCTGGGTTAACCGCCCCTg                                                                                                                               >  1:425419/4‑150 (MQ=255)
                     aaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAccgg                                                                                                                               <  1:810989/150‑3 (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)
                          gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTggg                                                                                                                                            >  1:876744/1‑132 (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)
                                                aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      <  1:1243877/132‑1 (MQ=255)
                                                aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      >  2:1243877/1‑132 (MQ=255)
                                                 gggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGc                                                                                                   <  2:730359/150‑1 (MQ=255)
                                                 gggCCATCATTTAAGGTTACGCCATTGAGCAAGGCTTCAAGTGCGTATTGTGTGGGAATACCTTCCACCTTCACATAATATATTTTTCCGGTGCGTTT‑ACACGGCTGGGTTAAACGCGACTGCAACGCGCCGTTATTGGTCAGCACAAGc                                                                                                   >  2:117000/1‑150 (MQ=255)
                                                  ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                  >  1:1448968/1‑150 (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)
                                                            tAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAACCGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCCACCCTTCGCt                                                                                        >  1:1629426/1‑150 (MQ=255)
                                                            tAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCt                                                                                        >  1:197903/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)
                                                                                 ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                   >  1:1119019/1‑150 (MQ=255)
                                                                                 ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                   <  1:954672/150‑1 (MQ=255)
                                                                                             gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                       >  1:1550485/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                                             >  1:1141116/1‑150 (MQ=255)
                                                                                                       ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                             <  2:514679/150‑1 (MQ=255)
                                                                                                       ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGGGCGGTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                             <  2:206437/150‑1 (MQ=255)
                                                                                                          aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATg                                          <  2:1360/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)
                                                                                                                  cACCTGCTCATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                  <  2:1448968/150‑1 (MQ=255)
                                                                                                                  cACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                  <  1:1158366/150‑1 (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)
                                                                                                                     cTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                               <  2:197903/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)
                                                                                                                                aGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCCGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                    >  1:705867/1‑150 (MQ=255)
                                                                                                                                  aTTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                  <  1:874758/150‑1 (MQ=255)
                                                                                                                                   tttttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCACCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                 <  1:1632142/150‑1 (MQ=255)
                                                                                                                                    ttttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGc                <  2:836592/150‑1 (MQ=255)
                                                                                                                                          ggTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc          >  1:496960/1‑150 (MQ=255)
                                                                                                                                          ggTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc          <  1:328800/150‑1 (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)
                                                                                                                                             gCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGt       <  1:1326385/150‑1 (MQ=255)
                                                                                                                                              cGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg      >  1:1122566/1‑150 (MQ=255)
                                                                                                                                               gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  2:1157773/1‑150 (MQ=255)
                                                                                                                                               gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:1312777/1‑150 (MQ=255)
                                                                                                                                               gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:1571986/1‑150 (MQ=255)
                                                                                                                                                 tt‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAc   >  2:1500172/1‑150 (MQ=255)
                                                                                                                                                  t‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt  >  1:1419802/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  NC_000913/1194499‑1194794

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: ‑