Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,372,696 C→A 7.0% L158L (CTG→CTT carA ← carbamoyl‑phosphate synthase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,372,6960CA 7.0% 256.1 / 6.4 84L158L (CTG→CTTcarAcarbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base A (4/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

TGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGCGATGG  >  NC_002947/5372548‑5372844
                                                                                                                                                     |                                                                                                                                                     
tGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGcc                                                                                                                                                       <  4:354995/149‑1 (MQ=255)
   gcAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAgg                                                                                                                                                    >  2:288766/1‑149 (MQ=255)
    cAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGt                                                                                                                                                   >  1:351145/1‑149 (MQ=255)
      ggATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTcc                                                                                                                                                 <  4:122922/149‑1 (MQ=255)
           ttCAGCTTGACGCCATAGTCGAAGGCAACCACGTCCTACCGCACCTCGGCAGCGTCGATGGTCGTGAGGCTGTCGGTTTACAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGTAGGAGACGACCTAAGCCAGATCCAtcc                                                                                                                                             >  1:309422/1‑146 (MQ=255)
             cAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGccc                                                                                                                                           <  2:272292/148‑1 (MQ=255)
              aGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCtt                                                                                                                                         <  7:6968/149‑1 (MQ=255)
               gCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCAt                                                                                                                                               <  6:97841/142‑1 (MQ=255)
               gCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCAt                                                                                                                                               >  5:97841/1‑142 (MQ=255)
                cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGATCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                       <  1:290308/149‑1 (MQ=255)
                cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                       <  3:212996/149‑1 (MQ=255)
                cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                       <  7:303579/149‑1 (MQ=255)
                   gACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGt                                                                                                                                                   >  1:85758/1‑134 (MQ=255)
                   gACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGt                                                                                                                                                   <  2:85758/134‑1 (MQ=255)
                      gCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCg                                                                                                                                 >  6:142932/1‑149 (MQ=255)
                      gCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCg                                                                                                                                 <  4:144122/149‑1 (MQ=255)
                        cATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCgg                                                                                                                                >  5:262636/1‑148 (MQ=255)
                         aTAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCggg                                                                                                                               >  3:45631/1‑148 (MQ=255)
                         aTAGTCGAAGGCAACCACGTGGTACGGCAGGTAGGAGGCGTCGATGGTCGGGTGGCAGTCGGTTTTCCGCTCCAACCCGCTGCCGCGCAACGCGAAACGTTACTAGCTGGAGACGACCTTGGCCACGTCCCTGCCCTTCAGGCCCGCGa                                                                                                                              >  4:154853/1‑149 (MQ=255)
                          tAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGaa                                                                                                                             <  2:192603/149‑1 (MQ=255)
                          tAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGaa                                                                                                                             <  7:41769/149‑1 (MQ=255)
                             tCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGcc                                                                                                                          >  8:277319/1‑149 (MQ=255)
                               gAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCcgc                                                                                                                        >  5:81384/1‑149 (MQ=255)
                                            tGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGAt                                                                                                           <  2:178318/149‑1 (MQ=255)
                                              gTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGTCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATg                                                                                                          <  4:284049/148‑1 (MQ=255)
                                                  ggCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGc                                                                                                     <  1:351644/149‑1 (MQ=255)
                                                       gTCGGCAGCGTCGATGGTCGGGTGGCTCTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGcttctt                                                                                                >  3:299581/1‑149 (MQ=255)
                                                                gTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGAtgt                                                                                       >  6:119515/1‑149 (MQ=255)
                                                                   gATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGAtgttgt                                                                                    <  3:152710/149‑1 (MQ=255)
                                                                      ggTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAc                                                                                 >  7:291960/1‑149 (MQ=255)
                                                                       gTCGTGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGGGTAGTCGACCTTGGCGATGTCCATGCCCTTCAGGCCCGTGAAGCCGCGGGCAGCTTCGATGGCAGCTTCTTCGCTGATGTTGTGAc                                                                                 <  3:154853/148‑1 (MQ=255)
                                                                          gggTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGc                                                                             >  3:211272/1‑149 (MQ=255)
                                                                               gCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGaa                                                                        >  2:103878/1‑149 (MQ=255)
                                                                                   tCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATg                                                                      >  6:193467/1‑147 (MQ=255)
                                                                                   tCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCa                                                                   <  4:299581/149‑1 (MQ=255)
                                                                                    cGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACTTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATCCAg                                                                  <  3:219494/149‑1 (MQ=255)
                                                                                    cGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAg                                                                  >  6:235501/1‑149 (MQ=255)
                                                                                     ggTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAg                                                                  <  2:9632/148‑1 (MQ=255)
                                                                                      gTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAccc                                                                               >  5:249469/1‑135 (MQ=255)
                                                                                      gTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAccc                                                                               <  6:249469/135‑1 (MQ=255)
                                                                                      gTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGcc                                                                >  7:91075/1‑149 (MQ=255)
                                                                                       ttttCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCg                                                                                              >  2:276998/1‑119 (MQ=255)
                                                                                       ttttCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCg                                                                                              <  1:276998/119‑1 (MQ=255)
                                                                                           cAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTc                                                            <  2:140095/148‑1 (MQ=255)
                                                                                                 ccACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGgcgc                                                     >  8:292211/1‑149 (MQ=255)
                                                                                                  cacaCGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGcc                                                    >  3:145454/1‑149 (MQ=255)
                                                                                                   acacGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCCGCGATGGCAGCTTCTGCGCTGATGTTGTCACCCGCCAGAATGCCGCGGGTGTGGGGGccc                                                   >  8:318352/1‑149 (MQ=255)
                                                                                                   acacGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGccc                                                   >  3:316267/1‑149 (MQ=255)
                                                                                                       gCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCtttt                                               >  6:275717/1‑149 (MQ=255)
                                                                                                        cTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCTGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAg                                                                  <  4:258311/129‑1 (MQ=255)
                                                                                                        cTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAg                                                                  >  3:258311/1‑129 (MQ=255)
                                                                                                         tGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGgca                                                                                                                   <  6:288964/80‑1 (MQ=255)
                                                                                                         tGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGgca                                                                                                                   >  5:288964/1‑80 (MQ=255)
                                                                                                             gcgcCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCa                                         <  5:235501/149‑1 (MQ=255)
                                                                                                                ccACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGGTGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGa                                      >  7:345407/1‑149 (MQ=255)
                                                                                                                   cTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATAc                                   <  4:211272/149‑1 (MQ=255)
                                                                                                                      gTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACggg                                >  3:310725/1‑149 (MQ=255)
                                                                                                                      gTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACgg                                 <  4:316267/148‑1 (MQ=255)
                                                                                                                       tAACGTTCCTTGGTGGAGACGACCTTGGGCCAGGTCCATCCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGAAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATCAAGCCGTTCTGGGCGCCCTTTTCa                                             >  1:212319/1‑136 (MQ=255)
                                                                                                                       tAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCa                                             <  2:212319/135‑1 (MQ=255)
                                                                                                                       tAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGt                               <  5:193467/149‑1 (MQ=255)
                                                                                                                         aCGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCa                             <  5:275717/149‑1 (MQ=255)
                                                                                                                          cGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                            <  6:81384/149‑1 (MQ=255)
                                                                                                                          cGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                            <  2:123549/149‑1 (MQ=255)
                                                                                                                                  gggggAGACGACCTTGGCCAGGTCCATGCCAGTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACg                      <  6:262636/144‑1 (MQ=255)
                                                                                                                                   gtggAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGgt                   <  1:103878/149‑1 (MQ=255)
                                                                                                                                    tggAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGgtg                  >  5:349694/1‑149 (MQ=255)
                                                                                                                                    tggAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGgtg                  >  1:341599/1‑149 (MQ=255)
                                                                                                                                    gggcatggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGAGGCCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATCCAGCCGTTCAGGGCGGCCTTTTCACGCACGATACGGCTCAGGCGACGGgtg                  >  7:223361/8‑149 (MQ=255)
                                                                                                                                      gagaCGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCATGATACGGGACAGGCGACGGGTGTc                >  6:32117/1‑149 (MQ=255)
                                                                                                                                      gagaCGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTc                <  5:119515/149‑1 (MQ=255)
                                                                                                                                       catggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                            <  6:269013/132‑1 (MQ=255)
                                                                                                                                       catggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                            >  5:269013/5‑136 (MQ=255)
                                                                                                                                       catggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTc                >  7:264380/5‑148 (MQ=255)
                                                                                                                                       agaCGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCg               <  5:142932/149‑1 (MQ=255)
                                                                                                                                        gacgacCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGa              >  4:298800/1‑149 (MQ=255)
                                                                                                                                        atggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTg                                                          <  3:276850/102‑1 (MQ=255)
                                                                                                                                        atggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTg                                                          >  4:276850/4‑105 (MQ=255)
                                                                                                                                         acgacCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGAt             >  3:233223/1‑149 (MQ=255)
                                                                                                                                           gacCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAc                                                                                 >  6:52126/1‑80 (MQ=255)
                                                                                                                                           gacCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAc                                                                                 <  5:52126/80‑1 (MQ=255)
                                                                                                                                           gacCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGCCACCCGCCAGAATGTGAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATg            <  1:288766/149‑1 (MQ=255)
                                                                                                                                              cTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCgg        <  4:45631/149‑1 (MQ=255)
                                                                                                                                              cTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCgg        >  1:197703/1‑149 (MQ=255)
                                                                                                                                                  ggCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGc       <  8:291960/147‑1 (MQ=255)
                                                                                                                                                     cAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGCGATgg  >  6:284742/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGCGATGG  >  NC_002947/5372548‑5372844

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: