Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,047,508 G→A 5.9% A146A (GCC→GCT gfcD ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,047,5080GA 5.9% 215.6 / 2.4 67A146A (GCC→GCTgfcDputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (35/28);  new base A (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

AAGCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCT  >  NC_000913/1047362‑1047655
                                                                                                                                                  |                                                                                                                                                   
aaGCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTg                                                                                                                                                     <  5:343194/147‑1 (MQ=255)
aaGCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTg                                                                                                                                                     <  6:80738/147‑1 (MQ=255)
aaGCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGc                                                                                                                                                   <  1:228716/149‑1 (MQ=255)
  gCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGcca                                                                                                                                                 <  1:193427/149‑1 (MQ=255)
  gCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGcca                                                                                                                                                 <  1:198982/149‑1 (MQ=255)
   cTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGccac                                                                                                                                                >  7:13912/1‑149 (MQ=255)
       tatCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGccac                                                                                                                                                >  8:258753/1‑145 (MQ=255)
       tatCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGa                                                                                                                                            <  1:292768/149‑1 (MQ=255)
         tCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGAt                                                                                                                                           <  4:212791/148‑1 (MQ=255)
          cGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCATCAGATAc                                                                                                                                         >  8:261382/1‑149 (MQ=255)
            cccGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTCCCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACt                                                                                                                                        >  6:303622/1‑148 (MQ=255)
              cGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATCGCCCCGCTTGCTTGCTGGCCACCAGATACTCAc                                                                                                                                     <  7:298706/149‑1 (MQ=255)
                      gagaCTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCaac                                                                                                                              >  5:335917/1‑147 (MQ=255)
                          ctctGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGAc                                                                                                                         <  2:355585/149‑1 (MQ=255)
                          ctctGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGAc                                                                                                                         <  8:344760/149‑1 (MQ=255)
                           tctGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGAcc                                                                                                                        >  8:143191/1‑149 (MQ=255)
                                     gACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCCCCGTCAAACAGACCCGTACCAGc                                                                                                               >  1:308484/1‑148 (MQ=255)
                                        aaTATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGccacc                                                                                                                                               <  7:330307/113‑1 (MQ=255)
                                        aaTATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGccacc                                                                                                                                               >  8:330307/1‑113 (MQ=255)
                                        aaTATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGCCCCGTACCAGCAata                                                                                                           >  4:145524/1‑149 (MQ=255)
                                        aaTATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAata                                                                                                           >  5:202538/1‑149 (MQ=255)
                                            tttATCGCTCACCCGGCAAAACGGGTTGCTAATATTGCACCCATAACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCAGCTTGCATGCTGGCCACCAGAAACACACCGACAAACAGACCCGTACCAGCAATATCAc                                                                                                       >  5:51982/1‑149 (MQ=255)
                                             ttATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCAGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACg                                                                                                      >  7:284629/1‑149 (MQ=255)
                                                    tCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTAcc                                                                                                >  3:236461/1‑148 (MQ=255)
                                                      aCCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCaaa                                                                                             >  8:168871/1‑149 (MQ=255)
                                                            cAAAACGGGTTGGTAATATTGCCCGCATTCCCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCAc                                                                                       <  6:308481/149‑1 (MQ=255)
                                                                aCGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAACCAGACACGTACCAGCAATATCACGTTTACCAAACGCCACTTGc                                                                                   <  6:335917/149‑1 (MQ=255)
                                                                aCGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGc                                                                                   >  8:152221/1‑149 (MQ=255)
                                                                   ggTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGt                                                                                <  1:140493/149‑1 (MQ=255)
                                                                    ggtGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCCGCAATATCACGTTTACCAAACGCCACTTGCgg                                                                                 >  3:232650/3‑147 (MQ=255)
                                                                    gTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCATCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTa                                                                               <  7:261382/149‑1 (MQ=255)
                                                                        gTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTAc                                                                                                                   >  8:253134/1‑109 (MQ=255)
                                                                        gTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTAc                                                                                                                   <  7:253134/109‑1 (MQ=255)
                                                                            tatTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa                                                                       <  8:284629/149‑1 (MQ=255)
                                                                                   cGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACcttctt                                                                <  4:236461/149‑1 (MQ=255)
                                                                                    gCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCttcttc                                                               <  3:169245/149‑1 (MQ=255)
                                                                                    gCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCttcttc                                                               >  6:138597/1‑149 (MQ=255)
                                                                                       ttACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca                                                            <  6:202538/149‑1 (MQ=255)
                                                                                           cGGCGTAGCCCCATGTCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCCCCAGATACTCACCGTCAAACAGCCCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCATCTTCCCACAgg                                                        >  5:278318/1‑149 (MQ=255)
                                                                                             gCGTAGCCCCATGCCATCCCGCGGGGGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCCCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCg                                                      <  5:303622/149‑1 (MQ=255)
                                                                                                    cccATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGttt                                                >  1:245047/1‑148 (MQ=255)
                                                                                                    cccATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGttt                                                >  1:250802/1‑148 (MQ=255)
                                                                                                         gCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACg                                                                                           <  6:370359/100‑1 (MQ=255)
                                                                                                         gCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACg                                                                                           >  5:370359/1‑100 (MQ=255)
                                                                                                             tCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAAt                                      <  1:111884/149‑1 (MQ=255)
                                                                                                               ccGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCGCCGTCAAACAGCCCCGTACCAGCAATATCACGTTTACCAAACGCGACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGa                                    >  3:70887/1‑149 (MQ=255)
                                                                                                                cGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAt                                   <  2:308484/149‑1 (MQ=255)
                                                                                                                cGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTAGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAt                                   <  7:40584/149‑1 (MQ=255)
                                                                                                                cGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTAGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAt                                   >  8:40584/1‑149 (MQ=255)
                                                                                                                cGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTAGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAt                                   >  8:34349/1‑149 (MQ=255)
                                                                                                                cGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTAGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAt                                   <  7:34349/149‑1 (MQ=255)
                                                                                                                  aGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGGACCAGCAATATCACGTTTCCCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt                                 >  1:177025/1‑149 (MQ=255)
                                                                                                                      tGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGCTTTGTCt                             <  7:258753/149‑1 (MQ=255)
                                                                                                                      tGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGGCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt                             <  5:47098/149‑1 (MQ=255)
                                                                                                                       gAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTa                                                                               >  6:160702/1‑98 (MQ=255)
                                                                                                                       gAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTa                                                                               <  5:160702/98‑1 (MQ=255)
                                                                                                                       gAAATCAAATGGCCCCCCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTATAAAATCGAATGATTTGTCtt                            >  2:29639/1‑149 (MQ=255)
                                                                                                                         aaTCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                          >  8:395384/1‑149 (MQ=255)
                                                                                                                          aTCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                          >  8:262562/1‑148 (MQ=255)
                                                                                                                             aaaTGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt                       <  5:138597/148‑1 (MQ=255)
                                                                                                                              aaTGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGa                     <  5:160257/149‑1 (MQ=255)
                                                                                                                                  gCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                 >  5:354778/1‑149 (MQ=255)
                                                                                                                                   ccccGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGAt                >  8:241682/1‑149 (MQ=255)
                                                                                                                                    cccGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATc               >  8:217162/1‑149 (MQ=255)
                                                                                                                                             ttgctGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGCCTGATCGCCGCTGa       >  4:307285/1‑148 (MQ=255)
                                                                                                                                               gctGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCAGAAAt    >  2:165347/1‑149 (MQ=255)
                                                                                                                                                ctGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATc   <  7:143191/149‑1 (MQ=255)
                                                                                                                                                cgggCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATc   <  1:289234/147‑1 (MQ=255)
                                                                                                                                                 tGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt  >  1:44276/1‑149 (MQ=255)
                                                                                                                                                 tGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt  >  7:147604/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCTGATATCGCCCGCATCGTGAGACTCTGCGCGATGACAATATTTATCGCTCACCCGGCAAAACGGGTTGGTAATATTGCCCGCATTACCGGCGTAGCCCCATGCCATCCCGAGGGTGAAATCAAATGGCCCCGCTTGCTTGCTGGCCACCAGATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCT  >  NC_000913/1047362‑1047655

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