Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,002,893 A→G 5.2% G238G (GGT→GGC) ‡ 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,8930AG 5.2% 300.9 / 4.0 96G238G (GGT→GGC) ‡fliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base G (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 = 9.29e-01

AACAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAA  >  NC_000913/2002749‑2003040
                                                                                                                                                |                                                                                                                                                   
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)
          ccGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCGGTGGTGAtaacga                                                                                                                                        <  8:59641/148‑6 (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                                                                                                                 <  2:394990/149‑1 (MQ=255)
                                cTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCa                                                                                                                 <  8:138050/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                                                                                                               <  5:63954/148‑1 (MQ=255)
                                   tagtagTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGt                                                                                                               <  2:73421/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                                                                                 >  5:381056/1‑149 (MQ=255)
                                                                tttGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGt                                                                                 >  4:261623/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                                                          >  4:326701/1‑148 (MQ=255)
                                                                                        gTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCtt                                                          <  3:326701/148‑1 (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                                                                                                                           >  5:470082/1‑55 (MQ=255)
                                                                                                                    cTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCAtt                                                                                                                           <  6:470082/55‑1 (MQ=255)
                                                                                                                     tGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTaa                            >  3:313745/1‑149 (MQ=255)
                                                                                                                       cGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAata                          <  2:392047/149‑1 (MQ=255)
                                                                                                                       cGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAata                          <  8:438938/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                                                                                          <  3:445368/77‑1 (MQ=255)
                                                                                                                               ttCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGc                                                                                          >  4:445368/1‑77 (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             >  3:148792/1‑149 (MQ=255)
                                                                                                                                    aTCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGtggtgg             <  6:381056/149‑1 (MQ=255)
                                                                                                                                    aTCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGtggtgg             <  2:228074/149‑1 (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)
                                                                                                                                                |                                                                                                                                                   
AACAGGTGTACCGCCTGAAGTGATAGTTGTAGCTTTAGTAGTATTTGCATCAGTTACAGTTGCATTTGCCGTTGCTCCAGTCGCCATTGTCACTGTACCATCATTAGCAACTGTTACTGCGTAATACTTCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAA  >  NC_000913/2002749‑2003040

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: