Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,423,575 C→G 5.3% P564A (CCA→GCA)  PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,423,5750CG 5.3% 241.6 / 4.1 75P564A (CCA→GCA) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base G (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCGTTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAT  >  NC_002947/3423428‑3423723
                                                                                                                                                   |                                                                                                                                                    
gcgTTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCa                                                                                                                                                    >  2:1046566/1‑150 (MQ=255)
gcgTTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCa                                                                                                                                                    >  1:1261411/1‑150 (MQ=255)
 cgTTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCAc                                                                                                                                                   <  2:70995/150‑1 (MQ=255)
   ttCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGt                                                                                                                                                 >  1:972418/1‑150 (MQ=255)
     ccAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCa                                                                                                                                               >  1:757579/1‑150 (MQ=255)
      cAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCAGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAg                                                                                                                                              >  1:289208/1‑150 (MQ=255)
       aaGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAgt                                                                                                                                             >  2:1142800/1‑150 (MQ=255)
         gCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAgt                                                                                                                                             >  1:711398/1‑148 (MQ=255)
         gCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAgt                                                                                                                                             <  2:711398/148‑1 (MQ=255)
         gCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGACTGACGTGGCACGTCactcg                                                                                                                                           <  1:963343/150‑5 (MQ=255)
            aTCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGtt                                                                                                                                        <  1:624384/150‑1 (MQ=255)
                gAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCTGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTcg                                                                                                                                    <  2:346867/150‑1 (MQ=255)
                 aGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGAAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTcgc                                                                                                                                   >  2:1021856/1‑150 (MQ=255)
                   gTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCgg                                                                                                                                 <  2:283528/150‑1 (MQ=255)
                   gTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCgg                                                                                                                                 >  1:283528/1‑150 (MQ=255)
                       gTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTac                                                                                                                             >  2:457608/1‑150 (MQ=255)
                          gATCGTCGAGCTGGGCGAATTGGACAGCTTCAAGAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAGAGCTATGGGCGAAGAACGAGTGACGTGCCACGCCAGTGCGTTTTGGAGGCTacaac                                                                                                                          >  2:324716/1‑150 (MQ=255)
                              gtcgAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGACATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCaccaac                                                                                                                      >  2:70547/1‑150 (MQ=255)
                                cgAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCaccaacca                                                                                                                    <  1:1142800/150‑1 (MQ=255)
                                    cTGGGCGAATTTGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGAt                                                                                                                <  1:1161556/150‑1 (MQ=255)
                                      gggCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGa                                                                                                              <  2:757579/150‑1 (MQ=255)
                                           aaTTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACNTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATAcc                                                                                                         <  1:808950/150‑1 (MQ=255)
                                              tGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTa                                                                                                      <  2:928907/150‑1 (MQ=255)
                                                gACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAg                                                                                                    <  1:915589/150‑1 (MQ=255)
                                                gACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACAAACCAGGATGAATACCTTAAg                                                                                                    >  1:563840/1‑150 (MQ=255)
                                                     cTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGGGTGTTTTCGCGGGTACCACCAACCAGAAAGAATAGCTTAAGGACGc                                                                                               >  1:290507/1‑150 (MQ=255)
                                                      ttCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGcc                                                                                              >  2:794491/1‑150 (MQ=255)
                                                         aacaaGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAg                                                                                                    >  2:647441/1‑141 (MQ=255)
                                                         aacaaGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAg                                                                                                    <  1:647441/141‑1 (MQ=255)
                                                            aaGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGt                                                                                        >  2:722312/1‑150 (MQ=255)
                                                                   agagCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGCCAGCTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAAt                                                                                                            <  1:1043800/123‑1 (MQ=255)
                                                                   agagCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAAt                                                                                                            >  2:1043800/1‑123 (MQ=255)
                                                                    gagCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCg                                                                                >  1:129103/1‑150 (MQ=255)
                                                                          cAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGt                                                                                                                                                 >  1:720440/1‑79 (MQ=255)
                                                                          cAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGt                                                                                                                                                 <  2:720440/79‑1 (MQ=255)
                                                                          cAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTg                                                                          >  1:659084/1‑150 (MQ=255)
                                                                            aGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGc                                                                        <  1:457608/150‑1 (MQ=255)
                                                                                cAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGt                                                                    <  2:129103/150‑1 (MQ=255)
                                                                                 aaGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAgg                                                                                                   >  2:559178/1‑118 (MQ=255)
                                                                                 aaGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAgg                                                                                                   >  2:559573/1‑118 (MQ=255)
                                                                                 aaGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAgg                                                                                                   <  1:559573/118‑1 (MQ=255)
                                                                                 aaGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAgg                                                                                                   <  1:559178/118‑1 (MQ=255)
                                                                                    cAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTg                                                                          >  2:1133925/1‑140 (MQ=255)
                                                                                    cAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTg                                                                          <  1:1133925/140‑1 (MQ=255)
                                                                                          ttctCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTc                                                          >  1:153775/1‑150 (MQ=255)
                                                                                             tcGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGGGGGGGGCATCTGCCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAg                                                       >  2:744409/1‑150 (MQ=255)
                                                                                                  aTCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGa                                                  <  2:289208/150‑1 (MQ=255)
                                                                                                   tCTATCGACACTTACCGCGAAAGATATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGAc                                                 >  1:33923/1‑150 (MQ=255)
                                                                                                      aTCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCtt                                              >  2:230415/1‑150 (MQ=255)
                                                                                                         gACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGaa                                           <  2:659084/150‑1 (MQ=255)
                                                                                                         gACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGaa                                           <  1:1046566/150‑1 (MQ=255)
                                                                                                         gACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACAAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGaa                                           <  2:563840/150‑1 (MQ=255)
                                                                                                          acacTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAg                                          >  1:158321/1‑150 (MQ=255)
                                                                                                          acacTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAg                                          >  1:1151601/1‑150 (MQ=255)
                                                                                                            acTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTaa                                                                                                     <  1:1068816/89‑1 (MQ=255)
                                                                                                            acTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTaa                                                                                                     >  2:1068816/1‑89 (MQ=255)
                                                                                                            acTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCa                                        <  1:794491/150‑1 (MQ=255)
                                                                                                             cTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCAt                                       <  1:1109446/150‑1 (MQ=255)
                                                                                                              ttACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTNANCGGCGTNATTGGCCGGTCGNTTGCGTCAAGGTGGACCTTGAAGCAtt                                      <  1:722312/150‑1 (MQ=255)
                                                                                                                   gcgAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGAc                                                                                                 <  2:554576/86‑1 (MQ=255)
                                                                                                                   gcgAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGAc                                                                                                 >  1:554576/1‑86 (MQ=255)
                                                                                                                           tATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAg                                          <  2:101890/133‑1 (MQ=255)
                                                                                                                           tATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAg                                          >  1:101890/1‑133 (MQ=255)
                                                                                                                             tGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTGAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGa                       >  2:686130/1‑150 (MQ=255)
                                                                                                                             tGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGa                       <  1:70547/150‑1 (MQ=255)
                                                                                                                                      aaCGAGTGACGTGCCACGGCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCtgtg              <  2:972418/150‑1 (MQ=255)
                                                                                                                                       aCGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGTTGAAGACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTgg             <  2:33923/150‑1 (MQ=255)
                                                                                                                                        acactGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTAtt                                                                           <  2:616190/83‑1 (MQ=255)
                                                                                                                                        acactGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTAtt                                                                           >  1:616190/5‑87 (MQ=255)
                                                                                                                                         gAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGc           <  1:1021856/150‑1 (MQ=255)
                                                                                                                                          actGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCt          >  2:950849/3‑150 (MQ=255)
                                                                                                                                            tGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGa        >  2:246289/1‑150 (MQ=255)
                                                                                                                                             gACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAg       <  2:1261411/150‑1 (MQ=255)
                                                                                                                                              aCGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAgg      >  1:492169/1‑150 (MQ=255)
                                                                                                                                                  gCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAt  <  2:158321/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAT  >  NC_002947/3423428‑3423723

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