Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,385,151 T→A 5.3% intergenic (‑121/+40) proS ← / ← oprD proline‑‑tRNA ligase/basic amino acid specific porin OprD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,385,1510TA 5.3% 292.4 / 4.7 94intergenic (‑121/+40)proS/oprDproline‑‑tRNA ligase/basic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base T (48/41);  new base A (3/2);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

GAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACG  >  NC_002947/1385007‑1385293
                                                                                                                                                 |                                                                                                                                              
gAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGcc                                                                                                                                            >  8:373298/1‑149 (MQ=255)
 aGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCg                                                                                                                                           <  6:200429/149‑1 (MQ=255)
 aGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCg                                                                                                                                           <  1:40116/149‑1 (MQ=255)
  gCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCgg                                                                                                                                          <  7:369595/149‑1 (MQ=255)
  gCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCgg                                                                                                                                          <  7:373298/149‑1 (MQ=255)
    aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                        >  6:302981/1‑149 (MQ=255)
    aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                        >  7:366840/1‑149 (MQ=255)
    aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                        >  3:149031/1‑149 (MQ=255)
    aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                        <  6:364912/149‑1 (MQ=255)
         ttGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt                                                                                                                                    <  8:62714/148‑1 (MQ=255)
           gACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTc                                                                                                                                  <  6:321190/148‑1 (MQ=255)
            aCTGGTGCGCATTGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAg                                                                                                                                <  6:267592/149‑1 (MQ=255)
            aCTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTct                                                                                                                                 <  2:238872/148‑2 (MQ=255)
              tGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGt                                                                                                                              >  2:389688/1‑149 (MQ=255)
                gggcgcATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGc                                                                                                                            <  7:147333/147‑1 (MQ=255)
                gTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGc                                                                                                                            >  2:366508/1‑149 (MQ=255)
                 tGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCGGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCg                                                                                                                           >  1:288578/1‑149 (MQ=255)
                    gcATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCgggg                                                                                                                        >  6:188310/1‑149 (MQ=255)
                       tGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCGGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCAc                                                                                                                     <  2:288578/149‑1 (MQ=255)
                         gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTg                                                                                                                   >  7:71185/1‑149 (MQ=255)
                                gtctCGATTGCCTAAGACGTTGAAATAGACCCCCATTGTACGGGGGGGAAACGAATGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAACAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTAAGGGGCAATCCAAGccc                                                                                                            >  4:331364/4‑148 (MQ=255)
                                     gATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCtt                                                                                                       <  8:366840/149‑1 (MQ=255)
                                     gATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCt                                                                                                        >  6:143989/1‑148 (MQ=255)
                                      aTTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTa                                                                                                      >  5:29988/1‑149 (MQ=255)
                                         gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                    <  6:319804/148‑1 (MQ=255)
                                         gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                    <  4:127962/148‑1 (MQ=255)
                                          ccTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAg                                                                                                   >  6:173958/1‑148 (MQ=255)
                                          ccTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAg                                                                                                   <  1:389688/148‑1 (MQ=255)
                                           cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCa                                                                                                            >  8:81232/1‑138 (MQ=255)
                                           cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCa                                                                                                            <  7:81232/138‑1 (MQ=255)
                                           cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa                                                                                                 >  7:87704/1‑149 (MQ=255)
                                           cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa                                                                                                 >  8:29793/1‑149 (MQ=255)
                                             aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCa                                                                                                                      >  5:200840/1‑126 (MQ=255)
                                             aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCa                                                                                                                      <  6:200840/126‑1 (MQ=255)
                                                  gTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAg                                                                                          >  7:219945/1‑149 (MQ=255)
                                                  gTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAg                                                                                          <  5:302981/149‑1 (MQ=255)
                                                  gTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGCCGTACAGCATGCCCCAGGGCGTGAAAATTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGCCTTTTTCAGGTGCGGGGAAATGCAACCCAGGATTACAGGATGCTCAg                                                                                          >  5:292423/1‑149 (MQ=255)
                                                   ttGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                         <  7:29793/149‑1 (MQ=255)
                                                         tAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGTCTTACAGGATGCTCAGCGGGTAc                                                                                   <  6:74093/149‑1 (MQ=255)
                                                         tAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTAc                                                                                   >  1:116039/1‑149 (MQ=255)
                                                              ccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCAc                                                                              >  4:239390/1‑149 (MQ=255)
                                                                        cGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCa                                                                                           >  2:144504/1‑126 (MQ=255)
                                                                        cGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCa                                                                                           <  1:144504/126‑1 (MQ=255)
                                                                             gggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtc                                                               >  6:337954/1‑149 (MQ=255)
                                                                              ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcg                                                              <  3:239390/149‑1 (MQ=255)
                                                                               gAAACGAAGGCGTACATGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgt                                                             <  8:71185/149‑1 (MQ=255)
                                                                               gAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgt                                                             <  4:149031/149‑1 (MQ=255)
                                                                                  aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcga                                                          >  2:119543/1‑149 (MQ=255)
                                                                                   cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgat                                                         >  5:370203/1‑149 (MQ=255)
                                                                                       ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTcgatcgat                                                     <  3:312027/149‑1 (MQ=255)
                                                                                       ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCAAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatcga                                                      <  8:87704/148‑1 (MQ=255)
                                                                                        gCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATg                                                    >  5:340124/1‑149 (MQ=255)
                                                                                          gTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATggg                                                  >  7:35613/1‑149 (MQ=255)
                                                                                          gTACAGGATGCCCCAGGGCGTGACACTTGTCAGCAAAGAAAAAACCCGGCGTGATTGCCGGGCTTTTTCAGGTGCTCGGCACTGCAAGCCAGGCTTACAGCATGCTCAGCGGGTACTCCAAGATCAGGCGCAGTTAGTCGATCGATggg                                                  >  5:246780/1‑149 (MQ=255)
                                                                                            aCAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATgggg                                                 >  5:53833/1‑148 (MQ=255)
                                                                                             cAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGat                                               >  1:307615/1‑148 (MQ=255)
                                                                                                     cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATggg                                                  >  7:54329/1‑138 (MQ=255)
                                                                                                     cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATggg                                                  <  8:54329/138‑1 (MQ=255)
                                                                                                      ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCTGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                         >  4:51932/1‑146 (MQ=255)
                                                                                                      ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                         >  3:50618/1‑146 (MQ=255)
                                                                                                      ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCACGCGCAGTTCGTCGCTCGATGGGGAATCgtag                                         >  5:302624/1‑146 (MQ=255)
                                                                                                      ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCAACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                         >  7:219497/1‑146 (MQ=255)
                                                                                                         gggCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGaca                                   <  6:53833/149‑1 (MQ=255)
                                                                                                            gctGACAAGTGTCACGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGacac                                  >  7:124502/3‑148 (MQ=255)
                                                                                                             gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCAt                               >  3:116398/1‑149 (MQ=255)
                                                                                                               gACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCg                             >  7:79452/1‑149 (MQ=255)
                                                                                                                  acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaa                          <  5:173958/149‑1 (MQ=255)
                                                                                                                  acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaa                          <  2:269468/149‑1 (MQ=255)
                                                                                                                         aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATca                                            <  2:124816/124‑2 (MQ=255)
                                                                                                                         aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATca                                            >  1:124816/1‑123 (MQ=255)
                                                                                                                         aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACgg                    <  8:219945/148‑1 (MQ=255)
                                                                                                                         aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGt                   <  5:188310/149‑1 (MQ=255)
                                                                                                                          cgCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGt                   <  8:124502/148‑1 (MQ=255)
                                                                                                                           gcaaaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTa                  <  6:292423/146‑1 (MQ=255)
                                                                                                                           gCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACgg                    >  1:383348/1‑147 (MQ=255)
                                                                                                                           gCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTa                  <  6:29988/149‑1 (MQ=255)
                                                                                                                             cAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGt                   >  3:89265/1‑146 (MQ=255)
                                                                                                                             cAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGt                   <  4:89265/146‑1 (MQ=255)
                                                                                                                             cAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAgg                <  1:119543/149‑1 (MQ=255)
                                                                                                                             cAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAgg                >  7:197218/1‑149 (MQ=255)
                                                                                                                              ctgAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGt               >  4:134959/3‑149 (MQ=255)
                                                                                                                              ctgAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGt               >  5:144035/3‑149 (MQ=255)
                                                                                                                              aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACg                     <  5:52942/143‑1 (MQ=255)
                                                                                                                              aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACg                     >  6:52942/1‑143 (MQ=255)
                                                                                                                              aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGAAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGAAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGt               >  2:73756/1‑149 (MQ=255)
                                                                                                                               aGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGgcg                                                                     <  8:334818/94‑1 (MQ=255)
                                                                                                                               aGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGgcg                                                                     >  7:334818/1‑94 (MQ=255)
                                                                                                                                gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg             >  8:361369/1‑149 (MQ=255)
                                                                                                                                  aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCTGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGc            <  3:51932/148‑1 (MQ=255)
                                                                                                                                  aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc           <  1:366508/149‑1 (MQ=255)
                                                                                                                                  aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc           <  3:331364/149‑1 (MQ=255)
                                                                                                                                    aaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCt          >  6:255685/1‑148 (MQ=255)
                                                                                                                                        ccGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGacg      >  7:242266/1‑148 (MQ=255)
                                                                                                                                           gCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGacgcacg  >  4:208466/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACG  >  NC_002947/1385007‑1385293

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