Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,145,939 A→G 8.6% intergenic (+40/+72) 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,9390AG 8.6% 233.1 / 10.6 81intergenic (+40/+72)rluC/yceF23S rRNA pseudouridine(955/2504/2580) synthase/m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (3/4);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCG  >  NC_000913/1145798‑1146084
                                                                                                                                             |                                                                                                                                                 
ttGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGttt                                                                                                                                           <  1:1118199/150‑1 (MQ=255)
 tGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTa                                                                                                                                          <  2:1262978/150‑1 (MQ=255)
 tGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTCGGCCTGATGCGGCGtaaa                                                                                                                                          <  1:9375/150‑4 (MQ=255)
    aGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTAcgc                                                                                                                                       >  1:364176/1‑150 (MQ=255)
     gTTTACTCATCCGGTGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTAcgcc                                                                                                                                      <  2:699552/150‑1 (MQ=255)
        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                                                                                                                             <  1:567212/150‑1 (MQ=255)
              tCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGgc                                                                                                                             >  2:890195/1‑150 (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                                                                                   <  1:609647/150‑1 (MQ=255)
                                                        tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                   >  2:1250110/1‑150 (MQ=255)
                                                        tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                   <  1:481268/150‑1 (MQ=255)
                                                        tgAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGgcg                                                                                   <  2:1021307/150‑1 (MQ=255)
                                                         gAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGAAACCACATCCGgcgc                                                                                  >  1:973207/1‑150 (MQ=255)
                                                           aGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCaa                                                                                >  1:544151/1‑150 (MQ=255)
                                                           aGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCaa                                                                                >  2:187380/1‑150 (MQ=255)
                                                            ggTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACg                                                                                                  >  2:731594/1‑131 (MQ=255)
                                                            ggTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACg                                                                                                  <  1:731594/131‑1 (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                                                      <  1:484960/150‑1 (MQ=255)
                                                                                     tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                      <  1:935276/150‑1 (MQ=255)
                                                                                     tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                      <  1:274479/150‑1 (MQ=255)
                                                                                     tGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTttc                                                      <  2:292380/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                  <  1:319965/150‑1 (MQ=255)
                                                                                                                         tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGc                  <  2:252146/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)
                                                                                                                                             |                                                                                                                                                 
TTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCG  >  NC_000913/1145798‑1146084

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: