Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,122,761 G→A 19.3% intergenic (+32/+115) PP_0983 → / ← PP_0984 membrane protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,122,7610GA19.3% 130.9 / 22.7 63intergenic (+32/+115)PP_0983/PP_0984membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (22/28);  new base A (5/7);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

CAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCA  >  NC_002947/1122613‑1122895
                                                                                                                                                    |                                                                                                                                      
cAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAg                                                                                                                                        >  2:210904/1‑149 (MQ=255)
cAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAg                                                                                                                                        >  5:228947/1‑149 (MQ=255)
      cTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCgg                                                                                                                                  >  3:247716/1‑149 (MQ=255)
         cTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCtt                                                                                                                               >  1:170829/1‑149 (MQ=255)
         cTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtt                                                                                                                               <  8:152510/149‑1 (MQ=255)
         cTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCt                                                                                                                                <  4:327029/148‑1 (MQ=255)
             ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTc                                                                                                                           <  5:70544/149‑1 (MQ=255)
             ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTc                                                                                                                           <  3:256872/149‑1 (MQ=255)
               ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTCGt                                                                                                                         <  2:338961/149‑1 (MQ=255)
                cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTtgtg                                                                                                                        <  6:5025/149‑1 (MQ=38)
                 gTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCTTgg                                                                                                                       <  2:170829/149‑1 (MQ=255)
                 gTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTgg                                                                                                                       <  6:345682/149‑1 (MQ=255)
                 gTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTg                                                                                                                        <  3:156035/148‑1 (MQ=255)
                    gAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGgt                                                                                                                     <  4:204629/148‑1 (MQ=255)
                    gAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtt                                                                                                                               <  8:289067/138‑1 (MQ=255)
                    gAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtt                                                                                                                               <  8:294091/138‑1 (MQ=255)
                    gAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtt                                                                                                                               >  7:289067/1‑138 (MQ=255)
                    gAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtt                                                                                                                               >  7:294091/1‑138 (MQ=255)
                       ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                >  6:198952/1‑134 (MQ=255)
                       ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                <  5:198952/134‑1 (MQ=255)
                         cgagcttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGAt                                                                                                               <  7:128215/144‑1 (MQ=255)
                         aGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTtgtgca                                                                                                                      >  1:216790/1‑140 (MQ=38)
                         aGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTtgtgca                                                                                                                      <  2:216790/142‑3 (MQ=38)
                            tATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTg                                                                                                            >  2:55058/1‑149 (MQ=255)
                            tATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTg                                                                                                            <  3:174161/149‑1 (MQ=255)
                              ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTgg                                                                                                           <  8:172301/148‑1 (MQ=255)
                                cGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGAt                                                                                                               >  8:128215/1‑142 (MQ=255)
                                 ggCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGGGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCg                                                                                                       <  7:295221/149‑1 (MQ=255)
                                   cTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGc                                                                                                     <  1:55058/149‑1 (MQ=255)
                                         gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTc                                                                                               <  1:210904/149‑1 (MQ=255)
                                             ctgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgg                                                                                           <  1:320318/149‑1 (MQ=255)
                                                 tgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGt                                                                                         >  1:122187/1‑147 (MQ=255)
                                                      gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGACCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTAGTGGGTGTTGATCTGGCCAGGCCTCTTCGCGTGTAATCCcg                                                                                  >  3:52853/1‑149 (MQ=255)
                                                         gtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCgcgc                                                                               >  4:341084/1‑149 (MQ=255)
                                                            aTACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTa                                                                            <  6:228947/149‑1 (MQ=255)
                                                             tACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTa                                                                            >  3:5117/1‑148 (MQ=255)
                                                                  gCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTa                                                                      >  2:324530/1‑149 (MQ=255)
                                                                        aTCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACcgcg                                                                 >  8:49938/1‑148 (MQ=255)
                                                                          cTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAg                                                              >  1:19854/1‑149 (MQ=255)
                                                                           tgtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGc                                                                                                     >  7:9133/1‑109 (MQ=255)
                                                                           tgtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGc                                                                                                     <  8:9133/109‑1 (MQ=255)
                                                                             tgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAg                                                              <  6:273828/146‑1 (MQ=255)
                                                                                cTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCAt                                                          >  2:320619/1‑147 (MQ=255)
                                                                                     ccGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGaaa                                                   <  2:19854/149‑1 (MQ=255)
                                                                                        gtgtgtGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAggg                                                >  2:25469/1‑149 (MQ=255)
                                                                                                  gttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGc                                       <  4:5117/148‑1 (MQ=255)
                                                                                                     gcgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACt                                   >  7:116691/1‑149 (MQ=255)
                                                                                                         cGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCaa                               <  3:341084/149‑1 (MQ=255)
                                                                                                              aGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCctc                          <  1:25469/149‑1 (MQ=255)
                                                                                                               ggttAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCctct                         <  4:74371/146‑1 (MQ=255)
                                                                                                                ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGccc                                             <  6:37583/128‑1 (MQ=255)
                                                                                                                  tAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGc                      >  3:318190/1‑149 (MQ=255)
                                                                                                                    aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCt                     >  1:347442/1‑148 (MQ=255)
                                                                                                                      aGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGc                                                                                     >  7:192895/1‑82 (MQ=255)
                                                                                                                      aGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGc                                                                                     <  8:192895/82‑1 (MQ=255)
                                                                                                                            aGCTGCACAAGAAAGCCGACCTCATGGTCGGCTTTTTCGAGCGTGTTGATCTGGCCCGGGCTCTTCGAGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc            <  4:52853/149‑1 (MQ=255)
                                                                                                                            aGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc            <  2:122187/149‑1 (MQ=255)
                                                                                                                             tcTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc            <  7:49938/147‑1 (MQ=255)
                                                                                                                                ccacgaaaaaGCCGACCCTGAGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCt        >  6:66160/8‑149 (MQ=255)
                                                                                                                                  acgaaaaaGCCGACCCTGAGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTAACGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGTCTTATTTCTTCTGCTTc      >  8:284507/6‑149 (MQ=255)
                                                                                                                                      gAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTAATTTCTTCTGCTTCGGCa  >  5:62295/1‑149 (MQ=255)
                                                                                                                                         aGCCGACCTGAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCGCGGCCTCTTGGCTTGTAAGCCCGCTGGTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCGTCTGGCTTATTTcttcttctt       <  2:196027/141‑1 (MQ=255)
                                                                                                                                         aGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGc                      >  8:144517/1‑126 (MQ=255)
                                                                                                                                         aGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGc                      <  7:144517/126‑1 (MQ=255)
                                                                                                                                         aGCCGACCTCAGGGTCGGCTTTTTCGTAGGTGTTGATCTGGCCCGGCCTCTTCGAGGGTAAGCCAGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCtt       >  1:196027/1‑141 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCA  >  NC_002947/1122613‑1122895

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: