Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,337,673 A→T 7.3% M105L (ATG→TTG)  PP_3812 → lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,337,6730AT 7.3% 185.2 / 3.0 64M105L (ATG→TTG) PP_3812lipase
Reads supporting (aligned to +/- strand):  ref base A (26/33);  new base T (2/3);  total (28/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

CGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGATCAAC  >  NC_002947/4337526‑4337819
                                                                                                                                                   |                                                                                                                                                  
cgcgcgGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCAt                                                                                                                                                   <  3:364126/149‑1 (MQ=255)
cgcgcgGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCAt                                                                                                                                                   <  2:40481/149‑1 (MQ=255)
       tACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                                                            <  8:306189/149‑1 (MQ=255)
            aCCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATCTTCTTCCAc                                                                                                                                         <  6:19786/147‑1 (MQ=255)
              cGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATg                                                                                                                                                  >  5:247105/1‑136 (MQ=255)
              cGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATg                                                                                                                                                  <  6:247105/136‑1 (MQ=255)
              cGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCAcggcg                                                                                                                                     >  8:259022/1‑149 (MQ=255)
                 tGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCAcggcggcg                                                                                                                                  <  7:86494/149‑1 (MQ=255)
                   cccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGt                                                                                                                                >  2:14053/1‑149 (MQ=255)
                      aGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCATTGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCATCCACGGCGGCGGTTgg                                                                                                                             <  1:23268/149‑1 (MQ=255)
                      aGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                                                            <  6:146453/134‑1 (MQ=255)
                      aGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                                                            >  5:146453/1‑134 (MQ=255)
                             ggTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTgg                                                                                                                      <  3:142359/149‑1 (MQ=255)
                               tAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTgg                                                                                                                             <  1:184290/140‑1 (MQ=255)
                               tAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTgg                                                                                                                             <  7:62214/140‑1 (MQ=255)
                               tAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTgg                                                                                                                             >  8:62214/1‑140 (MQ=255)
                               tAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTgg                                                                                                                             >  2:184290/1‑140 (MQ=255)
                                aaaGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGc                                                                                                                    >  3:127602/1‑148 (MQ=255)
                                    gTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCa                                                                                                                                          <  3:79329/122‑1 (MQ=255)
                                    gTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCa                                                                                                                                          >  4:79329/1‑122 (MQ=255)
                                    gTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCAcg                                                                                                                                        >  5:101934/1‑124 (MQ=255)
                                    gTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCAcg                                                                                                                                        <  6:101934/124‑1 (MQ=255)
                                       gACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTccc                                                                                                            <  1:14053/149‑1 (MQ=255)
                                         ccTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGa                                                                                                          >  1:284504/1‑149 (MQ=255)
                                            aTCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCAt                                                                                                                                                   >  8:204308/1‑105 (MQ=255)
                                            aTCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCAt                                                                                                                                                   <  7:204308/105‑1 (MQ=255)
                                            aTCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGAcc                                                                                                        >  2:122892/1‑148 (MQ=255)
                                                ggggTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATAAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGATGGGTGCTGGGCGACTTCCCGAcccca                                                                                                     >  2:259101/1‑145 (MQ=255)
                                                    tCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGa                                                                                               <  2:167830/149‑1 (MQ=255)
                                                           gAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGACGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATcc                                                                                        >  3:212292/1‑149 (MQ=255)
                                                                     aCGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCtggt                                                                              >  3:163220/1‑149 (MQ=255)
                                                                        aTTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCTTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTggtggt                                                                           >  4:168748/1‑149 (MQ=255)
                                                                        aTTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGGGATCTGCCGGTGTTCTTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTggtggt                                                                           >  4:163181/1‑149 (MQ=255)
                                                                         ttCAGGCCAATGGCCATTCGCTCAAGCTGCAAGGGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTggtggt                                                                           <  7:228230/148‑1 (MQ=255)
                                                                            aGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGc                                                                       >  3:225505/1‑149 (MQ=255)
                                                                             ggCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCtggtg                                                                             <  1:205194/142‑1 (MQ=255)
                                                                             ggCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCACGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGCTCCGCGACCtggtg                                                                             >  2:205194/1‑142 (MQ=255)
                                                                                cAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACATCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCgg                                                                   <  1:259101/149‑1 (MQ=255)
                                                                                  aTGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGt                                                                 <  1:152663/149‑1 (MQ=255)
                                                                                     gCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCTTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCg                                                              <  3:168748/149‑1 (MQ=255)
                                                                                     gCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCTTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCg                                                              <  3:163181/149‑1 (MQ=255)
                                                                                     gCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCg                                                              >  2:13929/1‑149 (MQ=255)
                                                                                         gTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTgg                                                          >  3:227193/1‑149 (MQ=255)
                                                                                             aTCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGt                                                      <  1:13929/149‑1 (MQ=255)
                                                                                             aTCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGt                                                      <  4:227193/149‑1 (MQ=255)
                                                                                                aaGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCt                                                    <  1:122892/148‑1 (MQ=255)
                                                                                                   cTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGTCGGCTGTTGGGTGCTGTGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTTGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgt                                                <  6:210223/149‑1 (MQ=255)
                                                                                                      cAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgtg                                               <  1:327001/147‑1 (MQ=255)
                                                                                                      cAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgt                                                <  4:212292/146‑1 (MQ=255)
                                                                                                          tggtgCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTac                                         <  4:163220/149‑1 (MQ=255)
                                                                                                               cgcCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACAccc                                     <  4:127602/148‑1 (MQ=255)
                                                                                                                gcCCGACTAATGTCAAAGGTCATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACGAGCGACTGATCCGCGACCTGGTGGTCGGCTAGGGTGCGGTGGCCGTCTATGTGGACTAAAcccc                                    >  7:292790/1‑148 (MQ=255)
                                                                                                                     aCTAATGTCAAAGGTGATCTGCCGGTGTTCTTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCAcc                              <  4:225505/149‑1 (MQ=255)
                                                                                                                       tAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCg                             >  4:102015/1‑148 (MQ=255)
                                                                                                                         aTGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGaa                          <  7:259022/149‑1 (MQ=255)
                                                                                                                         aTGTCAAAGGTCATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGaa                          <  2:284504/149‑1 (MQ=255)
                                                                                                                            tCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAgcg                       >  4:322767/1‑149 (MQ=255)
                                                                                                                              aaaGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCa                     >  6:331977/1‑149 (MQ=255)
                                                                                                                              aaaGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCa                     >  3:365642/1‑149 (MQ=255)
                                                                                                                                  gTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTAc                 <  3:322767/149‑1 (MQ=255)
                                                                                                                                            ggTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGa       <  3:102015/149‑1 (MQ=255)
                                                                                                                                             gtgtTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAGGCGCGCTACCCGACGGCgct      >  4:72206/1‑147 (MQ=255)
                                                                                                                                                  cATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGCCCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGATCAAc  >  5:304878/1‑148 (MQ=255)
                                                                                                                                                  cATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGATCAAc  <  6:304878/148‑1 (MQ=255)
                                                                                                                                                   aTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACgg          >  3:43525/1‑139 (MQ=255)
                                                                                                                                                   aTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACgg          <  4:43525/139‑1 (MQ=255)
                                                                                                                                                   aTGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGATCAAc  >  6:41728/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTACACCCCCTCACCGGAAGCGCACTACCCGACGGCGATCAAC  >  NC_002947/4337526‑4337819

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: ‑