Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,145,954 C→T 7.6% intergenic (+55/+57) rluC → / ← yceF 23S rRNA pseudouridine(955/2504/2580) synthase/m(7)GTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,145,9540CT 7.6% 237.4 / 7.1 79intergenic (+55/+57)rluC/yceF23S rRNA pseudouridine(955/2504/2580) synthase/m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (33/40);  new base T (3/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

TACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTC  >  NC_000913/1145806‑1146103
                                                                                                                                                    |                                                                                                                                                     
tACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCa                                                                                                                                                      <  1:746631/150‑1 (MQ=255)
tACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCa                                                                                                                                                      >  1:861123/1‑150 (MQ=255)
tACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCa                                                                                                                                                      >  2:1511317/1‑150 (MQ=255)
   tCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCa                                                                                                                                                   <  1:212860/150‑1 (MQ=255)
    cATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTACCGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAg                                                                                                                                                  >  1:118715/1‑150 (MQ=255)
    cATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATGCGGCGTAAACGCCTTATCAg                                                                                                                                                  <  1:616417/150‑1 (MQ=255)
    cATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAg                                                                                                                                                  >  1:959737/1‑150 (MQ=255)
     aTCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAgg                                                                                                                                                 <  1:804042/150‑1 (MQ=255)
      tCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATGCGGCGTAAACGCCTTATCAGgc                                                                                                                                                <  1:1414055/150‑1 (MQ=255)
      tCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATGCGGCGTAAACGCCTTATCAGgc                                                                                                                                                <  1:405542/150‑1 (MQ=255)
      tCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGgc                                                                                                                                                >  2:890195/1‑150 (MQ=255)
      tCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGgc                                                                                                                                                <  1:567212/150‑1 (MQ=255)
           ggACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCa                                                                                                                                           <  2:235843/150‑1 (MQ=255)
              ccGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCa                                                                                                                                        >  1:1339264/1‑150 (MQ=255)
                 gtgaggtgaTGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCg                                                                                                                                     >  2:891348/1‑150 (MQ=255)
                        gaTGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCg                                                                                                                              <  1:1260128/150‑1 (MQ=255)
                          tGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGa                                                                                                                            >  2:1178411/1‑150 (MQ=255)
                             atATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGATGCGTAACGAGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGc                                                                                                                         >  1:337402/2‑150 (MQ=255)
                              tATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGTCTGTTTTCATGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCTATTGCAGGATGCg                                                                                                                        <  2:1312693/150‑1 (MQ=255)
                                tCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGAc                                                                                                                      <  1:739784/150‑1 (MQ=255)
                                     gCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAAc                                                                                                                 <  2:930444/150‑1 (MQ=255)
                                      cTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAAcc                                                                                                                <  1:1112610/150‑1 (MQ=255)
                                      cTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAAcc                                                                                                                >  2:867161/1‑150 (MQ=255)
                                        ccGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACcac                                                                                                              <  1:890195/150‑1 (MQ=255)
                                         cGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACcaca                                                                                                             >  2:609647/1‑150 (MQ=255)
                                         cGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACcaca                                                                                                             >  2:368256/1‑150 (MQ=255)
                                              gatgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCgg                                                                                                        <  2:1030106/150‑1 (MQ=255)
                                                tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                                      <  2:1021307/150‑1 (MQ=255)
                                                tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                                      <  1:481268/150‑1 (MQ=255)
                                                tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                                      <  1:609647/150‑1 (MQ=255)
                                                tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                                      >  2:1250110/1‑150 (MQ=255)
                                                 gAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGAAACCACATCCGgcgc                                                                                                     >  1:973207/1‑150 (MQ=255)
                                                   aGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCaa                                                                                                   >  1:544151/1‑150 (MQ=255)
                                                   aGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCaa                                                                                                   >  2:187380/1‑150 (MQ=255)
                                                    ggTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACg                                                                                                                     <  1:731594/131‑1 (MQ=255)
                                                    ggTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACg                                                                                                                     >  2:731594/1‑131 (MQ=255)
                                                       atGAAGCGTTCTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCaaaaaa                                                                                               >  2:1004061/2‑150 (MQ=255)
                                                        tGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAAc                                                                                              <  1:652229/150‑1 (MQ=255)
                                                         gAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACt                                                                                             >  1:71570/1‑150 (MQ=255)
                                                                  tttGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCa                                                                                    <  1:891348/150‑1 (MQ=255)
                                                                    tGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGc                                                                                  >  1:234229/1‑150 (MQ=255)
                                                                     gCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCg                                                                                 >  1:139915/1‑150 (MQ=255)
                                                                      cAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCTCATCAGGCCGCCAGCACCGAGTTCCGGATGCGACTTAACGACATCCGGCGCAAAAAACTATCCCATCAGCgg                                                                                <  2:1468412/150‑1 (MQ=255)
                                                                         aaGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGtt                                                                             <  2:719051/150‑1 (MQ=255)
                                                                            cTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGttttt                                                                          <  1:56150/150‑1 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACACGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATACCATCAGCGGGttttta                                                                         >  1:917193/1‑149 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                                         <  1:1178411/150‑1 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                                         <  2:292380/150‑1 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                                         <  1:274479/150‑1 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                                         <  1:484960/150‑1 (MQ=255)
                                                                             tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                                         <  1:935276/150‑1 (MQ=255)
                                                                              gCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttcc                                                                        <  2:589402/150‑1 (MQ=255)
                                                                                 tAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttccctt                                                                     >  2:319965/1‑150 (MQ=255)
                                                                                   aCGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttcccttcc                                                                   <  2:544151/150‑1 (MQ=255)
                                                                                    cgcgcgCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttcccttccc                                                                  >  2:759516/1‑150 (MQ=255)
                                                                                    cgcgcgCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTAACCCATCAGCGGGTTTttcccttccc                                                                  >  1:413176/1‑150 (MQ=255)
                                                                                                cAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCt                                                      <  2:64450/150‑1 (MQ=255)
                                                                                                   gcttataAACGTGATAATGGGCGGTTGGCCTGATAAGGCGTTTAGGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACGACATCCGGCGCAAAAAACTATCCCATCATCGGGTTTTTCCCTTCCCGACGTAACATCTGac                                                   <  2:364176/147‑1 (MQ=255)
                                                                                                   gCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGac                                                   <  2:1339264/150‑1 (MQ=255)
                                                                                                     atatAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGacac                                                 >  1:274214/1‑150 (MQ=255)
                                                                                                      tataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGAcaca                                                >  1:647213/1‑150 (MQ=255)
                                                                                                       ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAg                                               >  2:1452029/1‑150 (MQ=255)
                                                                                                              tGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGAtt                                        <  2:71570/150‑1 (MQ=255)
                                                                                                                 tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGc                                     <  2:252146/150‑1 (MQ=255)
                                                                                                                 tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGc                                     <  1:319965/150‑1 (MQ=255)
                                                                                                                     gggCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTa                                 <  2:861123/150‑1 (MQ=255)
                                                                                                                         tggcggCCTGATGCGGCGTAAACGCCTTATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGATGTAACATCTGACACAGCGCGATTAGCGGTAaacc                             >  2:421735/5‑150 (MQ=38)
                                                                                                                             ggCCTGATGCGGCGTAAACGCCTTATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaaccaacc                         >  1:1537078/1‑150 (MQ=255)
                                                                                                                             ggCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaaccaacc                         <  1:187380/150‑1 (MQ=255)
                                                                                                                              gCCTGATGCGGCGTAAACGCCTTATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                        >  1:250186/1‑150 (MQ=255)
                                                                                                                               ccTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                       >  2:17757/1‑150 (MQ=255)
                                                                                                                               ccTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                       <  1:736376/150‑1 (MQ=255)
                                                                                                                               ccTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                       >  1:1267930/1‑150 (MQ=255)
                                                                                                                                cTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGAAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                      >  1:47729/1‑150 (MQ=255)
                                                                                                                                 tGATAAGGCTTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCg                     >  1:477994/1‑150 (MQ=255)
                                                                                                                                 tGATAAGGCGTTTACTCCGCATCAGGCCGCCAGCACCGATTGCCTGATGCGAAGTAACCACATCCGGCTCAAAAAAATATCCAATCAGCGGGTTTTTCCCTTCCCGACTTAACATATGACACAGCGCGATTAGCGGTAAACCAACCAGCg                     <  2:337402/150‑1 (MQ=255)
                                                                                                                                       ggCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAg               <  1:17757/150‑1 (MQ=255)
                                                                                                                                         cGTTTACGCCGCATCAGGCCGCAAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGa             >  2:1035202/1‑150 (MQ=255)
                                                                                                                                           tttACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATc           <  1:1004061/150‑1 (MQ=255)
                                                                                                                                               cgccgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAACCCAACCAGCGTGTTAGGATCACGc       >  2:1009049/1‑150 (MQ=255)
                                                                                                                                               cgccgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGc       >  1:401578/1‑150 (MQ=255)
                                                                                                                                                  cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCTACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGcccc    <  2:47729/150‑1 (MQ=255)
                                                                                                                                                  cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGcccc    <  2:1172436/150‑1 (MQ=255)
                                                                                                                                                  cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGcccc    <  2:234229/150‑1 (MQ=255)
                                                                                                                                                   gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCct   >  1:1111019/1‑150 (MQ=255)
                                                                                                                                                    cATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCctc  <  2:959737/150‑1 (MQ=255)
                                                                                                                                                    cATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGATGTAAACCAACCAGCGTGTTAGGATCACGCCCctc  <  2:118715/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTC  >  NC_000913/1145806‑1146103

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