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.

CCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  NC_000913/1194499‑1194794
                                                                                                                                                   |                                                                                                                                                     
 cATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAccc                                                                                                                                                   <  2:1084210/150‑1 (MQ=255)
  aTAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCg                                                                                                                                                  <  2:918763/150‑1 (MQ=255)
    aaCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGc                                                                                                                                                >  2:606717/1‑150 (MQ=255)
    aaCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGc                                                                                                                                                >  2:1308310/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)
         cgcgGGTTCGGCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACTCGGCTGGGt                                                                                                                                           <  2:375696/150‑1 (MQ=255)
          gcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGGAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGtt                                                                                                                                          <  1:616114/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:578251/1‑150 (MQ=255)
                 cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAcgcg                                                                                                                                   >  2:203667/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)
                     aaacAGCTCCGCGCCGGCGGGCAGGGTAGGGCCAACATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTCCACATAATAGATATTTCCGGTGCGTTTGCCCGGCTGGGTTAACCGCCCCTg                                                                                                                               >  1:425419/4‑150 (MQ=255)
                     aaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccgg                                                                                                                               <  1:810989/150‑3 (MQ=255)
                      aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAccgga                                                                                                                              <  2:508917/150‑4 (MQ=255)
                      aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                              >  1:1362462/1‑150 (MQ=255)
                          gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                            <  2:876744/132‑1 (MQ=255)
                          gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                            >  1:876744/1‑132 (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)
                                                aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      <  1:1243877/132‑1 (MQ=255)
                                                aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      >  2:1243877/1‑132 (MQ=255)
                                                 gggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGc                                                                                                   <  2:730359/150‑1 (MQ=255)
                                                 gggCCATCATTTAAGGTTACGCCATTGAGCAAGGCTTCAAGTGCGTATTGTGTGGGAATACCTTCCACCTTCACATAATATATTTTTCCGGTGCGTTTACACGGCTGGGTTAAACGCGACTGCAACGCGCCGTTATTGGTCAGCACAAGc                                                                                                   >  2:117000/1‑150 (MQ=255)
                                                  ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                  >  1:1448968/1‑150 (MQ=255)
                                                  ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                  >  2:898097/1‑150 (MQ=255)
                                                   gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                 >  1:40064/1‑150 (MQ=255)
                                                       tcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGAAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCt                                                                                             >  1:798279/1‑150 (MQ=255)
                                                            tAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAACCGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCCACCCTTCGCt                                                                                        >  1:1629426/1‑150 (MQ=255)
                                                            tAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCt                                                                                        >  1:197903/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)
                                                                                 ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                   >  1:1119019/1‑150 (MQ=255)
                                                                                 ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                   <  1:954672/150‑1 (MQ=255)
                                                                                             gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                       >  1:1550485/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                                             >  1:1141116/1‑150 (MQ=255)
                                                                                                       ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                             <  2:514679/150‑1 (MQ=255)
                                                                                                       ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGGGCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                             <  2:206437/150‑1 (MQ=255)
                                                                                                          aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATg                                          <  2:1360/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)
                                                                                                                  cACCTGCTCATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                  <  2:1448968/150‑1 (MQ=255)
                                                                                                                  cACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                  <  1:1158366/150‑1 (MQ=255)
                                                                                                                    ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                >  2:1662562/1‑150 (MQ=255)
                                                                                                                    ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                >  2:1410848/1‑150 (MQ=255)
                                                                                                                     cTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                               <  2:1500460/150‑1 (MQ=255)
                                                                                                                     cTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                               <  2:197903/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)
                                                                                                                                aGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCCGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                    >  1:705867/1‑150 (MQ=255)
                                                                                                                                  aTTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                  <  1:874758/150‑1 (MQ=255)
                                                                                                                                   tttttCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCACCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                 <  1:1632142/150‑1 (MQ=255)
                                                                                                                                    ttttCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGc                <  2:836592/150‑1 (MQ=255)
                                                                                                                                          ggTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc          >  1:496960/1‑150 (MQ=255)
                                                                                                                                          ggTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc          <  1:328800/150‑1 (MQ=255)
                                                                                                                                            tGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGGCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTGCATCgg        <  2:221721/150‑1 (MQ=255)
                                                                                                                                            ggCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATACACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGGCCGGCTTCATCgg        <  2:1661489/149‑1 (MQ=255)
                                                                                                                                            cgCGTTTAACCCAGCCGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCCCCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACGGGGATGAATTCTTTTAATGTTTTGCGTCCGGCATCATCgg        >  2:120085/2‑150 (MQ=255)
                                                                                                                                             gCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGt       <  1:1326385/150‑1 (MQ=255)
                                                                                                                                              cGTTTAACCCAGCCGGGTAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg      >  1:1122566/1‑150 (MQ=255)
                                                                                                                                               gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  2:1157773/1‑150 (MQ=255)
                                                                                                                                               gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:1312777/1‑150 (MQ=255)
                                                                                                                                               gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:1571986/1‑150 (MQ=255)
                                                                                                                                                 ttACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAc   >  2:1500172/1‑150 (MQ=255)
                                                                                                                                                  tACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt  >  1:1419802/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  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: