Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,372,692 T→G 6.8% K160Q (AAG→CAG)  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,6920TG 6.8% 258.1 / 6.1 86K160Q (AAG→CAG) carAcarbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base G (4/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGG  >  NC_002947/5372544‑5372838
                                                                                                                                                    |                                                                                                                                                    
aGCATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCtt                                                                                                                                                      <  3:347626/149‑1 (MQ=255)
 gCATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAGCCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTg                                                                                                                                                    <  2:218713/149‑1 (MQ=255)
  cATGCGCAGGATGTTCAGCTTGACGCGATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTgg                                                                                                                                                   >  1:9632/1‑149 (MQ=255)
    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                                                                                                                                         >  5:97841/1‑142 (MQ=255)
                   gCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCAt                                                                                                                                         <  6:97841/142‑1 (MQ=255)
                    cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGATCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                 <  1:290308/149‑1 (MQ=255)
                    cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                 <  7:303579/149‑1 (MQ=255)
                    cTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCa                                                                                                                                 <  3:212996/149‑1 (MQ=255)
                       gACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGt                                                                                                                                             <  2:85758/134‑1 (MQ=255)
                       gACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGt                                                                                                                                             >  1:85758/1‑134 (MQ=255)
                          gCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCg                                                                                                                           <  4:144122/149‑1 (MQ=255)
                          gCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCg                                                                                                                           >  6:142932/1‑149 (MQ=255)
                            cATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCgg                                                                                                                          >  5:262636/1‑148 (MQ=255)
                             aTAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCggg                                                                                                                         >  3:45631/1‑148 (MQ=255)
                             aTAGTCGAAGGCAACCACGTGGTACGGCAGGTAGGAGGCGTCGATGGTCGGGTGGCAGTCGGTTTTCCGCTCCAACCCGCTGCCGCGCAACGCGAAACGTTACTAGCTGGAGACGACCTTGGCCACGTCCCTGCCCTTCAGGCCCGCGa                                                                                                                        >  4:154853/1‑149 (MQ=255)
                              tAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGaa                                                                                                                       <  7:41769/149‑1 (MQ=255)
                              tAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGaa                                                                                                                       <  2:192603/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                                                                         <  6:249469/135‑1 (MQ=255)
                                                                                          gTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCAccc                                                                         >  5:249469/1‑135 (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                                                                                                             >  5:288964/1‑80 (MQ=255)
                                                                                                             tGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGgca                                                                                                             <  6:288964/80‑1 (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                      <  2:123549/149‑1 (MQ=255)
                                                                                                                              cGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                      <  6:81384/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                                                    >  4:276850/4‑105 (MQ=255)
                                                                                                                                            atggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTg                                                    <  3:276850/102‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
AGCATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGG  >  NC_002947/5372544‑5372838

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: