Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.7% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC17.7% 347.4 / 40.4 219intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (63/114);  new base C (6/36);  total (69/150)
Fisher's exact test for biased strand distribution p-value = 9.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                     |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  1:9527613/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                              >  1:811484/1‑150 (MQ=255)
         gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  1:8366361/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  1:7128917/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgg                                                                                                                                         >  1:5179780/1‑149 (MQ=255)
                gataTGAATGGGCGTTATCCCAAACATGTCTGGCCGGACGCCCCGGCAAATACGGCCCCGACGCGAGGGACTAAAAATTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtca                                                                                                                                      <  2:758786/147‑1 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGtggcc                                                                                                                                     >  1:8966587/1‑146 (MQ=255)
                  gaTGAATGGGCGTTATCCCAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTGGTATGGTTGGTTTTCTCCCTCGCCCTGGGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:9591041/150‑1 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccca                                                                                                                                  >  2:3251406/1‑149 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTggaccccc                                                                                                                                 >  1:7235657/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccccc                                                                                                                                 >  1:4988437/1‑148 (MQ=255)
                      aacgggCTTTATCCCAAACTTGTCGGGCCGCACGACCGGTAAATTACTGCACCGCGGCTAGTTATTAAAATTTTTTGGTATGGTTGTTTTTCTCCCGCTCCGGGGGGGAGTGGGCGGGGGTTTGGGCCTCAGCGCGCACGTtcaccctca                                                                                                                                <  2:7901307/147‑1 (MQ=255)
                       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                               <  2:7373564/150‑1 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacccacaccca                                                                                                                            >  1:8204305/1‑149 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGGCCCacacccacacccaca                                                                                                                    >  1:9688402/1‑147 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCAAACCcact                                                                                                                    >  1:7561753/1‑150 (MQ=255)
                                   ccAACCTTGTCTGGCCGGCCGCCCCGGCAAATACTGCACCGACGAAACGGACTAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:8602040/150‑1 (MQ=255)
                                   ccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCTAAGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTATCCCTCTCCCTGTGGGAGAGGGCCGCGGTGAGGGCATCAGCGCGCACGTTCACCAACACCCCAAcacacta                                                                                                                   >  1:7883847/1‑144 (MQ=37)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCGAACCCCCTcc                                                                                                                  >  2:9015910/1‑150 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAcccacaacc                                                                                                                 >  2:5475521/1‑144 (MQ=38)
                                         ttGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTATGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  1:9029079/149‑1 (MQ=38)
                                         aTGTCTGGCCGTACCCCCCGTAAAATACTCCACCGACGCGACGTACTAAAATTTTTTCGTATGGTTGTTTTCCTCCATCTCCCTGGGGGTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  1:6548703/150‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACTCCCACaa                                                                                                             >  1:9155738/1‑150 (MQ=255)
                                         aTGTCGGGCGGGACGACCGGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTACGTTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCttga                                                                                                             <  2:5297476/150‑4 (MQ=37)
                                          tgtgtgGCCGGACGACCCGCAAATTCCTCCACCGACGCCAGTCACTAAAAATTATACGTAAGTTTGTTTTTCACCCTCACCCTTGGGTTGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:346953/150‑1 (MQ=255)
                                          tGTCTGGCCGGACGTCCCGGAAATGACTGAACGGCCGCGATGGTCTAAAAAGTATTGGTAAGGTTGGTTTTGTGCCTCTCCCTGGGGGTGATTGCGGGGGTGTGGTCCTCAGCGATGACGTTCACCCTCACCCTAACCCTCTCCCTCatg                                                                                                            <  2:6150205/150‑3 (MQ=255)
                                          gggcTGGCGGGACGACCCGGCAAATCCGGCACGGACGCGCCGGACTAAAATTTATTCTTAGGGTGGGTTTACCCCCTCGCCCGGTGGTAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:2073354/147‑1 (MQ=255)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGAGCCCGCACCCCAACCCACTCCCACAAcggg                                                                                                         >  1:1206537/1‑149 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCACTACCACCAggggg                                                                                                        >  2:6258358/1‑148 (MQ=37)
                                              gggCCGGACCACCCGGAAAATAGTGCCCCTCCGCGACGTACTAAAAATTATTCGTTGGTTTGTTTTTCTCCCTCCCGCGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                        <  1:3149440/149‑1 (MQ=255)
                                               ggCCGGACGCCCCGGAAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                       <  1:9089758/150‑1 (MQ=255)
                                               gcccGCACCACCCCGCAAATACAGCACCGACCCGACGCACCAAAAACTATTCGTAAGGTTGTTTTACACCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                       <  2:9878060/148‑1 (MQ=37)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAAACACACCCACAAGGGagag                                                                                                      >  2:8344243/1‑150 (MQ=38)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCGAAACCACACCCCCAAGGgaggg                                                                                                      >  1:5044886/1‑148 (MQ=37)
                                                gCCGGACGACCCGGCAAATACCGCCCCGACGCGAGGGCCGAAAAAGTATTCGGATGGTTGGTTTTCGCCCTCCCCGGGGGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:8745086/150‑1 (MQ=255)
                                                gCCGCACGACCCGGAAAATACTGCACCGACGCGACGCACGAAAATTTTTTCGAAAGTTTGGTTTTCTCCCCCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                      <  2:8057976/150‑1 (MQ=37)
                                                 ccgcccgacCCGCAAAATACTGCACCGACGCGACGGACCAAAAATTTTTCGTATGGTTGTTTTTCCCCCGCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  1:4182843/145‑1 (MQ=37)
                                                 ccGGCCCACCCGCAAAATCCGGCCCCGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTCCCCCCCCCCCGTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  1:8866388/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCACTCCCACCACGGGGAgg                                                                                                     >  1:6678404/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAACCGCACCCACAAGGGGGAggg                                                                                                    >  2:2452351/1‑150 (MQ=37)
                                                     ccgccccGGCACATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTGGTTGTGCGCCCGCGCCGTGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                 <  2:7800579/146‑1 (MQ=255)
                                                      cgccccGGCAAATACGGCACCGCCGCGCCGGACGAAAAAGTATTCGTACGGTGGGTTTTCTCCCTCGCCCTGGGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                <  2:1620714/147‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:1522724/150‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCAAACCCCCACCCACCAAGGGGAGGGGGCCgcgc                                                                                            >  1:4725496/1‑147 (MQ=25)
                                                          ccGGCAAACACTCCACCGACGCGACGCACGAAAAAGTATTCTTACGTTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:2893290/150‑1 (MQ=37)
                                                            ggCAACTACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:6541638/150‑1 (MQ=255)
                                                                      gccccGACGCGCGTGACGAAAAATTATGCGTAAGGTTGTTTTTCTCCCTCTCCGTGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                <  1:4038195/147‑1 (MQ=255)
                                                                      gCACCGACGCGACGCACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGTGTGAGGGCCGGGGTGAGGGCATCAGCGCTCACGTTCACCCTCACCCCAACACTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                <  1:4038272/150‑1 (MQ=37)
                                                                        cccgccgGGACGGACTAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                              <  1:7339851/149‑1 (MQ=255)
                                                                         cggcggcgaggGACGAAAATGTTTTCTTTTGGTTGTTTTTCTCCCTCTCCGGGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:1223959/145‑1 (MQ=255)
                                                                         ccgcggcgccgGACGAAAAATTATTCGTATGGTGGGTTTTCCCCCGCCCCGTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:8896777/145‑1 (MQ=255)
                                                                         ccgccgCGACGTATTAAAATGTATTCTTCTGGTTGGTTTTCCCCCTCGCCGGGTGGGTTTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:1087648/150‑1 (MQ=255)
                                                                         ccGACGCGCGGGACTAAAATTTTTTGGTTTGTTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:9527613/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAATTTTTTCTTAAGGTTGGTTTTCTCCCACTCCCTGTGGTAGCGGGCCGGGTTGAGGGAATCAGCGCGCACGTTCACCCTCACCATAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:6988797/150‑1 (MQ=255)
                                                                         ccGACGCGACGGAGGAAAACTTATTCTTATGTTTGTTTTTCTCCCTATCCCTGTGGGAGTGGGCCGGGTGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:1981147/150‑1 (MQ=255)
                                                                         ccGACGCCACGTACCAAAACGTTTTCTTACGGTTGTCTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:4295399/150‑1 (MQ=255)
                                                                           gacgcgccgGCCTAAAAATTATTCGTATGGTTGGTTTGCGCCCTCTCCCTGTGGGCGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  1:8473619/150‑1 (MQ=255)
                                                                              gcgaggGACTAAAATGTTTTGGTTTGGTGGGTTTTCTCCCCCTCCGGGGGGGCGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                        <  2:3715400/150‑1 (MQ=255)
                                                                                gagggacgAAAAAGTTTTCGTATGTTTGGTTTACTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGGCCGATCGAGCACAAATTTTGAGAGatat                                                                      <  1:4880652/147‑1 (MQ=255)
                                                                                        aaaaaTTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                              <  1:2457573/150‑1 (MQ=255)
                                                                                             ttATTGGTATGGTTGTTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:5544567/149‑1 (MQ=39)
                                                                                             gttttCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:6041844/147‑1 (MQ=255)
                                                                                             gTATTCTCACGGTGGGTGTCCTCCCGCGCCCGGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:5219806/150‑1 (MQ=255)
                                                                                             gTATGCGTACGGTTGGTTTTCTCCCTCGCCGTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:663866/150‑1 (MQ=255)
                                                                                              ttttCGTACGTTTGTTTTTCTCCCTCTCCCTGTGGTTTTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  2:2277982/148‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  2:7407727/150‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCCCACCCAAACCCACGCCCACCAGGGCGAGGGGGCCGATCGCGCCCAAAAGTTGGGTGGTATCTTCTTCTGTCTTgggaccg                                                  >  2:2477422/1‑144 (MQ=16)
                                                                                                     aaGTTTGGTTTTATCCCTCTCCCTGTGGGAGAGGGCCGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  2:4090160/150‑1 (MQ=39)
                                                                                                      tggttggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  1:9934610/149‑1 (MQ=255)
                                                                                                           ggttTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                           <  1:5699035/150‑1 (MQ=255)
                                                                                                           ggttTCCTCCCCCGCCGGGGGGGTGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                           <  1:532533/150‑1 (MQ=255)
                                                                                                                  tccctcgccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                    <  2:3328099/150‑1 (MQ=255)
                                                                                                                  gccctctccgtGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                    <  1:3753142/149‑1 (MQ=255)
                                                                                                                   ccctctcccgGTGGTCGTGGGCCGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                   <  1:8439553/150‑1 (MQ=255)
                                                                                                                    catctccggGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                  <  1:3996844/148‑1 (MQ=255)
                                                                                                                     ctcgcccgGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                 <  2:5796363/150‑1 (MQ=255)
                                                                                                                          cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                            <  2:7014721/150‑1 (MQ=255)
                                                                                                                            tgtgGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          <  1:4051376/150‑1 (MQ=255)
                                                                                                                            ggtgGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          <  1:9742327/149‑1 (MQ=255)
                                                                                                                            ggtgGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          <  2:7424935/149‑1 (MQ=255)
                                                                                                                            ggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          <  2:723858/149‑1 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCCCACCCCAACCCACACCCTCAAGGGAGAGGGGCCCCAACCAGCCCAAAATTTGAGCGGTATATTCTTAGGTGTTGTGACCGAAGAACAGCAAAAAGGGCGGTTGGGcc                       >  2:4110753/1‑150 (MQ=255)
                                                                                                                                ggAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  1:7685612/150‑1 (MQ=255)
                                                                                                                                  tgcgggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                    <  1:2039896/147‑1 (MQ=255)
                                                                                                                                  agtgggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                    <  1:9488773/147‑1 (MQ=255)
                                                                                                                                   gtgggCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                   <  1:3604377/148‑1 (MQ=255)
                                                                                                                                     gggCCGGGGCGTGGGCATCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATGTTGAGAGTTATCTTCTTCTGGCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                 <  2:6627047/150‑1 (MQ=255)
                                                                                                                                      ggCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                <  1:4400507/150‑1 (MQ=255)
                                                                                                                                       gCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:3412935/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:5044886/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:19022/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:1586247/150‑1 (MQ=255)
                                                                                                                                            ggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:2853983/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  2:9155588/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:6302630/150‑1 (MQ=255)
                                                                                                                                             gTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag         <  2:1844673/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:4238710/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:5244797/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCACACCCCAAACCACAACCACAACGGCCAGGGGGACCAACGAGCACAAAATTTGAGAGATAAAATAATATGTCTTTTAAAAGAACAAAAAAAAAACCGGCGTTTTCGCCGGAAAAATTCGGAGaaaaa    >  2:7314424/1‑150 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:4683580/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729

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