Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,245 T→C 9.3% L375L (TTA→CTA) ‡ icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,2450TC 9.3% 249.8 / 12.4 86L375L (TTA→CTA) ‡icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base C (4/4);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATT  >  NC_000913/1196098‑1196392
                                                                                                                                                   |                                                                                                                                                   
gCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                                                                                                                     >  5:102006/1‑148 (MQ=255)
 cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACtt                                                                                                                                                    >  3:180760/1‑148 (MQ=255)
    aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                                                <  2:400730/149‑1 (MQ=255)
    aCATCGGTGACGAATGCCCACTGTTAGAAGCCAAACACGGAACTTCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTCGCTCTATTAATCACTCAGCAGAGATGATGCTGCGCCACATGGGTTGGACCGAACCGGCTGACTTAAtt                                                                                                                                                >  4:391590/1‑149 (MQ=255)
     cATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGGCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                                                                                                               <  4:117513/149‑1 (MQ=255)
          gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                                          <  5:193366/149‑1 (MQ=255)
           tGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCCGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                                                                                         <  5:340659/149‑1 (MQ=38)
           tGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                                         <  4:151684/149‑1 (MQ=255)
              cGAATGCGCCCTGTTTGAAGCCACCCTCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCCGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                                                                                                      <  5:255641/149‑1 (MQ=255)
               gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAACTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                                                                                                                     <  3:217300/149‑1 (MQ=255)
               gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                                                                                                                     <  6:97154/149‑1 (MQ=255)
               gAATGCGCCCTGTTGGAAGCCACCCACTGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                                                                                                                     <  3:222785/149‑1 (MQ=255)
                aaTGCGCACTGTTTGAAGCCACCAACGGTACTGCGACGAAATATGCCGGACAGGACAAACTAAATCCTGGATCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATCGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                                                                                                    >  8:365224/1‑149 (MQ=255)
                 aTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                                                                                                   <  3:168815/149‑1 (MQ=255)
                  tGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGa                                                                                                                                  <  6:102006/149‑1 (MQ=255)
                    cgcCCGGTTTGAATCCACACACGTTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGATATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                                                                                                                <  2:112677/149‑1 (MQ=38)
                     gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                                                                                                               >  3:3355/1‑149 (MQ=255)
                            ttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                                                                                                               <  1:37571/126‑1 (MQ=255)
                            ttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                                                                                                               >  2:37571/1‑126 (MQ=255)
                              gAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                                                                                                                      >  3:259464/1‑149 (MQ=255)
                               aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                                                                     >  8:382786/1‑149 (MQ=255)
                                    acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgcca                                                                                                                 >  8:166549/1‑146 (MQ=25)
                                    acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                                         >  1:2580/1‑124 (MQ=255)
                                    acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                                         <  2:2580/124‑1 (MQ=255)
                                      ccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCACCGCGaaaa                                                                                                              >  4:397197/1‑149 (MQ=255)
                                      ccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                                                                                              >  3:252807/1‑149 (MQ=255)
                                               tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAcc                                                                                                      >  6:306938/1‑148 (MQ=255)
                                                 cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                                   <  7:382786/149‑1 (MQ=255)
                                                   ccGAAATATGCCGGGCAGGACCTAGTAAATCGTGGATCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAATCGGCTGACTTAATTGTTAAAGGTATGGAAGTCGCAATCAACGCGAAAACCGTAACCTATGa                                                                                                 <  1:270307/149‑1 (MQ=255)
                                                    cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                                         >  2:55357/1‑108 (MQ=255)
                                                    cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                                         <  1:55357/108‑1 (MQ=255)
                                                    cGAAATATGCCGGACAGGACAAAGTAAATCCTGGCTCTACTATTCACTCCGCAGAGATGATGCTGAGCCACAAGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGTTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATTAc                                                                                                >  2:297264/1‑149 (MQ=255)
                                                        atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCg                                                                                            <  4:252807/149‑1 (MQ=255)
                                                          atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                                                          >  8:291202/1‑149 (MQ=255)
                                                           tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGc                                                                                         >  5:64825/1‑149 (MQ=255)
                                                              cGGTCAGTACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGCTTGGACCGAAGCGGCTGACTTAATTGTTAAAGTTATGGAAGTCGCAATCAACGCGAAAACCGTAACCTATGCCATCGAGCGTc                                                                                      >  2:155627/1‑149 (MQ=255)
                                                                gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTg                                                                                    <  4:131111/149‑1 (MQ=255)
                                                                   aGGTCAAAGTAAATCCTGGCTCTATTATTCTCTCCGATGATATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                                                  <  3:10870/148‑1 (MQ=255)
                                                                    ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatg                                                                                 >  2:201947/1‑148 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACCTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                                                                                                        <  5:230736/92‑1 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                                                                                                        >  6:230736/1‑92 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                                                      <  4:284957/142‑1 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                                                      >  3:284957/1‑142 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                                                      >  1:70597/1‑142 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                                                      <  2:70597/142‑1 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                               >  2:55132/1‑149 (MQ=255)
                                                                     gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                               >  2:86822/1‑149 (MQ=255)
                                                                       cAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAACTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                                                                                                             <  5:302027/117‑1 (MQ=255)
                                                                       cAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAACTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                                                                                                             >  6:302027/1‑117 (MQ=255)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                                                                                                      >  1:126513/1‑88 (MQ=18)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                                                                                                      <  2:126513/88‑1 (MQ=18)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                                                                              <  2:314596/128‑1 (MQ=14)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                                                                              <  2:308902/128‑1 (MQ=14)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                                                                              >  1:314596/1‑128 (MQ=14)
                                                                           gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                                                                              >  1:308902/1‑128 (MQ=14)
                                                                                  cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                                                  <  5:405327/133‑1 (MQ=255)
                                                                                  cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                                                  >  6:405327/1‑133 (MQ=255)
                                                                                  cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActg                                                                  <  7:365224/149‑1 (MQ=255)
                                                                                   tGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                                                                 >  1:419305/1‑149 (MQ=255)
                                                                                     gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgctg                                                               <  4:180760/149‑1 (MQ=255)
                                                                                       tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGa                                                              >  7:281181/1‑148 (MQ=255)
                                                                                       tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGa                                                              >  8:165499/1‑148 (MQ=255)
                                                                                           ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGt                                                         <  7:291202/149‑1 (MQ=255)
                                                                                           ttattCTCACCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGACGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGAGAAAACCGTAACCTATGACTTCGAGCGTCAGATGGATGGCGCTAACCTGCTGAAATg                                                          >  6:101126/1‑148 (MQ=21)
                                                                                              ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                                                        >  8:409544/1‑147 (MQ=255)
                                                                                                ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                                                      >  2:271225/1‑147 (MQ=255)
                                                                                                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCagag                                                   >  1:160471/1‑149 (MQ=255)
                                                                                                   tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCTAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                                                    <  7:175537/146‑1 (MQ=255)
                                                                                                   tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                                                    <  7:409544/146‑1 (MQ=255)
                                                                                                   tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt                                                 >  5:43754/1‑149 (MQ=255)
                                                                                                     cGCTGAGATGATGCTGCGCCACATTGGTTTGACCGAAGCGTCTGACTTAATTGTTAAAGGTATTGAAGGCGCAATCAACGAGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                               <  6:64825/149‑1 (MQ=18)
                                                                                                          agaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAGTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAc                                          >  7:375093/1‑149 (MQ=255)
                                                                                                          agaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTg                                            <  4:3355/147‑1 (MQ=255)
                                                                                                          agaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAc                                          >  8:209202/1‑149 (MQ=255)
                                                                                                             tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc                                        >  1:28048/1‑148 (MQ=255)
                                                                                                              gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                                                    <  7:116010/135‑1 (MQ=255)
                                                                                                              gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                                                    >  8:116010/1‑135 (MQ=255)
                                                                                                                     cgcCACATGGGTTGGACCGAATCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                 <  8:281181/147‑1 (MQ=255)
                                                                                                                     cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                 <  4:259464/147‑1 (MQ=255)
                                                                                                                     cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                 <  6:43754/147‑1 (MQ=255)
                                                                                                                       ccACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                             <  1:201947/149‑1 (MQ=255)
                                                                                                                        cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaaa                            >  6:271285/1‑149 (MQ=255)
                                                                                                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt                    <  2:419305/149‑1 (MQ=255)
                                                                                                                                     ccGAATCTGCTGACTTAATTGTTAATGGTATGGAAGGCGCAATCTACTCGAAAACCGTAACCTATGAATTCGAGCGGCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGt               <  1:297264/149‑1 (MQ=18)
                                                                                                                                      cGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTa              <  1:86822/149‑1 (MQ=255)
                                                                                                                                             gcagACTTAATTCTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAAGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGCGTTTGGTGACGCGATCATCGAAAACATGAAAAGCCGTAGTTTGt        >  4:323723/4‑148 (MQ=18)
                                                                                                                                             gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt       <  5:271285/149‑1 (MQ=255)
                                                                                                                                              cttaCTTAATTGTTAAAGGTATGGAAGGCGCAATAAACGCGTAATACTTAACCTATGACTTCGTGCTTCTGATGGATGGCGCTATACTGCTGAAATGTTCAGAGTTTTGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt       <  3:323723/145‑1 (MQ=18)
                                                                                                                                                   ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt  >  4:275704/1‑148 (MQ=255)
                                                                                                                                                   ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt  >  7:192540/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATT  >  NC_000913/1196098‑1196392

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