Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,362,934 C→A 14.2% P43H (CCT→CAT)  yhcD → putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,9340CA14.2% 264.1 / 11.8 99P43H (CCT→CAT) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base C (32/51);  new base A (13/3);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCA  >  NC_000913/3362790‑3363081
                                                                                                                                                |                                                                                                                                                   
tGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTa                                                                                                                                                   <  3:18959/147‑1 (MQ=255)
tGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTa                                                                                                                                                   <  6:81968/147‑1 (MQ=255)
  ccTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAat                                                                                                                                                 <  7:191158/147‑1 (MQ=255)
    tGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAatata                                                                                                                                              >  3:1473/1‑148 (MQ=255)
    tGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAatata                                                                                                                                              >  1:33106/1‑148 (MQ=255)
      gAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAatata                                                                                                                                              <  8:321853/146‑1 (MQ=255)
           aCAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTaa                                                                                                                                                  >  3:330014/1‑137 (MQ=255)
           aCAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTaa                                                                                                                                                  <  4:330014/137‑1 (MQ=255)
               ggATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTAttt                                                                                                                                   >  7:200526/1‑148 (MQ=255)
                    ttAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTccc                                                                                                                              >  7:234765/1‑148 (MQ=255)
                     tAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCAGCTCTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATACTAATATATCGCATTATTTTTCCCg                                                                                                                             >  7:318567/1‑148 (MQ=255)
                       aaaaaaaCGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGtt                                                                                                                           <  5:40346/148‑1 (MQ=255)
                          aaaaCGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAg                                                                                                                        <  8:200526/148‑1 (MQ=255)
                           aaaCGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGc                                                                                                                       >  5:97177/1‑148 (MQ=255)
                                   cTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGttttt                                                                                                               <  6:97177/148‑1 (MQ=255)
                                   cTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGtttt                                                                                                                <  1:194239/147‑1 (MQ=255)
                                          acacCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAat                                                                                                                                                 <  3:192140/107‑1 (MQ=255)
                                          acacCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAat                                                                                                                                                 >  4:192140/1‑107 (MQ=255)
                                          acacCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATATCGCATTATTTTTCAAGTTCAGCCCGTTTTTTGCCAgg                                                                                                        >  4:235799/1‑148 (MQ=255)
                                               aTTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGaa                                                                                                    >  3:295174/1‑147 (MQ=255)
                                               aTTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAat                                                                                                   <  6:192247/148‑1 (MQ=255)
                                                ttGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATTACCCATTAAATTTAAAAAAAAACCCGTTTTTTGCCAGGTGAata                                                                                                  >  2:214632/1‑148 (MQ=255)
                                                 tGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGTGAAACATTAAAATCACTAGGTGAAGATCAAAATATAACGCATAAAATTAAACAATCAACCCGTTTTTTTACAGGTGAaaat                                                                                                 >  2:107393/1‑148 (MQ=255)
                                                   gTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTc                                                                                               <  4:1473/148‑1 (MQ=255)
                                                     tttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTc                                                                                               <  4:60644/146‑1 (MQ=255)
                                                     tttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCAc                                                                                             >  3:165682/1‑148 (MQ=255)
                                                     tttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCACGTTTTTTGCCAGGTGAATATTc                                                                                               >  3:60644/1‑146 (MQ=255)
                                                     tttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCAAATATATCGCATTAATTTTCACGTTAAGCCCGTTTTTTGCCAGGTGAATATTCAc                                                                                             >  6:94001/1‑148 (MQ=255)
                                                     tttGCCTTTTCTCCCCCAGACAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCAGTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCAc                                                                                             >  5:316403/1‑148 (MQ=255)
                                                      ttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATATAATTAAAAATAAAACCCGTTTTTTGCCAGGTGATTATTCACt                                                                                            >  8:188832/1‑148 (MQ=255)
                                                       ggCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTCCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTg                                                                                           <  6:81143/146‑1 (MQ=255)
                                                         ccTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGa                                                                                          >  3:214841/1‑147 (MQ=255)
                                                         ccTTTTATCCCCCAGCCAATGCAGATTGTATAGAGATAGCCGCTGTTGATTTTGATCGGGAAACAATAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGAt                                                                                         <  4:225642/148‑1 (MQ=255)
                                                           ttttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAATAAACACGTTTTTTGCCAGGTGAATATTCACTGATAg                                                                                       >  2:74352/1‑148 (MQ=255)
                                                            tttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTAttttt                                                                                                                                 <  1:102148/105‑1 (MQ=255)
                                                            tttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTAttttt                                                                                                                                 >  2:102148/1‑105 (MQ=255)
                                                            tttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCACGTTTTTTGCCAGGTGAATATTCAATGATAg                                                                                       >  5:85689/1‑147 (MQ=255)
                                                            tttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATATCGCATTATTTTTAACTTTAAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAg                                                                                       >  2:289530/1‑147 (MQ=255)
                                                             ttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAg                                                                                       <  8:237256/146‑1 (MQ=255)
                                                             ttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGt                                                                                      <  3:235799/147‑1 (MQ=255)
                                                              tctcCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                     <  8:202609/147‑1 (MQ=255)
                                                              tctcCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATTAATTTTAAATATCAGCACGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                     >  8:31948/1‑147 (MQ=255)
                                                                  ccccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATc                                                                                   <  5:363110/145‑1 (MQ=255)
                                                                  cccaAGCCAATGCAGATGGTATAGAGATTGGAGCGTTTGATTTTGATCGGGAAACTTTAATATGACTAGGTGTAGATCCTAATATATCGCATTTTTTTTCCCGTTCAGGCCTTTTTTTGCCAGGTGAATATTCACTGATAGTATc                                                                                   <  6:85689/145‑1 (MQ=255)
                                                                   cccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGtt                                                                                                                  >  7:25532/1‑113 (MQ=255)
                                                                   cccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGtt                                                                                                                  <  8:25532/113‑1 (MQ=255)
                                                                    ccAGCCAATGCAGATGGAATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAACTCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaa                                                                              <  7:210818/148‑1 (MQ=255)
                                                                     cAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTTATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATAACGCATTATATTTAAAAAAAAAAACGTTATTTGCCAGGTGAATATTAACTGATAGTATCAGTaaa                                                                             >  2:362874/1‑148 (MQ=255)
                                                                     cAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTTAAAATTCAACCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaaa                                                                             >  4:252030/1‑148 (MQ=255)
                                                                       gCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATg                                                                           <  1:214632/148‑1 (MQ=255)
                                                                          aaTGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTg                                                                                                             >  1:237317/1‑111 (MQ=255)
                                                                          aaTGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTg                                                                                                             <  2:237317/111‑1 (MQ=255)
                                                                            tGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGa                                                                       >  5:261435/1‑147 (MQ=255)
                                                                             gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAgg                                                                                                        >  1:222387/1‑113 (MQ=255)
                                                                             gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAgg                                                                                                        <  2:222387/113‑1 (MQ=255)
                                                                             gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                     >  1:15030/1‑132 (MQ=255)
                                                                             gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                     <  2:15030/132‑1 (MQ=255)
                                                                                 aTGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                 <  5:113232/148‑1 (MQ=255)
                                                                                  tGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaaa                                                                <  7:367724/148‑1 (MQ=255)
                                                                                    gTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAgg                                                              <  2:375504/148‑1 (MQ=255)
                                                                                    gTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAg                                                               <  7:31948/147‑1 (MQ=255)
                                                                                     tataGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                 <  2:107516/144‑1 (MQ=255)
                                                                                     tataGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                 >  1:107516/1‑144 (MQ=255)
                                                                                        agagATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATAGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCaa                                                           >  5:21583/1‑147 (MQ=255)
                                                                                         gtgaTTGCCGGTGTTTTTTTTGATCGGGATACATTAAAATCACTAGGTGTAGATCCTAATATTTCGCATTATTTTTCACGTTCAGCCCGTTTTTTGCCAGGTGTATATTCACTGATATTATCAGTAAATGGCGAAAAAAAAGGCAAc                                                          <  3:96844/145‑1 (MQ=255)
                                                                                         gagaTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACa                                                         >  7:99672/1‑148 (MQ=255)
                                                                                         gagaTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCACGTTTTTTGCAAGGTGAATATTCACTGATAGTATCAGTAAATGTCGAAAAAAAATGCAACa                                                         >  3:241810/1‑148 (MQ=255)
                                                                                          agaTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                     <  3:20278/119‑1 (MQ=255)
                                                                                           gaTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCAATTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTAt                                                                                    >  4:20278/1‑119 (MQ=255)
                                                                                             ttGCAGCTGTTGATTTTCATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaa                                                                              <  1:237950/123‑1 (MQ=255)
                                                                                              tGCCGCTGTTGATTTTCATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATTACCCATAAAATTAAAAAAAAAAACCGTTTTTAGCCAGGTGAATATTCACTGATAGTATCAGTaaa                                                                             >  2:237950/1‑123 (MQ=255)
                                                                                               ggcggtGTTGATTTTGATCGGGAAACATTAAAATCACTTGGTGTAGATCATAATATATCGCATTTTTTTTGCCGTTGAGGCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACTTTGCt                                                    <  8:318567/142‑1 (MQ=255)
                                                                                                   cTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCTTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaa                                                                              <  3:92381/117‑1 (MQ=255)
                                                                                                    tGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATATCGCATAAATTTAACAAAAAAACCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaaa                                                                             >  4:92381/1‑117 (MQ=255)
                                                                                                      ttGATTTTGATCGGGAAACTTTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTg                                                                                           <  3:131515/101‑1 (MQ=255)
                                                                                                       tGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGAt                                                                                         <  6:57858/102‑1 (MQ=255)
                                                                                                       tGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGAt                                                                                         >  5:57858/1‑102 (MQ=255)
                                                                                                       tGATATTGATCGCGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTTAAAAAAAAAACCGTTTTTTGCCAGGTGAATATTCACTGa                                                                                          >  4:131515/1‑101 (MQ=255)
                                                                                                         aTTTTGATAGGGAAACATTAAAATCACTAGGTGTAGATACTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTTAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTtga                                         <  1:362874/148‑1 (MQ=255)
                                                                                                          ttttGATCGGGAAACATTAAAATCACTTGGTGTAGATCCTAATATAACGCATTATTTTTCCCGTTCATCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTtga                                         <  1:289530/147‑1 (MQ=255)
                                                                                                          ttttGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTtga                                         <  4:214841/147‑1 (MQ=255)
                                                                                                              gATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGaa                                     <  2:33106/147‑1 (MQ=255)
                                                                                                               attgggAAAAATTAAAATAACTAGGTGTAGATACTAATATATCGAATTATTTTTCACGTTCAGACGGTTTTTTGTCAGGTTAATATTTAATGATAGTATCAGTAAATGGCGAAAAAAAAGGCTATATTGGTATTAGATTTGATGaaa                                    <  4:295174/144‑1 (MQ=255)
                                                                                                               aTCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTg                                                                                           <  7:31966/92‑1 (MQ=255)
                                                                                                               aTCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGaaa                                    <  8:195075/147‑1 (MQ=255)
                                                                                                                tCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAAt                                  <  4:165682/148‑1 (MQ=255)
                                                                                                                tCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATATCCCATTATTTTTAAAATTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGa                                                                                          >  8:31966/1‑92 (MQ=255)
                                                                                                                 cGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACAt                                                        <  6:391104/125‑1 (MQ=255)
                                                                                                                 cGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACAt                                                        >  5:391104/1‑125 (MQ=255)
                                                                                                                                  aCTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGa                 <  6:261435/147‑1 (MQ=255)
                                                                                                                                      ggTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCa              <  5:94001/146‑1 (MQ=255)
                                                                                                                                        tgtAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCa          >  1:387416/1‑148 (MQ=255)
                                                                                                                                         gtAGATCCAAATATAACGCATAATATTTCCCTTAAAGACCGTTTTTTGGCAGGTGAATATTCACTGATAGTATCAGTAATTGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAATTGACTTCTGCAGGCAt         >  8:389538/1‑148 (MQ=255)
                                                                                                                                          tAGATCCTAATTTATCGAATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAtt        <  8:234765/148‑1 (MQ=255)
                                                                                                                                          tAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATAATCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAtt        >  6:231151/1‑148 (MQ=255)
                                                                                                                                           aGATCATAATATATCGCATTATTTTTCACTTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAtt        <  3:252030/147‑1 (MQ=255)
                                                                                                                                            gATCATAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt       <  1:374029/147‑1 (MQ=255)
                                                                                                                                              tCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCTGTAAATTGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCtt    >  1:276515/1‑147 (MQ=255)
                                                                                                                                                cTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTTCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca  <  1:107393/148‑1 (MQ=255)
                                                                                                                                                cTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca  <  6:391978/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCA  >  NC_000913/3362790‑3363081

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: