Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,139 C→T 5.2% S92S (AGC→AGT hemB → 5‑aminolevulinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,313,1390CT 5.2% 230.6 / 2.4 77S92S (AGC→AGThemB5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (33/40);  new base T (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ATCGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGACCATG  >  NC_002947/3312990‑3313278
                                                                                                                                                     |                                                                                                                                           
aTCGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGc                                                                                                                                             >  1:977503/1‑150 (MQ=255)
aTCGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGc                                                                                                                                             <  1:884280/150‑1 (MQ=255)
aTCGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGc                                                                                                                                             <  1:783164/150‑1 (MQ=255)
  cGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGa                                                                                                                                           <  1:665004/150‑1 (MQ=255)
    atgatTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGaca                                                                                                                                         >  2:436416/1‑150 (MQ=255)
     tgatTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGacac                                                                                                                                        <  2:977503/150‑1 (MQ=255)
                     aCCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCCGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGc                                                                                                                        >  1:2234721/1‑150 (MQ=255)
                          tATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGt                                                                                                                   <  2:2064451/150‑1 (MQ=255)
                           aTGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGtt                                                                                                                  <  1:1124211/150‑1 (MQ=255)
                            tGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTg                                                                                                                 >  2:1872678/1‑150 (MQ=255)
                              cccGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTgcg                                                                                                               >  2:1402692/1‑150 (MQ=255)
                               ccGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTgcgc                                                                                                              <  1:1004868/150‑1 (MQ=255)
                               ccGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTgcgc                                                                                                              >  2:1964667/1‑150 (MQ=255)
                                     tCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTgc                                                                                                                >  1:1040718/1‑142 (MQ=255)
                                     tCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTgc                                                                                                                <  2:1040718/142‑1 (MQ=255)
                                      cAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGGCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTc                                                                                                       <  2:1143087/150‑1 (MQ=255)
                                       aaCCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTcg                                                                                                      >  1:2079516/1‑150 (MQ=255)
                                         ccGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGAGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTcgcg                                                                                                    >  2:433859/1‑150 (MQ=255)
                                          cGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGTGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGt                                                                                                   >  2:2181579/1‑150 (MQ=255)
                                          cGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGt                                                                                                   >  1:2254913/1‑150 (MQ=255)
                                              tCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCt                                                                                               <  1:1656084/150‑1 (MQ=255)
                                                cccGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGc                                                                                             >  2:1438984/1‑150 (MQ=255)
                                                      aaGCTGCTGGCCCAGGAAAGCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGac                                                                                       <  1:792111/150‑1 (MQ=255)
                                                        gctgctGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGacac                                                                                     <  2:2079516/150‑1 (MQ=255)
                                                        gctgctGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGacac                                                                                     >  2:895660/1‑150 (MQ=255)
                                                         ctgctgGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACAcc                                                                                    <  1:239335/150‑1 (MQ=255)
                                                               gCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCg                                                                              <  1:341426/150‑1 (MQ=255)
                                                                 ccAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTg                                                                                              >  1:1161832/1‑132 (MQ=255)
                                                                 ccAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGAGTATCTg                                                                                              <  2:1161832/132‑1 (MQ=255)
                                                                  cAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGaa                                                                           >  1:29082/1‑150 (MQ=255)
                                                                  cAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGaa                                                                           >  1:350757/1‑150 (MQ=255)
                                                                   aGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                <  1:1326750/144‑1 (MQ=255)
                                                                   aGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                >  2:1326750/1‑144 (MQ=255)
                                                                   aGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                <  1:157227/144‑1 (MQ=255)
                                                                   aGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                >  2:157227/1‑144 (MQ=255)
                                                                        aTCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtggtg                                                                     <  2:2234721/150‑1 (MQ=255)
                                                                        aTCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtggtg                                                                     <  2:2228004/150‑1 (MQ=255)
                                                                         tCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtggtgg                                                                    <  2:206795/150‑1 (MQ=255)
                                                                         tCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtggtgg                                                                    <  1:1402692/150‑1 (MQ=255)
                                                                          cGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATggtggt                                                                   <  1:818350/150‑1 (MQ=255)
                                                                          cGAACGGTATGCCCGCGCCGGGATCAAGTCGGTGGTGAGCTTTGGCGTATGGCACCACCTGGATGCCGATGGCAGCGACACCGGGAACGAAAACGGCCGGGTTGCGCGCATGGCGCGTATCGGCAAGGACACCGTGCCGGAAATggtggt                                                                   <  2:829491/150‑1 (MQ=255)
                                                                            aGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGa                                                                 >  2:786719/1‑150 (MQ=255)
                                                                                cTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTc                                                             >  1:277280/1‑150 (MQ=255)
                                                                                 tATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACCCCGTGCCGGAAATGGTGGTGATGTcc                                                            <  1:1438984/150‑1 (MQ=255)
                                                                                 tATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTcc                                                            <  1:1623792/150‑1 (MQ=255)
                                                                                 tATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTcc                                                            <  1:1964667/150‑1 (MQ=255)
                                                                                   ggCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGa                                                          <  2:632738/149‑1 (MQ=255)
                                                                                        gcgcCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCt                                                     >  1:951884/1‑150 (MQ=255)
                                                                                          gcCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGc                                                   <  2:951884/150‑1 (MQ=255)
                                                                                              gTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCt                                               >  1:2127855/1‑150 (MQ=255)
                                                                                                aTCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGc                                             <  2:277280/150‑1 (MQ=255)
                                                                                                aTCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGc                                             <  2:155252/150‑1 (MQ=255)
                                                                                                       cGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATaca                                      <  1:436416/150‑1 (MQ=255)
                                                                                                            atgaCCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTg                                                                                                                                     >  1:1177989/1‑50 (MQ=255)
                                                                                                            atgaCCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTg                                                                                                                                     <  2:1177989/50‑1 (MQ=255)
                                                                                                                ccTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGc                                             >  1:314454/1‑134 (MQ=255)
                                                                                                                ccTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGc                                             <  2:314454/134‑1 (MQ=255)
                                                                                                                 cTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtg                                                                        <  1:1384289/106‑1 (MQ=255)
                                                                                                                 cTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAAtg                                                                        >  2:1384289/1‑106 (MQ=255)
                                                                                                                  tttGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAACCGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAACTGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGc                           >  1:166292/1‑150 (MQ=255)
                                                                                                                  tttGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGc                           <  2:2127855/150‑1 (MQ=255)
                                                                                                                  tttGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGc                           >  1:549563/1‑150 (MQ=255)
                                                                                                                      gCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAg                                  <  2:812498/139‑1 (MQ=255)
                                                                                                                      gCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAg                                  >  1:812498/1‑139 (MQ=255)
                                                                                                                       cGTATCGCACCACCTGGATGCCGATGGCAGTGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                >  1:1923162/1‑92 (MQ=255)
                                                                                                                       cGTATCGCACCACCTGGATGCCGATGGCAGTGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGc                                                                                <  2:1923162/92‑1 (MQ=255)
                                                                                                                        gTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTgc                     <  1:433859/150‑1 (MQ=255)
                                                                                                                         tATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACAcc                                    <  2:907404/134‑1 (MQ=255)
                                                                                                                         tATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACAcc                                    >  1:907404/1‑134 (MQ=255)
                                                                                                                          aTCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTgcgg                   <  1:596582/150‑1 (MQ=255)
                                                                                                                                  accTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGa                                                          >  2:1774938/1‑103 (MQ=255)
                                                                                                                                  accTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGa                                                          <  1:1774938/103‑1 (MQ=255)
                                                                                                                                    cTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATaca                                      >  1:2132595/1‑121 (MQ=255)
                                                                                                                                    cTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATaca                                      <  2:2132595/121‑1 (MQ=255)
                                                                                                                                    cTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCAt         >  1:55213/1‑150 (MQ=255)
                                                                                                                                      ggATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGa       >  1:597098/1‑150 (MQ=255)
                                                                                                                                       gATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGAc      >  2:1041661/1‑150 (MQ=255)
                                                                                                                                         tGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGACCa    >  1:1392965/1‑150 (MQ=255)
                                                                                                                                           ccGATGGCAGTGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGACCATg  <  1:2181579/150‑1 (MQ=255)
                                                                                                                                           ccGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGACCATg  <  2:2254913/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
ATCGATGATTTCGTGCCGATCACCAGTATGCCCGGTGTCAACCGTATCCCCGAGAAGCTGCTGGCCCAGGAAATCGAGCGCTATGCCCGCGCCGGTATCAAGTCGGTGATGACCTTTGGCGTATCGCACCACCTGGATGCCGATGGCAGCGACACCTGGAACGAAAACGGCCTGGTTGCGCGCATGTCGCGTATCTGCAAGGACACCGTGCCGGAAATGGTGGTGATGTCCGACACCTGCTTCTGCGAATACACCAGCCACGGCCACTGCGGCGTGCTGCATGACCATG  >  NC_002947/3312990‑3313278

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