Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,329,641 C→A 6.7% L92M (CTG→ATG)  PP_4690 → Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,329,6410CA 6.7% 232.0 / 3.1 72L92M (CTG→ATG) PP_4690Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base A (2/3);  total (36/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 = 1.00e+00

GCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGGCAGCGCGCAGCGCCTGG  >  NC_002947/5329493‑5329789
                                                                                                                                                    |                                                                                                                                                     
gCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGcccc                                                                                                                                                       >  2:264629/1‑149 (MQ=255)
  aCCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTg                                                                                                                                                     <  1:246215/149‑1 (MQ=255)
   ccGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGTAGCCCCTgg                                                                                                                                                    <  5:58577/149‑1 (MQ=255)
    cGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGTAGCCCATggg                                                                                                                                                   >  6:225092/1‑149 (MQ=255)
    cGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTggg                                                                                                                                                   >  8:36032/1‑149 (MQ=255)
      ggACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCt                                                                                                                                                 >  6:291743/1‑149 (MQ=255)
        aCATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGc                                                                                                                                               <  1:269397/149‑1 (MQ=255)
        aCATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTAATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGc                                                                                                                                               >  5:53077/1‑149 (MQ=255)
          aTCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCAGGGGCTGcaa                                                                                                                                             <  2:301633/149‑3 (MQ=255)
          aTCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCAGGGGCTGcaa                                                                                                                                             <  5:56226/149‑3 (MQ=255)
           tCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCAc                                                                                                                                           >  7:137174/1‑149 (MQ=255)
             ccGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACg                                                                                                                                          >  4:120365/1‑148 (MQ=255)
               gCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCATGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACg                                                                                                                                          >  2:138740/1‑146 (MQ=255)
                cTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGa                                                                                                                                      <  8:126345/149‑1 (MQ=255)
                       gggCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGTAGCCCCTGGGCTGCCACGAAGACAGCCAg                                                                                                                               <  5:225092/149‑1 (MQ=255)
                        ggCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAg                                                                                                                               <  4:298708/148‑1 (MQ=255)
                        ggCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAg                                                                                                                               <  3:4483/148‑1 (MQ=255)
                          cAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGc                                                                                                                            <  2:153690/149‑1 (MQ=255)
                                  gcgTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGtt                                                                                                                    <  1:212646/149‑1 (MQ=255)
                                  gcgTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGtt                                                                                                                    >  5:217233/1‑149 (MQ=255)
                                    gTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTc                                                                                                                   >  7:94851/1‑148 (MQ=255)
                                       tCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAg                                                                                                                <  6:217233/148‑1 (MQ=255)
                                       tCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAg                                                                                                                <  5:291743/148‑1 (MQ=255)
                                           gATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAggg                                                                                                              <  4:20810/146‑1 (MQ=255)
                                                 aGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGc                                                                                                                                                  <  7:82687/105‑1 (MQ=255)
                                                 aGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGc                                                                                                                                                  >  8:82687/1‑105 (MQ=255)
                                                 aGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTcagc                                                                                                        >  3:287314/1‑146 (MQ=255)
                                                 aGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTcagc                                                                                                        <  4:287314/146‑1 (MQ=255)
                                                      cAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCgg                                                                                                <  8:137174/149‑1 (MQ=255)
                                                      cAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCgg                                                                                                <  4:60655/149‑1 (MQ=255)
                                                         ccTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGc                                                                                               >  2:108196/1‑147 (MQ=255)
                                                          cTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAg                                                                                            >  3:6799/1‑149 (MQ=255)
                                                            cATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATc                                                                                                                         <  4:178476/118‑1 (MQ=255)
                                                            cATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATc                                                                                                                         >  3:178476/1‑118 (MQ=255)
                                                             aTCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAgcg                                                                                          <  5:72613/148‑1 (MQ=255)
                                                              tCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAgcggc                                                                                        >  3:70612/1‑149 (MQ=255)
                                                                 aCTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGTCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGcgat                                                                                     <  8:94851/149‑1 (MQ=255)
                                                                       cccATCATGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGa                                                                               <  1:138740/149‑1 (MQ=255)
                                                                       cccATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGa                                                                               >  5:134076/1‑149 (MQ=255)
                                                                        ccATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGAc                                                                              >  8:40415/1‑149 (MQ=255)
                                                                          aTCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGa                                                                            >  1:321184/1‑149 (MQ=255)
                                                                               ggTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATAccc                                                                       >  6:80796/1‑149 (MQ=255)
                                                                                gTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAgg                                                                                                                              <  2:273024/93‑1 (MQ=255)
                                                                                gTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAgg                                                                                                                              >  1:273024/1‑93 (MQ=255)
                                                                                   cgcTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCgg                                                                   >  5:9230/1‑149 (MQ=255)
                                                                                      tGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGt                                                                >  6:244013/1‑149 (MQ=255)
                                                                                          ccTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCa                                                             <  4:292843/148‑1 (MQ=255)
                                                                                          ccTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCa                                                             <  6:9230/148‑1 (MQ=255)
                                                                                             gATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAccc                                                          >  8:156393/1‑148 (MQ=255)
                                                                                             gATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGcc                                                          >  1:175003/1‑148 (MQ=255)
                                                                                             gATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGcc                                                          >  3:62628/1‑148 (MQ=255)
                                                                                                 gAAAGCGGCGAATGCGTGGCTGGCCCCTGTGAGTGTGAGCAGTTGCAGCACCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGt                                                      <  6:53077/148‑1 (MQ=255)
                                                                                                  aaaGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGc                                                    <  2:321184/149‑1 (MQ=255)
                                                                                                   aaGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGcc                                                   <  6:103166/149‑1 (MQ=255)
                                                                                                       ggcgAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAg                                               <  3:188884/149‑1 (MQ=255)
                                                                                                                   gCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTc                                    >  1:52276/1‑148 (MQ=255)
                                                                                                                    cAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGc                                             <  4:70963/138‑1 (MQ=255)
                                                                                                                    cAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGc                                             >  3:70963/1‑138 (MQ=255)
                                                                                                                        ccccTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCAcccc                              <  6:134076/149‑1 (MQ=255)
                                                                                                                              tgAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGc                        >  6:84128/1‑149 (MQ=255)
                                                                                                                                 gggTGAGCAGTTGCAGCCCCTGGGCTGCCAGCGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTa                      >  1:194874/1‑149 (MQ=255)
                                                                                                                                 gggTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTAc                     <  5:244013/149‑1 (MQ=255)
                                                                                                                                 gggTGAGCAGTTGCAGCCCCTGGGCTACCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTAc                     >  1:322344/1‑149 (MQ=255)
                                                                                                                                    tGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgc                  >  2:67851/1‑149 (MQ=255)
                                                                                                                                       gCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGGGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGACAGCACCTCCACCCCGTCCGCTACGgcagc               >  5:187218/1‑149 (MQ=255)
                                                                                                                                          gTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcg              >  7:299167/1‑147 (MQ=255)
                                                                                                                                          gTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcg              <  8:299167/147‑1 (MQ=255)
                                                                                                                                           ttGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcgcg            >  2:267718/1‑148 (MQ=255)
                                                                                                                                           tggCAGCCCAGGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGc                                             >  8:298991/3‑115 (MQ=255)
                                                                                                                                           tgccAGCCCAGGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGc                                             <  7:298991/112‑1 (MQ=255)
                                                                                                                                             gCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcgcgca          <  4:231747/148‑1 (MQ=255)
                                                                                                                                             gCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcgcg            >  3:231747/1‑146 (MQ=255)
                                                                                                                                              cAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGgcagcgcgcagc        <  4:70612/149‑1 (MQ=255)
                                                                                                                                                 ccccTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGGCAGCGCGCAGCGcc     >  7:281216/1‑149 (MQ=255)
                                                                                                                                                    cTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGGCAGCGCGCAGCGCCTgg  <  7:40415/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGCGCCAGCCTCATCCACTGCGCCCATCACGGTGCGCTGTTCCTGATAGAAAGCGGCGAATGCGTGGCAGGCCCCTGTGAGGGTGAGCAGTTGCAGCCCCTGGGCTGCCACGAAGACAGCCAGGGCATCTGGTTCAAGGGGTCAGCAGCACCGGCAAGCGGCGATCGATGACGATACCCTCGGCGTCCAGCCGAGTGCCATAGGCCAGCACCTCCACCCCGTCCGCTACGGCAGCGCGCAGCGCCTGG  >  NC_002947/5329493‑5329789

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