Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,002,895 C→T 5.2% G238S (GGT→AGT) ‡ fliC ← flagellar filament structural protein (flagellin)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,002,8950CT 5.2% 308.3 / 4.5 96G238S (GGT→AGT) ‡fliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base C (46/45);  new base T (2/3);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 101 total aligned reads displayed.

AACAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAG  >  NC_000913/2002749‑2003041
                                                                                                                                                  |                                                                                                                                                  
aaCAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCgg                                                                                                                                                  <  6:218284/149‑1 (MQ=255)
aaCAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCg                                                                                                                                                   <  7:186395/148‑1 (MQ=255)
        tACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTc                                                                                                                                          >  6:454330/1‑149 (MQ=255)
        tACCGCCTGAAGTGATAGTGGTAGCTTTAGAAGTATTTGCATCAGTTACAGTTGCATTTGCCGGTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTc                                                                                                                                          <  4:19055/149‑1 (MQ=255)
          ccGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCg                                                                                                                                         >  4:108174/1‑148 (MQ=255)
          ccGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGc                                                                                                                                        <  1:293438/149‑1 (MQ=255)
           cGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCa                                                                                                                                       <  5:454330/149‑1 (MQ=255)
            gCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCa                                                                                                                                       >  2:390171/1‑148 (MQ=255)
             ccTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATa                                                                                                                                     <  6:342531/149‑1 (MQ=255)
                   gTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAAt                                                                                                                                >  1:421367/1‑148 (MQ=255)
                       tAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTa                                                                                                                           <  2:296303/149‑1 (MQ=255)
                         gTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTAcc                                                                                                                         <  3:188534/149‑1 (MQ=255)
                           tGTACCTTTATTATTTATTGCAACAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCa                                                                                                                        <  7:121460/148‑1 (MQ=255)
                            gTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCAt                                                                                                                       >  4:460166/1‑148 (MQ=255)
                               gCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTa                                                                                                                           >  5:96284/1‑141 (MQ=255)
                               gCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTa                                                                                                                           <  6:96284/141‑1 (MQ=255)
                                cTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCa                                                                                                                  <  8:138050/149‑1 (MQ=255)
                                cTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCa                                                                                                                  <  2:394990/149‑1 (MQ=255)
                                 tttAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAg                                                                                                                 >  6:63954/1‑149 (MQ=255)
                                  ttAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATa                                                                                                                                     >  3:171469/1‑128 (MQ=255)
                                  ttAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATa                                                                                                                                     <  4:171469/128‑1 (MQ=255)
                                   tagtagTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGt                                                                                                                <  2:73421/148‑1 (MQ=255)
                                   tagtagTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGt                                                                                                                <  5:63954/148‑1 (MQ=255)
                                      tagtaTTTGCCTCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGtata                                                                                                             >  4:132594/1‑148 (MQ=255)
                                           tttGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACAcc                                                                                                       <  7:448637/149‑1 (MQ=255)
                                                         gcagccactgatactggcggaactaacccagcttcaattgagggtgtttatactgataatggtaatgattactatgcgaaaATCACCGGTGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGa                                                                                              >  6:361322‑M1/82‑144 (MQ=255)
                                                         gcagccactgatactggcggaactaacccagcttcaattgagggtgtttatactgataatggtaatgattactatgcgaaaATCACCGGTGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGa                                                                                              <  5:361322‑M1/63‑1 (MQ=255)
                                                           ttGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTGCCATCGTTATCACCACCTGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGTAGCTggg                                                                                       <  1:390171/149‑1 (MQ=255)
                                                           ttGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTggg                                                                                       >  1:392047/1‑149 (MQ=255)
                                                           ttGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTggg                                                                                       >  4:165021/1‑149 (MQ=255)
                                                              cATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCa                                                                                                                  <  7:118451/119‑1 (MQ=255)
                                                              cATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCa                                                                                                                  >  8:118451/1‑119 (MQ=255)
                                                                tttGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGt                                                                                  >  4:261623/1‑149 (MQ=255)
                                                                tttGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGt                                                                                  >  5:381056/1‑149 (MQ=255)
                                                                    ccGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTcc                                                                               <  2:421367/148‑1 (MQ=255)
                                                                    ccGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCg                                                                              >  7:344171/1‑149 (MQ=255)
                                                                     cGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGc                                                                             <  3:460166/149‑1 (MQ=255)
                                                                      gTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTAAGTTCCGc                                                                             >  4:420374/1‑148 (MQ=255)
                                                                       ttGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGaa                                                                                             >  3:251691/1‑131 (MQ=255)
                                                                       ttGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGaa                                                                                             <  4:251691/131‑1 (MQ=255)
                                                                       ttGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAAGTGaa                                                                                             >  3:257945/1‑131 (MQ=255)
                                                                           tCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTAt                                                                       >  7:48508/1‑149 (MQ=255)
                                                                           tCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTAt                                                                       >  5:368592/1‑149 (MQ=255)
                                                                             cAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCa                                                                     <  3:108174/149‑1 (MQ=255)
                                                                                        gTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCtt                                                           <  3:326701/148‑1 (MQ=255)
                                                                                        gTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCtt                                                           >  4:326701/1‑148 (MQ=255)
                                                                                          cACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCg                                                        >  5:20945/1‑149 (MQ=255)
                                                                                           aCTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCg                                                        >  6:149012/1‑148 (MQ=255)
                                                                                              gTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATc                                                                                                                               <  8:13698/74‑1 (MQ=255)
                                                                                              gTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATc                                                                                                                               >  7:13698/1‑74 (MQ=255)
                                                                                              gTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAg                                                     >  7:438938/1‑148 (MQ=255)
                                                                                                 ccATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGaaa                                                  <  6:368592/148‑1 (MQ=255)
                                                                                                 ccATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAg                                                 >  3:372585/1‑149 (MQ=255)
                                                                                                  catcatTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAgg                                                >  1:228074/1‑149 (MQ=255)
                                                                                                   atcatTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGTGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAggg                                               >  3:351516/1‑149 (MQ=255)
                                                                                                   atcatTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAggg                                               >  8:108073/1‑149 (MQ=255)
                                                                                                    tcatTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGt                                              <  8:344171/149‑1 (MQ=255)
                                                                                                    tcatTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGt                                              <  3:261623/149‑1 (MQ=255)
                                                                                                       ttAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAAt                                           >  8:12744/1‑149 (MQ=255)
                                                                                                          gCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTc                                         >  7:154351/1‑148 (MQ=255)
                                                                                                            actgataatggtaatgattactatgcgaaaATCACCGGTGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAg                                      <  8:178182‑M1/119‑1 (MQ=255)
                                                                                                              cTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGt                                     >  6:403080/1‑148 (MQ=255)
                                                                                                                gTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAg                                  >  4:252570/1‑149 (MQ=255)
                                                                                                                 ttACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAg                                  <  5:403080/148‑1 (MQ=255)
                                                                                                                    cTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCAtt                                                                                                                            <  6:470082/55‑1 (MQ=255)
                                                                                                                    cTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCAtt                                                                                                                            >  5:470082/1‑55 (MQ=255)
                                                                                                                     tGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTaa                             >  3:313745/1‑149 (MQ=255)
                                                                                                                       cGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAata                           <  8:438938/149‑1 (MQ=255)
                                                                                                                       cGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAata                           <  2:392047/149‑1 (MQ=255)
                                                                                                                        gTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAGGTTTAatat                          >  6:145651/1‑149 (MQ=255)
                                                                                                                         tAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAtt                         <  6:20945/149‑1 (MQ=255)
                                                                                                                          aaTACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGGATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGGTTAATAttg                        >  3:316094/1‑149 (MQ=255)
                                                                                                                           aTACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgt                       <  6:165225/149‑1 (MQ=255)
                                                                                                                           aTACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgt                       <  3:252570/149‑1 (MQ=255)
                                                                                                                            tACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgtt                      <  4:372585/149‑1 (MQ=255)
                                                                                                                            attactatgcgaaaATCACCGGTGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgt                       >  7:178182‑M1/15‑148 (MQ=255)
                                                                                                                               ttCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGc                                                                                           >  4:445368/1‑77 (MQ=255)
                                                                                                                               ttCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGc                                                                                           <  3:445368/77‑1 (MQ=255)
                                                                                                                               ttCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgtttgt                   >  5:485195/1‑149 (MQ=255)
                                                                                                                               ttCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgtttgt                   >  7:301044/1‑149 (MQ=255)
                                                                                                                                 cccATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTtgtg                  <  8:48508/148‑1 (MQ=255)
                                                                                                                                   cATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACCCCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTgtggtg               >  3:22984/1‑149 (MQ=255)
                                                                                                                                   cATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTgtggtg               >  5:291533/1‑149 (MQ=255)
                                                                                                                                    aTCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGtggtgg              <  6:381056/149‑1 (MQ=255)
                                                                                                                                    aTCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGtggtgg              <  2:228074/149‑1 (MQ=255)
                                                                                                                                    aTCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGtggtgg              >  3:148792/1‑149 (MQ=255)
                                                                                                                                       gTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTa            <  8:154351/148‑1 (MQ=255)
                                                                                                                                        ttATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGTGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGc          <  4:351516/149‑1 (MQ=255)
                                                                                                                                         tATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCa         >  6:326253/1‑149 (MQ=255)
                                                                                                                                           TCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCAcc       >  4:407623‑M1/1‑149 (MQ=255)
                                                                                                                                           TCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTACCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCAc        <  7:67032‑M1/148‑1 (MQ=255)
                                                                                                                                             ACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCa      <  8:380306‑M1/148‑1 (MQ=255)
                                                                                                                                               CACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaaa   >  3:132594‑M1/1‑149 (MQ=255)
                                                                                                                                               CACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaa     <  6:323320‑M1/147‑1 (MQ=255)
                                                                                                                                               CACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaa     <  2:113623‑M1/147‑1 (MQ=255)
                                                                                                                                                  cGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgtt                      <  3:395230/127‑1 (MQ=255)
                                                                                                                                                  cGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATAttgtt                      >  4:395230/1‑127 (MQ=255)
                                                                                                                                                  cGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaaa   >  3:275637/1‑146 (MQ=255)
                                                                                                                                                  cGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaaa   <  4:275637/146‑1 (MQ=255)
                                                                                                                                                  cGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCAGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAg  >  1:230446/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAG  >  NC_000913/2002749‑2003041

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