Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,362,956 C→A 13.2% S50S (TCC→TCA 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,9560CA13.2% 289.0 / 4.4 105S50S (TCC→TCAyhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base C (39/49);  new base A (14/3);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 6.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCC  >  NC_000913/3362810‑3363102
                                                                                                                                                  |                                                                                                                                                  
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)
                      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)
                               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)
                                        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)
                                                   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                                                                                                                             <  2:222387/113‑1 (MQ=255)
                                                         gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAgg                                                                                                                             >  1:222387/1‑113 (MQ=255)
                                                         gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                                          <  2:15030/132‑1 (MQ=255)
                                                         gCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTa                                                                                                          >  1:15030/1‑132 (MQ=255)
                                                             aTGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                      <  5:113232/148‑1 (MQ=255)
                                                              tGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaaa                                                                                     <  7:367724/148‑1 (MQ=255)
                                                                gTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAg                                                                                    <  7:31948/147‑1 (MQ=255)
                                                                 tataGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                      >  1:107516/1‑144 (MQ=255)
                                                                 tataGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                      <  2:107516/144‑1 (MQ=255)
                                                                    agagATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATAGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCaa                                                                                >  5:21583/1‑147 (MQ=255)
                                                                     gtgaTTGCCGGTGTTTTTTTTGATCGGGATACATTAAAATCACTAGGTGTAGATCCTAATATTTCGCATTATTTTTCACGTTCAGCCCGTTTTTTGCCAGGTGTATATTCACTGATATTATCAGTAAATGGCGAAAAAAAAGGCAAc                                                                               <  3:96844/145‑1 (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                                                                                                              >  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)
                                                                                                              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)
                                                                                                                             tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaa                      >  4:8883/1‑147 (MQ=255)
                                                                                                                               atatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaac                     >  2:294289/1‑147 (MQ=255)
                                                                                                                                    tagCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGcc                <  5:116882/145‑1 (MQ=255)
                                                                                                                                    tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCg               >  2:374944/1‑148 (MQ=255)
                                                                                                                                     cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCgg              >  6:373157/1‑148 (MQ=255)
                                                                                                                                       cATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTCAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTAGATGAAAATGGCTACATTTGTATTGATAAGGAATATCTGCAACAAGCCGGtt            >  5:219634/1‑148 (MQ=255)
                                                                                                                                        attattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGtt            >  7:38224/1‑147 (MQ=255)
                                                                                                                                         ttattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGttt           >  2:66335/1‑147 (MQ=255)
                                                                                                                                         ttattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGAAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATAAGGAATTTCTGCAACAAGCCGGttt           >  3:67003/1‑147 (MQ=255)
                                                                                                                                          tattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGaaa                                                         >  8:395848/1‑100 (MQ=255)
                                                                                                                                          tattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGaaa                                                         <  7:395848/100‑1 (MQ=255)
                                                                                                                                          tattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTa          >  5:69739/1‑147 (MQ=255)
                                                                                                                                          tattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTa          >  2:137561/1‑147 (MQ=255)
                                                                                                                                           attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTGTGCAAAAAGGCGGTTTaa         >  3:82121/1‑147 (MQ=255)
                                                                                                                                             ttttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaaaa      >  4:274502/1‑148 (MQ=255)
                                                                                                                                               ttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                              <  4:134541/106‑1 (MQ=255)
                                                                                                                                               ttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                              >  3:134541/1‑106 (MQ=255)
                                                                                                                                                 cccGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGTAACAAGCCGGTTTAAAAATTc   <  2:387416/147‑1 (MQ=255)
                                                                                                                                                 cccGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTcc  >  5:195595/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCC  >  NC_000913/3362810‑3363102

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: