Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 60,259 T→G 5.6% S30R (AGC→CGC)  rluA ← dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091360,2590TG 5.6% 290.6 / 5.3 100S30R (AGC→CGC) rluAdual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (52/42);  new base G (4/2);  total (56/44)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 102 total aligned reads displayed.

CGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAACCTGC  >  NC_000913/60110‑60408
                                                                                                                                                     |                                                                                                                                                      
cgcgGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCt                                                                                                                                                        >  1:1359959/1‑150 (MQ=255)
 gcgGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATTTCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTc                                                                                                                                                       >  1:1468565/1‑150 (MQ=255)
 gcgGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTc                                                                                                                                                       >  2:318937/1‑150 (MQ=255)
 gcgGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCGCTTCCAGACGACCCGGCACTGACAACAAACCGCTc                                                                                                                                                       <  1:1631991/150‑1 (MQ=255)
     cTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCt                                                                                                                                                   >  2:1758281/1‑150 (MQ=255)
     cTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCt                                                                                                                                                   >  1:611169/1‑150 (MQ=255)
      tttGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCtt                                                                                                                                                  <  2:352248/150‑1 (MQ=255)
        tGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgt                                                                                                                                                <  2:1803399/150‑1 (MQ=255)
        tGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgt                                                                                                                                                >  1:1810050/1‑150 (MQ=255)
          gTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTGGGCttgttg                                                                                                                                              >  1:1648178/1‑150 (MQ=255)
          gTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                              >  2:1649202/1‑150 (MQ=255)
          gTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                              >  2:1278924/1‑150 (MQ=255)
          gTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCGGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                              <  1:969442/150‑1 (MQ=255)
              gcgcCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAGCCGAGCGGTttgttgtcag                                                                                                                                          <  2:812708/150‑5 (MQ=255)
                     tACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAACCCGCGCGGCTTGTTGACCCCCataat                                                                                                                                   >  1:1624630/1‑150 (MQ=255)
                      aCAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCataata                                                                                                                                  <  2:410997/150‑1 (MQ=255)
                       cAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGCCCACCATAatat                                                                                                                                 >  1:276108/1‑150 (MQ=255)
                           cACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTc                                                                                                                             <  2:959836/150‑1 (MQ=255)
                              gccgcTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAACCCGCTCGGCTTGTTGACCACCATAATATGGTCATc                                                                                                                          >  2:912900/1‑150 (MQ=255)
                                 ccTGGTTGCCATTTCAAGACGATGCACCTATTCTGACTGCGGATAATCACGCTGAATTCGTGTCATCACGCTGTCTTTGTGCTCATCCAGACGACCCTGCACTGACAACTAACCGCTCGGCTTGTTGACGACCATAATATGGTCTTCCTg                                                                                                                       <  2:313272/149‑1 (MQ=255)
                                  gtGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACACTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGa                                                                                                                      <  1:80100/149‑1 (MQ=255)
                                     gTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATAc                                                                                                                   <  1:722200/150‑1 (MQ=255)
                                     gTAGCCAAATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAAAGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTCTTGACCACCATAATATGGTCATCCTGATAc                                                                                                                   >  2:550248/1‑150 (MQ=255)
                                       aGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAg                                                                                                                 <  1:21708/150‑1 (MQ=255)
                                          cATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGGCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACTAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGAt                                                                                                              <  1:1429487/150‑1 (MQ=255)
                                             atCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaac                                                                                                           >  1:1212175/1‑150 (MQ=255)
                                              tCCAGACGATGCCCCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACTACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaacc                                                                                                          >  1:1656444/1‑150 (MQ=255)
                                               ccAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaacca                                                                                                         >  1:215801/1‑150 (MQ=255)
                                                    cGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAg                                                                                                    <  1:1725556/150‑1 (MQ=255)
                                                    cGATGCACCGATTCTGCCAGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAg                                                                                                    <  2:329005/150‑1 (MQ=255)
                                                         cACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTc                                                                                               >  1:1587310/1‑150 (MQ=255)
                                                         cACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTc                                                                                               >  1:477581/1‑150 (MQ=255)
                                                         cACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACAGACAACAAACCGCTCGGCTTGTTGAACACCATAATAAGGTCATCCTGATACAGGATAACCAACCAGGGTTc                                                                                               >  1:1373173/1‑150 (MQ=255)
                                                          aCCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTcc                                                                                              <  2:926879/150‑1 (MQ=255)
                                                            cGATTCTGCCTGCGGAAAATCACGCTGAATGCGAGTCATCAAGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTg                                                                                            <  2:1659958/150‑1 (MQ=255)
                                                             gATTCTGCCTGCGGATAATCACGCTGAATGCTCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGAACCGGCACTGACAACAATCCGCTCGGCTTGTTTACCACCAAAATATGGTCAACCTGATACAGTATAAACAACCAGGGTACCTGc                                                                                           >  2:1120802/1‑150 (MQ=255)
                                                              aTTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCg                                                                                          <  2:787847/150‑1 (MQ=255)
                                                               ttCTGCCTGCGGATAATCACGCTGAATGCGCCTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAACCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCgg                                                                                         <  1:1147890/150‑1 (MQ=255)
                                                                 ctgcctgcGGTTTATCACGCTGAATGCGGGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCTGCACTGACAACACACGGGTCGGCTTGTTGACCACCTTAATATGGTCATCCTGATACAGGATAACCATCCAGGGTTCCTGGGGTg                                                                                       <  2:1792581/150‑1 (MQ=255)
                                                                 ctgcctgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGTCCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTg                                                                                       <  1:1683867/150‑1 (MQ=255)
                                                                 ctgcctgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTg                                                                                       <  1:224922/150‑1 (MQ=255)
                                                                 ctgcctgcGGATAATCACGCTGAATGCGCGTCATCACGCTGGCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTg                                                                                       <  1:1163281/150‑1 (MQ=255)
                                                                    cctgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGAt                                                                                    >  1:1391741/1‑150 (MQ=255)
                                                                      tgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTg                                                                                  <  2:1787757/150‑1 (MQ=255)
                                                                         ggATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAg                                                                               >  1:649079/1‑150 (MQ=255)
                                                                          gATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGt                                                                              <  2:474614/150‑1 (MQ=255)
                                                                              aTCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTc                                                                          >  2:644071/1‑150 (MQ=255)
                                                                                cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCa                                                                        <  1:1477838/150‑1 (MQ=255)
                                                                                cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCa                                                                        <  1:15462/150‑1 (MQ=255)
                                                                                cACGCTGAATGCGCGTCATCACGCAGTCTTTGTGGTCTTACAGAAGACCCGGCACTGACAACAAACCGCTCGGAATGTTGACCACAAAAATATTGTCATCCAGCTACATGATAACCAAACAGCGTTCCTGAGGTGGATAGAAGTTTTCCa                                                                        >  2:1349542/1‑150 (MQ=255)
                                                                                   gCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATcc                                                                     >  1:374748/1‑150 (MQ=255)
                                                                                    cTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCTCTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATccc                                                                    >  1:44639/1‑150 (MQ=255)
                                                                                    cTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAGCCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTCCGGAGGATTGTAGTTTTCCATccc                                                                    >  1:1147073/1‑150 (MQ=255)
                                                                                            gcgTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCg                                                           <  2:1349552/150‑1 (MQ=255)
                                                                                            gcgGCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCg                                                           <  2:938601/150‑1 (MQ=255)
                                                                                             cgTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCgg                                                          >  1:790028/1‑150 (MQ=255)
                                                                                              gTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGc                                                         >  1:1403783/1‑150 (MQ=255)
                                                                                               gcatcaCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCt                                                        <  2:1628748/149‑1 (MQ=255)
                                                                                                  tcaCGCTTGATTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCGCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGAGGTGGATTGTAGTTTTCCATCCCCAATTTTCGGCTcc                                                      <  2:1517346/150‑1 (MQ=255)
                                                                                                  acaCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCg                                                     >  1:1098507/2‑150 (MQ=255)
                                                                                                    aCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGtt                                                   >  2:1574351/1‑150 (MQ=255)
                                                                                                     cGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTa                                                  <  1:196205/150‑1 (MQ=255)
                                                                                                      gCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTAc                                                 >  1:1621933/1‑150 (MQ=255)
                                                                                                       cTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACt                                                >  2:488043/1‑150 (MQ=255)
                                                                                                        tGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTg                                               <  2:1815718/150‑1 (MQ=255)
                                                                                                          tCTTTGTGCTCTTCCAGACGACCCGGCACTGACAAGAAACCTCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGAt                                             <  1:1812291/150‑1 (MQ=255)
                                                                                                          tCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGAt                                             >  1:447366/1‑150 (MQ=255)
                                                                                                            tttGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGc                                           <  2:1587310/150‑1 (MQ=255)
                                                                                                                  ctctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACa                                     >  2:276299/1‑150 (MQ=255)
                                                                                                                  ctctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACa                                     >  1:1562632/1‑150 (MQ=255)
                                                                                                                  ctctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCGCCGTTACTGATGCGTTACa                                     >  2:1006126/1‑150 (MQ=255)
                                                                                                                    ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                                   >  1:1761346/1‑150 (MQ=255)
                                                                                                                    ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                                   >  2:631015/1‑150 (MQ=255)
                                                                                                                    ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                                   <  1:318937/150‑1 (MQ=255)
                                                                                                                    ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGAAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                                   >  1:62040/1‑150 (MQ=255)
                                                                                                                      tCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGCTGCGTTACAACGa                                 >  1:1160691/1‑150 (MQ=255)
                                                                                                                      tCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGa                                 >  2:1572572/1‑150 (MQ=255)
                                                                                                                        cAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCACCCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATc                               >  1:1272892/1‑150 (MQ=255)
                                                                                                                        cAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATc                               <  2:1359959/150‑1 (MQ=255)
                                                                                                                        cAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATc                               <  2:43534/150‑1 (MQ=255)
                                                                                                                         aGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCa                              >  1:1244015/1‑150 (MQ=255)
                                                                                                                            cgacCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATAAGGTCATCCTGATACAGGATAACCAACCAGGGTGCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGAGCAAAc                           >  2:1524505/1‑150 (MQ=255)
                                                                                                                                 cGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAgg                      <  2:1098507/150‑1 (MQ=255)
                                                                                                                                       tGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATcc                <  2:269562/150‑1 (MQ=255)
                                                                                                                                       tGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATcc                <  2:477581/150‑1 (MQ=255)
                                                                                                                                       ggtcaacaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCCTcc                >  1:1616417/4‑150 (MQ=255)
                                                                                                                                        gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCa               >  1:1116178/1‑150 (MQ=255)
                                                                                                                                         tcaacaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTACCTGCGGTGGATTGTAGTTTTCGATCCCCATTATCGGCTCCGTTACTGATGCGTTTCAACGATCAAACGCAGGGCATCCAg              >  1:1583782/2‑150 (MQ=255)
                                                                                                                                         tcaacaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGATACTGCGGTGGAATGTAGTTATCCATACCCATTTTCGGCTCCGTTACTGATGCGATACAACGATCAAACGCATCGCATACat              >  2:1673513/2‑149 (MQ=255)
                                                                                                                                         acaacaAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAg              <  2:1648178/150‑1 (MQ=255)
                                                                                                                                          caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGa             >  2:1782073/1‑150 (MQ=255)
                                                                                                                                          caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGa             >  2:656808/1‑150 (MQ=255)
                                                                                                                                             aaaaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGc          <  1:921663/149‑1 (MQ=255)
                                                                                                                                                aCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACCACGATCAAACGCAGGGCATCCAGACGCCaa       >  2:112981/1‑150 (MQ=255)
                                                                                                                                                aCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCaa       >  2:1703071/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGAAGCGTTACAACGATCAAACGCAGGGCATACAGACGCCAAc      >  2:78787/1‑150 (MQ=255)
                                                                                                                                                  cGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGGTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcc     <  2:1548331/150‑1 (MQ=255)
                                                                                                                                                   gCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcct    >  2:1779982/1‑150 (MQ=255)
                                                                                                                                                    cTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTAACAACGATCAAACGCAGGGCATCCAGACGCCAAcct  >  1:322638/1‑149 (MQ=255)
                                                                                                                                                     tCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcctgc  >  1:489730/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAACCTGC  >  NC_000913/60110‑60408

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