Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,278 G→T 7.4% E1012* (GAA→TAA) ‡ rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,2780GT 7.4% 332.8 / 13.3 95E1012* (GAA→TAA) ‡rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base T (4/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

TCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCC  >  NC_000913/4184141‑4184425
                                                                                                                                         |                                                                                                                                                   
tCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgctg                                                                                                                                             >  2:165012/1‑146 (MQ=255)
tCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgctg                                                                                                                                             <  1:165012/146‑1 (MQ=255)
   tCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGa                                                                                                                                        >  2:724959/1‑148 (MQ=255)
         cGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatt                                                                                                                                          >  2:1094262/1‑137 (MQ=255)
         cGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatt                                                                                                                                          <  1:1094262/140‑4 (MQ=255)
                  tCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACt                                                                                                                         >  3:514487/1‑148 (MQ=255)
                     gCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGaa                                                                                                                      >  1:18850/1‑148 (MQ=255)
                      ccGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgat                                                                                                                                           <  4:817510/126‑3 (MQ=255)
                            ccGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAg                                                                                                               <  4:514487/148‑1 (MQ=255)
                              gtgctgtgctGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGtt                                                                                                             >  1:662518/1‑148 (MQ=255)
                                 ctgtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCga                                                                                                          <  1:724959/148‑1 (MQ=255)
                                   gtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgaga                                                                                                        <  2:18850/148‑1 (MQ=255)
                                     gctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgc                                                                                                                                               >  3:871880/1‑107 (MQ=255)
                                     gctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgc                                                                                                                                               <  4:871880/107‑1 (MQ=255)
                                           aGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTa                                                                                                                                  <  2:900721/114‑1 (MQ=255)
                                           aGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTa                                                                                                                                  >  1:900721/1‑114 (MQ=255)
                                                gTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGATGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAacac                                                                                                                  <  1:964368/125‑1 (MQ=255)
                                                gTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAacac                                                                                                                  >  2:964368/1‑125 (MQ=255)
                                                gTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACg                                                                                                                 >  2:818328/1‑126 (MQ=255)
                                                gTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACg                                                                                                                 <  1:818328/126‑1 (MQ=255)
                                                  ggCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGa                                                                                                                               <  2:932280/110‑1 (MQ=255)
                                                  ggCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGa                                                                                                                               >  1:932280/1‑110 (MQ=255)
                                                    cgtgtaaGCTGAGATGCTCGACAAAAGACCGCGCGATGGCTGGTGGTAGATGGGCCTTACAGACTAAGTAATACAAAATAAGCTAGAACAGCTGACTGTGGACTACGTCGAACTGAAACACTAGCTGGAGAAGAAACTCGAAGCCaaa                                                                                       >  4:754747/6‑148 (MQ=255)
                                                    cGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGa                                                                                                                       >  3:29602/1‑116 (MQ=255)
                                                    cGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGa                                                                                                                       <  4:29602/116‑1 (MQ=255)
                                                      ttGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatt                                                                                                                                          >  1:605418/1‑92 (MQ=255)
                                                       tGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGtt                                                                                                             <  1:1099779/123‑1 (MQ=255)
                                                       tGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGtt                                                                                                             >  2:1099779/1‑123 (MQ=255)
                                                       tGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGagaag                                                                                                      >  1:924729/1‑130 (MQ=255)
                                                       tGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGagaag                                                                                                      <  2:924729/130‑1 (MQ=255)
                                                       tGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcgc                                                                                    >  4:229206/1‑148 (MQ=255)
                                                            cTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgat                                                                                                                                           >  1:1003245/1‑86 (MQ=255)
                                                            cTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgat                                                                                                                                           <  2:1003245/88‑3 (MQ=255)
                                                             tGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAACATCAGCTGGAACAGCTGGCTGa                                                                                                                                        >  1:777952/1‑90 (MQ=255)
                                                             tGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGa                                                                                                                                        <  2:777952/90‑1 (MQ=255)
                                                               agaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAgaagaa                                                                                                    <  3:926975/124‑1 (MQ=255)
                                                               agaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAgaagaa                                                                                                    >  4:926975/1‑124 (MQ=255)
                                                                gaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgaga                                                                                                        >  2:426641/1‑119 (MQ=255)
                                                                gaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgaga                                                                                                        <  1:426641/119‑1 (MQ=255)
                                                                    cTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGaa                                                                                             >  4:815703/1‑126 (MQ=255)
                                                                    cTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGaa                                                                                             <  3:815703/126‑1 (MQ=255)
                                                                        aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctg                                                                                                                                             >  4:559551/1‑74 (MQ=255)
                                                                        aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCg                                                                                                           <  2:530671/108‑1 (MQ=255)
                                                                        aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCg                                                                                                           >  1:530671/1‑108 (MQ=255)
                                                                        aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAGGCGAAACGCCGCAAAATCa                                                                          >  4:354322/1‑141 (MQ=255)
                                                                        aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAGGCGAAACGCCGCAAAATCa                                                                          <  3:354322/141‑1 (MQ=255)
                                                                         cAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGaa                                                                                        <  4:1026298/126‑1 (MQ=255)
                                                                         cAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGaa                                                                                        >  3:1026298/1‑126 (MQ=255)
                                                                         cAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGc                                                                  >  3:718161/1‑148 (MQ=255)
                                                                           aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcga                                                                >  4:554682/1‑148 (MQ=255)
                                                                           aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcga                                                                <  3:229206/148‑1 (MQ=255)
                                                                            aCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgag                                                                                                         >  1:522672/1‑106 (MQ=255)
                                                                            aCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgag                                                                                                         <  2:522672/106‑1 (MQ=255)
                                                                                  gcgcgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTg                                                         <  2:662518/148‑1 (MQ=255)
                                                                                     cgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGa                                                                                                                            <  3:1044631/78‑1 (MQ=255)
                                                                                     cgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGa                                                                                                                            >  4:1044631/1‑78 (MQ=255)
                                                                                       aTCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCa                                                      >  4:1144425/1‑146 (MQ=255)
                                                                                       aTCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCa                                                      <  3:1144425/146‑1 (MQ=255)
                                                                                                 ggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCa                                                                          <  1:42555/116‑1 (MQ=255)
                                                                                                 ggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCa                                                                          >  2:42555/1‑116 (MQ=255)
                                                                                                   aGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGTTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTg                                                         >  4:1078390/1‑131 (MQ=255)
                                                                                                   aGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGCCTGATCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTg                                                         <  3:1078390/131‑1 (MQ=255)
                                                                                                          ccTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGtt                                 >  4:262957/1‑148 (MQ=255)
                                                                                                            tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                               <  3:262957/148‑1 (MQ=255)
                                                                                                              aCAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                               >  1:412753/1‑146 (MQ=255)
                                                                                                              aCAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                               <  2:412753/146‑1 (MQ=255)
                                                                                                               cAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTg                                                                                                                        >  2:335056/1‑56 (MQ=255)
                                                                                                               cAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTg                                                                                                                        <  1:335056/56‑1 (MQ=255)
                                                                                                                 gACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGtat                          >  3:17934/1‑148 (MQ=255)
                                                                                                                     aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGa                                         >  1:674548/1‑129 (MQ=255)
                                                                                                                     aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGa                                         <  2:674548/129‑1 (MQ=255)
                                                                                                                        agaAACAAAATCAGCTGGAACAGCTGGCTGAGCCGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTg                                             >  2:99798/1‑122 (MQ=255)
                                                                                                                        agaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTg                                             <  1:99798/122‑1 (MQ=255)
                                                                                                                          aaaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCgg                 <  4:718161/148‑1 (MQ=255)
                                                                                                                           aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGc                                            <  1:132192/120‑1 (MQ=255)
                                                                                                                           aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGc                                            >  2:132192/1‑120 (MQ=255)
                                                                                                                           aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGt                <  4:17934/148‑1 (MQ=255)
                                                                                                                            aCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCAcc                                                    >  3:726913/1‑111 (MQ=255)
                                                                                                                            aCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCAcc                                                    <  4:726913/111‑1 (MQ=255)
                                                                                                                             cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                               >  1:44207/1‑131 (MQ=255)
                                                                                                                             cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                               <  2:44207/131‑1 (MQ=255)
                                                                                                                             cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTa              >  2:612417/1‑148 (MQ=255)
                                                                                                                                aaTCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTg                     <  1:187452/138‑1 (MQ=255)
                                                                                                                                aaTCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTg                     >  2:187452/1‑138 (MQ=255)
                                                                                                                                   cAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCATCCAGGGCGACGATCTGGCACTggg                                                 >  1:1038752/1‑107 (MQ=255)
                                                                                                                                   cAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCggg                                                 <  2:1038752/107‑1 (MQ=255)
                                                                                                                                    aGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCg       <  3:554682/148‑1 (MQ=255)
                                                                                                                                     gCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGt      >  2:722225/1‑148 (MQ=255)
                                                                                                                                      cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt                                                                            <  2:610979/77‑1 (MQ=255)
                                                                                                                                      cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt                                                                            >  1:610979/1‑77 (MQ=255)
                                                                                                                                      cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                          >  4:495459/1‑111 (MQ=255)
                                                                                                                                      cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                          <  3:495459/111‑1 (MQ=255)
                                                                                                                                       tGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAg                              <  3:683004/122‑1 (MQ=255)
                                                                                                                                       tGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAg                              >  4:683004/1‑122 (MQ=255)
                                                                                                                                       tGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTAt    <  1:612417/148‑1 (MQ=255)
                                                                                                                                         gAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATcc  <  1:722225/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
TCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCC  >  NC_000913/4184141‑4184425

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