Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,372,694 G→C 7.2% A159G (GCC→GGC)  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,6940GC 7.2% 261.0 / 9.7 84A159G (GCC→GGC) carAcarbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base C (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 = 9.12e-01

CATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGC  >  NC_002947/5372546‑5372839
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                          <  6:97841/142‑1 (MQ=255)
                 gCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCAt                                                                                                                                          >  5:97841/1‑142 (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                                                                                                                        <  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                                                                                                              >  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                       <  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             >  1:341599/1‑149 (MQ=255)
                                                                                                                                      tggAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGgtg             >  5:349694/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                       >  5:269013/5‑136 (MQ=255)
                                                                                                                                         catggacCTGGCCAAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAg                       <  6:269013/132‑1 (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)
                                                                                                                                                    ggCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGc  <  8:291960/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CATGCGCAGGATGTTCAGCTTGACGCCATAGTCGAAGGCAACCACGTGGTACGGCAGGTCGGCAGCGTCGATGGTCGGGTGGCTGTCGGTTTTCAGCTCCCACACGCTGGAGCGCCACTCGTAACGTTCCTTGGTGGAGACGACCTTGGCCAGGTCCATGCCCTTCAGGCCCGGGAAGCCGCGGGCAGCAGCGATGGCAGCTTCTTCGCTGATGTTGTCACCCGCCAGAATGCAGCCGTTCTGGGCGCCCTTTTCACGCAGGATACGGGTCAGGCGACGGGTGTCGATGCCGGC  >  NC_002947/5372546‑5372839

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: