Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,650 G→A 8.3% P101L (CCG→CTG)  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,6500GA 8.3% 361.5 / 15.0 123P101L (CCG→CTG) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base G (58/53);  new base A (4/6);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 124 total aligned reads displayed.

ATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGC  >  NC_000913/1194501‑1194796
                                                                                                                                                      |                                                                                                                                                  
aTAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCg                                                                                                                                                    <  2:918763/150‑1 (MQ=255)
  aaCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGc                                                                                                                                                  >  2:606717/1‑150 (MQ=255)
      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)
        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                                                                                                                                        <  1:742574/144‑1 (MQ=255)
                  cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                  <  1:1435613/150‑2 (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)
                      cAGCCCCGCGCCGGCGGGCAGGGCAGGGCCATCATTTAAGGATACGCCATTGCGCAAGGCCTCAAGAGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAccggaaa                                                                                                                              <  1:295988/150‑6 (MQ=255)
                        gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                              >  1:876744/1‑132 (MQ=255)
                        gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTggg                                                                                                                                              <  2:876744/132‑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)
                                    ggggCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGTCTTCTAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                <  2:765997/149‑1 (MQ=255)
                                    cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                <  2:1362462/150‑1 (MQ=255)
                                    cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                <  2:1593065/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)
                                                ggCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                    >  2:617140/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)
                                                              gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                      >  1:1721606/1‑150 (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)
                                                                            cAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCAAATAATAGATTTTTCCGGTGCGTTTACCCGCCTGGGTTAAACGCGCATGCAACGCGCAGATATTGGTCAGCACCAGCAACCCTTCTCTATCTCGGTCAAGGCGAcc                                                                        >  1:389739/1‑150 (MQ=255)
                                                                               ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                     >  1:1119019/1‑150 (MQ=255)
                                                                               ggCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                     <  1:954672/150‑1 (MQ=255)
                                                                                     aaGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATa                                                               <  2:737515/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)
                                                                                                     ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGGGCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                               <  2:206437/150‑1 (MQ=255)
                                                                                                        aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATg                                            <  2:1360/150‑1 (MQ=255)
                                                                                                          aCCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCTCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGaa                                          >  2:84717/1‑150 (MQ=255)
                                                                                                           ccTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGCTCCGGGATGAAt                                         >  2:309012/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:1410848/1‑150 (MQ=255)
                                                                                                                   cTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                                 <  2:197903/150‑1 (MQ=255)
                                                                                                                      cacaTAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                              >  1:332884/1‑150 (MQ=255)
                                                                                                                         ataataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                           <  2:798279/150‑1 (MQ=255)
                                                                                                                           aataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGc                         >  1:1416647/1‑150 (MQ=255)
                                                                                                                             taGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                       >  2:1427192/1‑150 (MQ=255)
                                                                                                                              aGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                      >  1:157329/1‑150 (MQ=255)
                                                                                                                                aTTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                    <  1:874758/150‑1 (MQ=255)
                                                                                                                                  ttttCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGc                  <  2:836592/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)
                                                                                                                                        ggTGCGATTACCCGGCTCGGTTAAACGAGCCTGCAACGCGCCGTTATTGGTCAGCACCACCAACCCTACGCAAACGCGGTCAAGGCGACCAGCTGCATAAACACCCAGAACCCGGATGAATTCTTTTAATGTATTGCGTCCGGCTTCATc            >  2:48865/1‑150 (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       >  1:1312777/1‑150 (MQ=255)
                                                                                                                                             gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa       >  1:1571986/1‑150 (MQ=255)
                                                                                                                                             gTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa       >  2:1157773/1‑150 (MQ=255)
                                                                                                                                               ttACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAc     >  2:1500172/1‑150 (MQ=255)
                                                                                                                                                tACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt    >  1:1419802/1‑150 (MQ=255)
                                                                                                                                                 aaCCCAGCCGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg   >  2:352680/2‑150 (MQ=255)
                                                                                                                                                 aaCCCAGCCGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg   >  2:1120655/2‑150 (MQ=255)
                                                                                                                                                  aCCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg   <  2:40064/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGC  >  NC_000913/1194501‑1194796

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: