Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,948,276 G→A 5.2% C64Y (TGC→TAC)  yifE → UPF0438 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,948,2760GA 5.2% 255.9 / 2.4 76C64Y (TGC→TAC) yifEUPF0438 family protein
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base A (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AAACATTATCCACGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAA  >  NC_000913/3948128‑3948424
                                                                                                                                                    |                                                                                                                                                    
aaaCATTATCCACGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATg                                                                                                                                                      <  3:554025/149‑1 (MQ=255)
      tATCCACGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTgg                                                                                                                                                <  4:18403/149‑1 (MQ=255)
          cACGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGaa                                                                                                                                            >  4:564661/1‑149 (MQ=255)
           aCGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAAc                                                                                                                                           >  3:99304/1‑149 (MQ=255)
             gTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGt                                                                                                                                         >  1:93287/1‑149 (MQ=255)
               ggATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGa                                                                                                                                       <  3:281390/149‑1 (MQ=255)
               ggATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGa                                                                                                                                       <  3:460048/149‑1 (MQ=255)
                gATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAg                                                                                                                                      <  3:319669/149‑1 (MQ=255)
                  ttCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGcc                                                                                                                                    >  1:274864/1‑149 (MQ=255)
                  ttCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGcc                                                                                                                                    >  1:438819/1‑149 (MQ=255)
                     tctcGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGt                                                                                                                                 <  4:162435/149‑1 (MQ=255)
                     tctcGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGt                                                                                                                                 <  2:238354/149‑1 (MQ=255)
                     tctcGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGt                                                                                                                                 >  3:501534/1‑149 (MQ=255)
                     tctcGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGt                                                                                                                                 >  3:197109/1‑149 (MQ=255)
                       tcGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGa                                                                                                                               <  4:434991/149‑1 (MQ=255)
                            aTGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGAGCTTGGCAAACGCGAGCCGTTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGTCAGaa                                                                                                                          <  4:484223/149‑1 (MQ=255)
                             tGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGt                                                                                                                                         >  3:393827/1‑133 (MQ=255)
                             tGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGt                                                                                                                                         <  4:393827/133‑1 (MQ=255)
                               gTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                                                                                                       >  1:502083/1‑149 (MQ=255)
                                       aCCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgt                                                                                                               <  1:141151/149‑1 (MQ=255)
                                          aTCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAACCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                                                                                            <  1:425434/149‑1 (MQ=255)
                                           tCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTgg                                                                                                           >  2:495263/1‑149 (MQ=255)
                                                   gCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATACCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                                                                                   >  1:252924/1‑149 (MQ=255)
                                                           gCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgc                                                                                           <  2:64191/149‑1 (MQ=255)
                                                             ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGt                                                                                         <  1:313532/149‑1 (MQ=255)
                                                             ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGt                                                                                         >  3:519351/1‑149 (MQ=255)
                                                                    tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc                                                                                  >  2:581600/1‑149 (MQ=255)
                                                                    tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc                                                                                  <  1:548333/149‑1 (MQ=255)
                                                                    tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTCGTCCAAGTATATGACGCGTATTAAGc                                                                                  <  3:475766/149‑1 (MQ=255)
                                                                          ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                                                                                            <  1:521892/117‑1 (MQ=255)
                                                                          ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                                                                                            >  2:521892/1‑117 (MQ=255)
                                                                          ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaa                                                                            >  2:443430/1‑149 (MQ=255)
                                                                          ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATACCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaa                                                                            >  2:223622/1‑149 (MQ=255)
                                                                              gCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACg                                                                        >  1:115026/1‑149 (MQ=255)
                                                                              gCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACg                                                                        >  2:360768/1‑149 (MQ=255)
                                                                                cTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTc                                                                               >  1:13345/1‑140 (MQ=255)
                                                                                cTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTc                                                                               <  2:13345/140‑1 (MQ=255)
                                                                                cTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtt                                                                      >  4:547316/1‑149 (MQ=255)
                                                                                cTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtt                                                                      >  1:268279/1‑149 (MQ=255)
                                                                                 ttCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt                                                                     >  2:142832/1‑149 (MQ=255)
                                                                                 ttCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt                                                                     >  4:372323/1‑149 (MQ=255)
                                                                                 ttCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt                                                                     >  4:207455/1‑149 (MQ=255)
                                                                                  tCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt                                                                    <  3:551340/149‑1 (MQ=255)
                                                                                    aaTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTca                                                                  >  1:157411/1‑149 (MQ=255)
                                                                                          ttGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCt                                                            >  2:406853/1‑149 (MQ=255)
                                                                                                cTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCcgg                                                      <  3:232029/149‑1 (MQ=255)
                                                                                                 ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTGTCACACCCTTTCcggc                                                     >  2:171253/1‑149 (MQ=255)
                                                                                                   ggCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCcggcgg                                                   <  2:274864/149‑1 (MQ=255)
                                                                                                          gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCg                                            >  1:58325/1‑149 (MQ=255)
                                                                                                            gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCa                                          <  2:383059/149‑1 (MQ=255)
                                                                                                             aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGGTTTCACACCCTTTCCGGCGGTAAACCGCAg                                         <  1:581600/149‑1 (MQ=255)
                                                                                                              gCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAgg                                        <  4:197109/149‑1 (MQ=255)
                                                                                                                 ggTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTg                                     >  3:45137/1‑149 (MQ=255)
                                                                                                                   ttACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGaa                                   <  1:349860/149‑1 (MQ=255)
                                                                                                                    tACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAg                                  <  3:207455/149‑1 (MQ=255)
                                                                                                                    tACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAg                                  <  2:58325/149‑1 (MQ=255)
                                                                                                                    tACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAg                                  <  2:22544/149‑1 (MQ=255)
                                                                                                                    tACCGAGGAAGAGAAACTCTTCGTAGCAGTATACCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAg                                  <  1:223622/149‑1 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTg                               <  2:577637/147‑1 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTg                               >  1:577637/1‑147 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCt                             >  1:34492/1‑149 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCt                             <  3:241205/149‑1 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCt                             >  1:555381/1‑149 (MQ=255)
                                                                                                                         aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGCTAAACCGCAGGTTGAAGGTGCt                             >  4:466398/1‑149 (MQ=255)
                                                                                                                           gAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGa                           <  4:87504/149‑1 (MQ=255)
                                                                                                                            aaGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGaa                          <  4:507999/149‑1 (MQ=255)
                                                                                                                              gagaAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGa                        <  3:527714/149‑1 (MQ=255)
                                                                                                                               agaAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGAc                       <  1:561432/149‑1 (MQ=255)
                                                                                                                                   aCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTaca                   >  4:500187/1‑149 (MQ=255)
                                                                                                                                    ctctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTacac                  <  3:372323/149‑1 (MQ=255)
                                                                                                                                    ctctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTacac                  >  1:218334/1‑149 (MQ=255)
                                                                                                                                     tctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACAcc                 >  1:192021/1‑149 (MQ=255)
                                                                                                                                      ctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCg                <  2:93287/149‑1 (MQ=255)
                                                                                                                                        gcttAGCAGTATGCCGTTGCGATCTTGTGCCAGTGACAGAAGCAGCACGCGCGTGGTTCACGTATATGACGCGTATTAAGCCTCCAAAACTTTTTCACGCCCTTCCCTGCGGTAAACCGAAGGTTGAAGGTTCTGTAGAGTAAATCGAt              <  4:66126/146‑1 (MQ=255)
                                                                                                                                         cGTAGCAGTATACCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGAtt             <  2:252924/149‑1 (MQ=255)
                                                                                                                                            aGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCACAGGTTGAAGGTGCTGAAGACTACACCGATTCTg          >  3:244420/1‑149 (MQ=255)
                                                                                                                                             gCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGa         >  3:155468/1‑149 (MQ=255)
                                                                                                                                                    gCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATtaa  <  1:492630/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAACATTATCCACGTGGATTCTCTCGTCATGGTGATTTCACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAA  >  NC_000913/3948128‑3948424

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: ‑