Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,216 A→G 7.6% D122G (GAT→GGT)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2160AG 7.6% 511.9 / 23.6 157D122G (GAT→GGT) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (74/71);  new base G (5/7);  total (79/78)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

Only 100 of 162 total aligned reads displayed.

ACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAAC  >  NC_000913/1327069‑1327363
                                                                                                                                                   |                                                                                                                                                     
aCAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                                                                       <  4:856493/148‑1 (MQ=255)
  aaaCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGGGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATg                                                                                                                                                     >  2:477395/1‑148 (MQ=255)
    aCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCggtggt                                                                                                                                                   >  3:586097/1‑148 (MQ=255)
    aCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                                                                   <  3:559196/148‑1 (MQ=255)
         aaGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGa                                                                                                                                                 >  4:676546/1‑145 (MQ=255)
               aGTTGTGTACCCGTGACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCaaa                                                                                                                                      >  4:113591/1‑148 (MQ=255)
                  tgtgTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                   <  1:820055/148‑1 (MQ=255)
                    tgtACCCGTGACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                                 <  3:113591/148‑1 (MQ=255)
                       aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                   <  3:517480/143‑1 (MQ=255)
                       aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                   >  4:517480/1‑143 (MQ=255)
                       aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                              >  3:120264/1‑148 (MQ=255)
                          cGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                              >  3:1029350/1‑145 (MQ=255)
                                   ccGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGcc                                                                                                                  >  1:905499/1‑148 (MQ=255)
                                    cGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                              >  3:244527/1‑135 (MQ=255)
                                    cGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                              <  4:244527/135‑1 (MQ=255)
                                     ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                  <  3:413689/130‑1 (MQ=255)
                                     ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                  >  4:413689/1‑130 (MQ=255)
                                      gAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                               >  4:108130/1‑148 (MQ=255)
                                          gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                                                                      <  2:513758/107‑1 (MQ=255)
                                          gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                                                                      >  1:513758/1‑107 (MQ=255)
                                              tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGttt                                                                                                                                <  1:637107/123‑1 (MQ=255)
                                              tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGttt                                                                                                                                >  2:637107/1‑123 (MQ=255)
                                              tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                                                                       <  4:338756/148‑1 (MQ=255)
                                                   cGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTCGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAc                                                                                                                                               <  4:739052/105‑1 (MQ=255)
                                                      tgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                              >  1:429509/1‑133 (MQ=255)
                                                         ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                                 <  2:204412/111‑1 (MQ=255)
                                                         ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                                 >  1:204412/1‑111 (MQ=255)
                                                         ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                              >  3:779657/1‑114 (MQ=255)
                                                         ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                      <  3:85070/122‑1 (MQ=255)
                                                         ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                      >  4:85070/1‑122 (MQ=255)
                                                            aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcg                                                                                          <  3:586749/147‑1 (MQ=255)
                                                            aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcg                                                                                          >  2:1108719/1‑147 (MQ=255)
                                                            aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcg                                                                                          >  4:586749/1‑147 (MQ=255)
                                                             ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                        <  2:905499/148‑1 (MQ=255)
                                                                      aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAttt                                                                               >  3:554786/1‑148 (MQ=255)
                                                                       aaaCTGCGTGGCGCACGCTGGATTACCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                      >  4:66850/1‑140 (MQ=255)
                                                                       aaaCTGCGTGGCGCACGCTGGATTACCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                      <  3:66850/140‑1 (MQ=255)
                                                                        aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGt                                                                                                    <  4:315634/125‑1 (MQ=255)
                                                                        aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGt                                                                                                    >  3:315634/1‑125 (MQ=255)
                                                                               tGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                             >  4:1080442/1‑125 (MQ=255)
                                                                                ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                                                                     <  4:554786/148‑1 (MQ=255)
                                                                                    cACGCTGGATTGCCGTGGGTCGCTTGGACGTGAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                                                                       <  3:971415/64‑1 (MQ=255)
                                                                                      cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg                                                               <  2:946321/148‑1 (MQ=255)
                                                                                        cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                      >  3:710000/1‑91 (MQ=255)
                                                                                        cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                      <  4:710000/91‑1 (MQ=255)
                                                                                              tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                       <  3:108130/148‑1 (MQ=255)
                                                                                              tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                       >  4:81630/1‑148 (MQ=255)
                                                                                              tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                       <  4:120264/148‑1 (MQ=255)
                                                                                                 cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                             <  2:207600/139‑1 (MQ=255)
                                                                                                       tCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                              <  1:215523/148‑1 (MQ=255)
                                                                                                        cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                                                                        <  3:696897/121‑1 (MQ=255)
                                                                                                        cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGa                                             <  1:946308/148‑1 (MQ=255)
                                                                                                                cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                     >  2:745318/1‑84 (MQ=255)
                                                                                                                cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                     <  1:745318/84‑1 (MQ=255)
                                                                                                                cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                       >  1:13188/1‑130 (MQ=255)
                                                                                                                cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                       <  2:13188/130‑1 (MQ=255)
                                                                                                                  ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                                                                            >  4:1034222/1‑107 (MQ=255)
                                                                                                                  ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                         >  2:42123/1‑110 (MQ=255)
                                                                                                                   tAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAATGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                  >  1:497005/1‑148 (MQ=255)
                                                                                                                          tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTa                                                                                                                               <  2:675221/48‑1 (MQ=255)
                                                                                                                          tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTa                                                                                                                               >  1:675221/1‑48 (MQ=255)
                                                                                                                          tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTgg                           >  2:800855/1‑148 (MQ=255)
                                                                                                                                ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                  >  1:664493/1‑39 (MQ=255)
                                                                                                                                ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAg                        <  4:57933/145‑1 (MQ=255)
                                                                                                                                  gctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                        >  1:271723/1‑79 (MQ=255)
                                                                                                                                  gctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAg                                            >  2:813323/1‑123 (MQ=255)
                                                                                                                                  gctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACATGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGt                              <  3:430835/137‑1 (MQ=255)
                                                                                                                                  gctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACATGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGt                              >  4:430835/1‑137 (MQ=255)
                                                                                                                                  gctgctgTTCACCACCGATGGTGAAATG‑‑GCAAACCGTTTAATGCACCCTAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                        <  2:271723/79‑1 (MQ=255)
                                                                                                                                    tgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                      <  2:667682/79‑1 (MQ=255)
                                                                                                                                     gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCg                >  1:418635/1‑148 (MQ=255)
                                                                                                                                       tgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                               <  2:224592/51‑1 (MQ=255)
                                                                                                                                       tgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                               >  1:224592/1‑51 (MQ=255)
                                                                                                                                         ttCACCATCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAg            >  3:299229/1‑148 (MQ=255)
                                                                                                                                         ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                             <  2:89308/67‑1 (MQ=255)
                                                                                                                                         ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                             >  1:89308/1‑67 (MQ=255)
                                                                                                                                         ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                                                                     >  3:1075694/1‑91 (MQ=255)
                                                                                                                                          tCACCATCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                      >  3:232952/1‑73 (MQ=255)
                                                                                                                                          tCACCATCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                      <  4:232952/73‑1 (MQ=255)
                                                                                                                                          tCACCATCGGTGGTGAACAGCAGCAGACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGAt                      <  1:759337/139‑1 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                                                                                >  4:952208/1‑63 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                                                                                <  3:952208/63‑1 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                         <  1:565157/86‑1 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                         >  2:565157/1‑86 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg                                                               <  3:842850/96‑1 (MQ=255)
                                                                                                                                          tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGc           >  1:1030414/1‑148 (MQ=255)
                                                                                                                                           caccatCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCt          <  4:299229/148‑1 (MQ=255)
                                                                                                                                            accaccGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                     >  2:377013/1‑56 (MQ=255)
                                                                                                                                            accaccGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                     <  1:377013/56‑1 (MQ=255)
                                                                                                                                             ccatcGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                 <  3:35860/107‑1 (MQ=255)
                                                                                                                                             ccatcGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                 <  2:437980/107‑1 (MQ=255)
                                                                                                                                             ccatcGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                 >  1:437980/1‑107 (MQ=255)
                                                                                                                                               accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                              <  3:200121/108‑1 (MQ=255)
                                                                                                                                               accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                              >  4:200121/1‑108 (MQ=255)
                                                                                                                                               accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa      >  1:1050210/1‑148 (MQ=255)
                                                                                                                                               accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGAGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                              <  2:903245/108‑1 (MQ=255)
                                                                                                                                               accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGAGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                              >  1:903245/1‑108 (MQ=255)
                                                                                                                                                  gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                      <  2:134817/113‑1 (MQ=255)
                                                                                                                                                  gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                      >  1:134817/1‑113 (MQ=255)
                                                                                                                                                  gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaaa   >  3:1126646/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAAC  >  NC_000913/1327069‑1327363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: