Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 908,438 T→C 6.1% N286S (AAT→AGT)  ltaE ← L‑allo‑threonine aldolase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,4380TC 6.1% 201.2 / 3.3 66N286S (AAT→AGT) ltaEL‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/2);  total (35/31)
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.73e-01

CTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGC  >  NC_000913/908292‑908586
                                                                                                                                                  |                                                                                                                                                    
cTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAttt                                                                                                                                                    >  6:15713/1‑149 (MQ=255)
cTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAttt                                                                                                                                                    >  7:152762/1‑149 (MQ=255)
   aaCGCGCCAGGACTGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGATCGCGCGAGACGTCAAGCGGCGGCACCAGGCGGACAATCGGCGAAGCGTTAATCAGCACGTTTCTCGCTTTCATGTAATCGCCTCCCGCGGCAGCATTttc                                                                                                                                                  <  6:255742/148‑1 (MQ=255)
   aaCGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttct                                                                                                                                                 >  6:131129/1‑149 (MQ=255)
   aaCGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttc                                                                                                                                                  <  5:15713/148‑1 (MQ=255)
   aaCGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttc                                                                                                                                                  <  3:249510/148‑1 (MQ=255)
   aaCGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGGATTAGCCTACCGCGGCAGCATTttct                                                                                                                                                 >  1:316817/1‑149 (MQ=255)
    aCGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttctt                                                                                                                                                <  2:138671/149‑1 (MQ=255)
      gcgcCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTcc                                                                                                                                              <  5:219931/149‑1 (MQ=255)
      gcgcCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTAGCGCGAGACGTCAAGATGCGTCACCAGGAGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTcc                                                                                                                                              >  7:350042/1‑149 (MQ=255)
           agtaaTGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCATCATTTTCTTCCCCGAc                                                                                                                                         <  3:407941/146‑1 (MQ=255)
           aGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGGTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGAc                                                                                                                                         >  1:358914/1‑149 (MQ=255)
             gAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATTCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTGTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGAcgc                                                                                                                                       <  7:306533/149‑1 (MQ=255)
                tGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaa                                                                                                                                    >  7:153136/1‑149 (MQ=255)
                      ccAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaaaca                                                                                                                              <  5:273287/149‑1 (MQ=255)
                            ggcggcGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCata                                                                                                                         <  8:68930/148‑1 (MQ=255)
                              cggcgACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTg                                                                                                                      <  6:377574/149‑1 (MQ=255)
                                     ttCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCt                                                                                                               >  2:320554/1‑149 (MQ=255)
                                       ccGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGCCCTGa                                                                                                             >  3:116480/1‑149 (MQ=255)
                                        cGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGAc                                                                                                            <  5:131129/149‑1 (MQ=255)
                                         cccAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAAACTGCGAGGCGTTAATCAGCACGGTTCTCGCTTTCATGAATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACg                                                                                                           >  5:251869/2‑149 (MQ=255)
                                          ccAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGc                                                                                                          <  5:167542/149‑1 (MQ=255)
                                                 ttCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATcaca                                                                                                    <  2:1379/148‑1 (MQ=255)
                                                 ttCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACAt                                                                                                   <  8:152762/149‑1 (MQ=255)
                                                 ttCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACAt                                                                                                   >  7:42265/1‑149 (MQ=255)
                                                   cgcgcgAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttcttc                                                                                                                                               <  6:230762/103‑1 (MQ=255)
                                                   cgcgcgAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttcttc                                                                                                                                               >  5:230762/1‑103 (MQ=255)
                                                         agaCGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAGTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCt                                                                                           >  7:408355/1‑149 (MQ=255)
                                                              tCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTTGTGTCCTGACGCATCACATACGCGCCTGCtt                                                                                       >  7:152293/1‑148 (MQ=255)
                                                                        tCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGgt                                                                                                                    >  5:85045/1‑109 (MQ=255)
                                                                        tCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGgt                                                                                                                    <  6:85045/109‑1 (MQ=255)
                                                                          aCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTc                                                                                                                 >  4:341456/1‑110 (MQ=255)
                                                                          aCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTc                                                                                                                 <  3:341456/110‑1 (MQ=255)
                                                                            cAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCg                                                                        >  8:288530/1‑149 (MQ=255)
                                                                            cAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCg                                                                        >  2:245281/1‑149 (MQ=255)
                                                                             aGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGc                                                                       <  3:8754/149‑1 (MQ=255)
                                                                               gCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCa                                                                     <  3:174231/149‑1 (MQ=255)
                                                                                cGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCAt                                                                    >  4:331254/1‑149 (MQ=255)
                                                                                   aCAATCGGCGAGGCGTTAATCAGCCCGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATc                                                                   >  3:81192/1‑147 (MQ=255)
                                                                                            gAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGtt                                                        >  7:234591/1‑149 (MQ=255)
                                                                                                 gTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGt                                                   >  4:96651/1‑149 (MQ=255)
                                                                                                  ttAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTg                                                  <  3:331254/149‑1 (MQ=255)
                                                                                                  ttAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTg                                                  <  8:234591/149‑1 (MQ=255)
                                                                                                  ttAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTg                                                  <  1:320554/149‑1 (MQ=255)
                                                                                                   tAATCAGCACGTGTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGCCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTg                                                  <  4:116480/148‑1 (MQ=255)
                                                                                                      tCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCt                                              >  8:222234/1‑149 (MQ=255)
                                                                                                         gCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTcc                                           <  7:288530/149‑1 (MQ=255)
                                                                                                              tttCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTTGTGTCCTGACGCATCACATCAGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTACTGCaa                                      <  8:152293/149‑1 (MQ=255)
                                                                                                                  tcGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGACCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGAGGTCTTCCTGAAAgcgc                                  >  1:188865/1‑149 (MQ=255)
                                                                                                                    gCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAgcgcg                                 >  4:190842/1‑148 (MQ=255)
                                                                                                                    gCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAgcg                                   <  1:245281/146‑1 (MQ=255)
                                                                                                                    gCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATACGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAgcgcg                                 >  7:263694/1‑148 (MQ=255)
                                                                                                                          aTGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGt                           >  8:385203/1‑148 (MQ=255)
                                                                                                                            gTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTAt                        >  1:434184/1‑149 (MQ=255)
                                                                                                                             tATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTAtt                       >  1:244071/1‑149 (MQ=255)
                                                                                                                                  gCCTAACGCGGCTGCATTTTCTTCCGCGAGGCGATCATACAGCATATTGGTGTCCTGACGCATCTCATCCGCGCCTGCTTCTCGCAGCTGCCCCGCCATCCTTGCAGGGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCa                  <  6:251869/149‑1 (MQ=255)
                                                                                                                                   ccTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAg                 >  2:378676/1‑149 (MQ=255)
                                                                                                                                      aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTcc                                           <  1:220282/120‑1 (MQ=255)
                                                                                                                                      aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTcc                                           >  2:220282/1‑120 (MQ=255)
                                                                                                                                       aCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGc              <  2:434184/148‑1 (MQ=255)
                                                                                                                                         gcgGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCata           >  7:180934/1‑149 (MQ=255)
                                                                                                                                         gcgGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCata           >  3:233819/1‑149 (MQ=255)
                                                                                                                                          cgGCAGCATTTTCTTCCCCGACGAGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCatat          <  2:358914/149‑1 (MQ=255)
                                                                                                                                            gcagcaTTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCatatat        >  3:298751/1‑149 (MQ=255)
                                                                                                                                            gcagcaTATTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCatatat        <  5:362697/149‑1 (MQ=255)
                                                                                                                                               gcaTTTTCTTCCCCGACTCGAACAAACAGCATATTGTTGTCCTGACGCATCACATCAGCGCCTGCTTCACGCAGCTGCTCCGACATACAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATccc     >  5:192749/1‑149 (MQ=255)
                                                                                                                                               gcaTTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATccc     >  5:186281/1‑149 (MQ=255)
                                                                                                                                                caTTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCcg    >  3:181218/1‑149 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTAttttt                    >  7:138924/1‑131 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTAttttt                    <  8:138924/131‑1 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCa                  >  7:364363/1‑133 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCa                  <  8:364363/133‑1 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCcggc  <  3:96651/149‑1 (MQ=255)
                                                                                                                                                  ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCcggc  >  3:304836/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGC  >  NC_000913/908292‑908586

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