Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,414 T→G 7.0% intergenic (‑143/+308) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,4140TG 7.0% 139.6 / 4.2 63intergenic (‑143/+308)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (24/33);  new base G (3/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

TGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC  >  NC_000913/4094267‑4094562
                                                                                                                                                   |                                                                                                                                                    
tGCCTGTTGTGTGTTCCGGCAAATTTGGGGGTGTCGGATGCACGACCCGTCGGCCGGATAAGGCGCTCGCGCCGCAGCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTg                                                                                                                                                     <  2:216463/149‑1 (MQ=18)
    tGTTCTGTGTTCCTGGAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                                                                                                 <  1:639944/149‑1 (MQ=255)
     gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                                                                                                                <  1:400915/149‑1 (MQ=255)
       tCTGTGTTCCTGTAAATTGGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCGGGCAATCAATGCCTGATGCGACGCTGTCGGGTCTTATCAGGCCTACCACTATTGGCGCCGTGTAGGCCTGATCAGGc                                                                                                                                              >  1:108411/1‑149 (MQ=21)
         tgtgTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCGACCGGCAATCAATGCCTGATGCGACGCTGTCGAGGCTTGACAGGCATACAACTATTGCCGACCTGTAGGACGAATAAgcgcc                                                                                                                                            >  2:41665/1‑145 (MQ=18)
          gtgtTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACt                                                                                                                                           >  1:145794/1‑149 (MQ=255)
           tgtTCCTGTAAATTCGGGGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACtt                                                                                                                                          <  1:655107/149‑1 (MQ=255)
               ccTCTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTCTCGCGTCTTAGCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                                                                                                                      <  2:666294/149‑1 (MQ=255)
                cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                     <  1:228701/149‑1 (MQ=255)
                cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                     >  1:578340/1‑149 (MQ=255)
                 tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGTTTTTACGCTTGCGCCGGCTCCGGCGATCAAAGCTTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCg                                                                                                                                    >  2:517663/1‑149 (MQ=17)
                 tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTTGGCCTGATAAGGCGCTCGTGCCg                                                                                                                                    >  2:565108/1‑149 (MQ=255)
                       ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                                                                                                              >  1:702045/1‑149 (MQ=255)
                            tgttgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCAGGcatt                                                                                                                         <  1:613404/149‑3 (MQ=255)
                                 tCGGATGCACGACCCGTAGGCCGGCTAAGGCGCTCGCGCCGCATCCGGCGATCAGTACCCAATGCAACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAAGCCTGATAAGGCACTTGTGCCGGATCCGGCAATCAATg                                                                                                                    >  1:571162/1‑149 (MQ=17)
                                 ggggAGGCACGACCCGGAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                                    <  2:383835/147‑1 (MQ=255)
                                  cGGATGGACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAACGCCTGATGCGACGCTGTCGCGTCCTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                                                                                                   <  2:309357/149‑1 (MQ=255)
                                     gtGCACGACCCGTAGGCCGGAGAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTg                                                                                                                <  1:302502/148‑1 (MQ=255)
                                         aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGGCATCAATGCCTGATGc                                                                                                            >  1:237883/1‑149 (MQ=255)
                                                gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTg                                                                                                     <  2:106352/149‑1 (MQ=39)
                                                    gtcggtcaaGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                                                 <  2:73664/142‑1 (MQ=255)
                                                     ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                >  1:610804/1‑149 (MQ=25)
                                                      cGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCCTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                               >  1:87085/1‑149 (MQ=255)
                                                       ggATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                              >  1:487099/1‑149 (MQ=255)
                                                        gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                             <  2:459028/149‑1 (MQ=17)
                                                        gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                             <  2:459079/149‑1 (MQ=2)
                                                        gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                             <  1:426102/149‑1 (MQ=39)
                                                        gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCCGATAAGGCACGTGTGCCGCATCCGGCCATCAATGCCTGATGCGACGCTGTCGCGACt                                                                                             >  2:45778/1‑149 (MQ=14)
                                                          tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                           >  1:585612/1‑149 (MQ=255)
                                                           gtggCGCGCGCGCCGCATCCGGGAATGCATGCCTGATGCGACGCTGTCGCGGCTTATCAGGCCTCCAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                          <  2:398687/147‑1 (MQ=11)
                                                           aaGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                          <  2:145794/149‑1 (MQ=40)
                                                            aGGCGCTCGCGCCGCATCCGGCAGTCACTGCCTGATGCGACGCTGTCGCGTCTTATCCGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                         <  2:462703/149‑1 (MQ=255)
                                                             ggCGCTCGCGCCGCATCCGGAAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                                                        <  1:185569/149‑1 (MQ=255)
                                                              gcgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                       <  1:258730/149‑1 (MQ=38)
                                                               cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGCCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                      >  1:370878/1‑149 (MQ=38)
                                                                gcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCCTCCGGCGATCAATGCATGATGGGACGCTGTCGCGTCTTAttcggc                                                                                     >  2:214942/1‑144 (MQ=11)
                                                                             ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGc                                                                        <  1:642292/149‑1 (MQ=32)
                                                                               ggCAATCCATGCCTGATGCGACGCTGTCGCGTCGTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGTCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                      <  1:324258/149‑1 (MQ=12)
                                                                               ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                                         >  2:457178/1‑114 (MQ=32)
                                                                                    tcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAcccg                                                                 >  1:554139/1‑147 (MQ=14)
                                                                                      aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTa                                                               >  2:527049/1‑149 (MQ=34)
                                                                                          ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCGATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCTCCGACCTGTAGGcc                                                           >  2:72132/1‑149 (MQ=255)
                                                                                           cTGGTGCGACGCTGTCGCCTCTTTTCAGGCCTACAACTATCGCCGACCTGCAGGCCTGATAAGGCACTTGTGCCGCCCCCGTCAATCAATGCCCGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCg                                                          <  1:245789/149‑1 (MQ=11)
                                                                                             gATGCGACGCCGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTCTTATCAGGACTACAACGGTTGCTGAACTGTCGGCCGGa                                                        >  1:347018/1‑149 (MQ=11)
                                                                                              aTGCGACGCGGTCGCGTCGTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGAt                                                       <  1:665520/149‑1 (MQ=25)
                                                                                                gcgcCGCGGTCGGGTCTTATCAGGCCTCCAACGATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                     <  2:458315/149‑1 (MQ=11)
                                                                                                 cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                >  2:686561/1‑105 (MQ=33)
                                                                                                 cGACGCGGTCGCGTCTTACCCGGCCTACAACTATGGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAg                                                    <  1:285903/149‑1 (MQ=255)
                                                                                                          cgcgGCTTATCAGGCCTACAACTATGGCCGGCCTGTAGGCCTGATAAGGCACTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgtg                                           <  2:152271/149‑1 (MQ=255)
                                                                                                            cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                         <  2:390356/149‑1 (MQ=38)
                                                                                                                    cAGGCCTACAATTATTGCCGACCTGTAGACCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCgg                                 <  2:242687/149‑1 (MQ=255)
                                                                                                                      ggCCTACAACTTTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGGTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGca                               <  1:233538/149‑1 (MQ=21)
                                                                                                                          tACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATACGGACATCTATTCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca                           >  1:144195/1‑149 (MQ=11)
                                                                                                                           ccAACTATTCCCGACCTGTAGGCCTGATAAGGCCCTTGGGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaa                          <  2:29210/148‑1 (MQ=255)
                                                                                                                           aCAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaa                          <  2:551633/149‑1 (MQ=37)
                                                                                                                            cAACTATTGCCGACCTGTAGGCCTGATAAGGCCATTGTCACGCACCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat                         <  1:565108/149‑1 (MQ=11)
                                                                                                                                    gCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                 <  1:294691/149‑1 (MQ=34)
                                                                                                                                      cGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCg               >  2:130077/1‑149 (MQ=35)
                                                                                                                                       gACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAGGCGGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATCAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGaggaga              >  2:550768/1‑144 (MQ=255)
                                                                                                                                        aCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc             >  2:573880/1‑149 (MQ=34)
                                                                                                                                               ggCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCCCTCGAGCCGCATCCGGCAATCAATGCCTGCTGCGCCGCCGTcg      >  2:222679/1‑149 (MQ=11)
                                                                                                                                                gCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTACGCCGGAAAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc     >  2:485253/1‑149 (MQ=21)
                                                                                                                                                gCCTGATAAGGCACTCGGGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc     <  2:398746/149‑1 (MQ=25)
                                                                                                                                                   tGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc  <  1:527049/149‑1 (MQ=33)
                                                                                                                                                   |                                                                                                                                                    
TGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC  >  NC_000913/4094267‑4094562

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