Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,213 C→T 6.2% T121I (ACC→ATC)  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,2130CT 6.2% 480.1 / 15.8 145T121I (ACC→ATC) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base C (68/68);  new base T (5/4);  total (73/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 156 total aligned reads displayed.

TACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTA  >  NC_000913/1327068‑1327359
                                                                                                                                                 |                                                                                                                                                    
tACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCg                                                                                                                                                    <  2:720630/148‑1 (MQ=255)
 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)
       cGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGa                                                                                                                                             <  4:586097/148‑1 (MQ=255)
                   tgtgTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                               <  1:820055/148‑1 (MQ=255)
                     tgtACCCGTGACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                             <  3:113591/148‑1 (MQ=255)
                        aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg‑‑g                                                                                                                                      <  2:682430/136‑1 (MQ=255)
                        aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                               <  3:517480/143‑1 (MQ=255)
                        aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          >  3:120264/1‑148 (MQ=255)
                            gTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          <  2:393483/144‑1 (MQ=255)
                            gTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          >  1:393483/1‑144 (MQ=255)
                            gTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCa                                                                                                                      <  1:477395/148‑1 (MQ=255)
                                     cGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          <  4:244527/135‑1 (MQ=255)
                                     cGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          >  3:244527/1‑135 (MQ=255)
                                      ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                              >  4:413689/1‑130 (MQ=255)
                                       gAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                           >  4:108130/1‑148 (MQ=255)
                                           gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                                                                  >  1:513758/1‑107 (MQ=255)
                                               tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                               >  2:59703/1‑120 (MQ=255)
                                               tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGttt                                                                                                                            <  1:637107/123‑1 (MQ=255)
                                               tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                                                                   <  4:338756/148‑1 (MQ=255)
                                                       tgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                          >  1:429509/1‑133 (MQ=255)
                                                       tgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                          <  2:429509/133‑1 (MQ=255)
                                                          ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                             <  2:204412/111‑1 (MQ=255)
                                                          ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          >  3:779657/1‑114 (MQ=255)
                                                          ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                          <  4:779657/114‑1 (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)
                                                                       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                                                                                                >  3:315634/1‑125 (MQ=255)
                                                                         aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGt                                                                                                <  4:315634/125‑1 (MQ=255)
                                                                                tGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGGTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                         <  3:1080442/125‑1 (MQ=255)
                                                                                 ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                                                                 <  4:554786/148‑1 (MQ=255)
                                                                                  gcgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                                                                            <  3:881497/136‑1 (MQ=255)
                                                                                     cACGCTGGATTGCCGTGGGTCGCTTGGACGTGAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                                                                   <  3:971415/64‑1 (MQ=255)
                                                                                       cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg                                                           >  3:405990/1‑148 (MQ=255)
                                                                                       cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg                                                           <  2:946321/148‑1 (MQ=255)
                                                                                         cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑TCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGa                                                                                              <  3:225575/111‑1 (MQ=255)
                                                                                         cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                  <  4:710000/91‑1 (MQ=255)
                                                                                         cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                  >  3:710000/1‑91 (MQ=255)
                                                                                               tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                   <  4:521708/148‑1 (MQ=255)
                                                                                               tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                   <  3:108130/148‑1 (MQ=255)
                                                                                               tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                   <  4:120264/148‑1 (MQ=255)
                                                                                               tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                   >  4:81630/1‑148 (MQ=255)
                                                                                                  cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                         <  2:207600/139‑1 (MQ=255)
                                                                                                  cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                         >  1:207600/1‑139 (MQ=255)
                                                                                                  cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt                                                >  2:946308/1‑148 (MQ=255)
                                                                                                    tGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCAAAACTGCGTGa                                              <  3:81630/148‑1 (MQ=255)
                                                                                                        tCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                          <  1:215523/148‑1 (MQ=255)
                                                                                                         cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                                                                    >  4:696897/1‑121 (MQ=255)
                                                                                                         cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                                                                    <  3:696897/121‑1 (MQ=255)
                                                                                                         cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTga                                                                >  3:252674/1‑125 (MQ=255)
                                                                                                         cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGa                                         <  1:946308/148‑1 (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‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                                                                        <  3:1034222/107‑1 (MQ=255)
                                                                                                                   ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                     <  1:42123/110‑1 (MQ=255)
                                                                                                                    tAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAATGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                              >  1:497005/1‑148 (MQ=255)
                                                                                                                           tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTa                                                                                                                           >  1:675221/1‑48 (MQ=255)
                                                                                                                           tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGa                                              >  2:232021/1‑125 (MQ=255)
                                                                                                                           tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGa                                              <  1:232021/125‑1 (MQ=255)
                                                                                                                            gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCg                                     <  1:905347/133‑1 (MQ=255)
                                                                                                                            gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCg                                     >  2:905347/1‑133 (MQ=255)
                                                                                                                              ggTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAg                    <  2:505280/148‑1 (MQ=255)
                                                                                                                                 ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                              <  2:664493/39‑1 (MQ=255)
                                                                                                                                 ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                              >  1:664493/1‑39 (MQ=255)
                                                                                                                                 ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGa                                                                                                    >  1:879517/1‑65 (MQ=255)
                                                                                                                                 ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGa                                                                                                    <  2:879517/65‑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)
                                                                                                                                     tgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                  <  2:667682/79‑1 (MQ=255)
                                                                                                                                      gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCa                                                                      >  4:1135176/1‑90 (MQ=255)
                                                                                                                                        tgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                           <  2:224592/51‑1 (MQ=255)
                                                                                                                                          ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                                                                 <  4:1075694/91‑1 (MQ=255)
                                                                                                                                          ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                             >  1:174077/1‑111 (MQ=255)
                                                                                                                                           tCACCATCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                  <  4:232952/73‑1 (MQ=255)
                                                                                                                                           tCACCATCGGTGGTGAACAGCAGCAGACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGAt                  <  1:759337/139‑1 (MQ=255)
                                                                                                                                           tCACCATCGGTGGTGAACAGCAGCAGACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGAt                  >  2:759337/1‑139 (MQ=255)
                                                                                                                                           tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                                                                            >  4:952208/1‑63 (MQ=255)
                                                                                                                                           tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                     >  2:565157/1‑86 (MQ=255)
                                                                                                                                           tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                     <  1:565157/86‑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)
                                                                                                                                              ccatcGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                             <  3:35860/107‑1 (MQ=255)
                                                                                                                                                accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                          >  4:200121/1‑108 (MQ=255)
                                                                                                                                                accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                          <  3:200121/108‑1 (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)
                                                                                                                                                 |                                                                                                                                                    
TACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTA  >  NC_000913/1327068‑1327359

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: