Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,066,342 G→A 7.2% A86A (GCC→GCT)  PP_0924 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,066,3420GA 7.2% 201.8 / 4.3 69A86A (GCC→GCT) PP_0924hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (3/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGAA  >  NC_002947/1066194‑1066488
                                                                                                                                                    |                                                                                                                                                  
gTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTggcg                                                                                                                                                    <  5:103048/149‑1 (MQ=255)
      cccAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCATCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTgg                                                                                                                                              >  6:128201/1‑149 (MQ=255)
        cAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgggc                                                                                                                                            <  8:322234/149‑3 (MQ=255)
          aaGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGAt                                                                                                                                          <  1:219914/149‑1 (MQ=255)
          aaGCACGGCGATGGAACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGAt                                                                                                                                          <  8:321813/149‑1 (MQ=255)
           aGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                         <  3:392610/149‑1 (MQ=255)
           aGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                         <  2:7947/149‑1 (MQ=255)
            gCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATc                                                                                                                                         >  7:49730/1‑148 (MQ=255)
             cACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGa                                                                                                                                       <  8:294929/149‑1 (MQ=255)
                  cGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAgg                                                                                                                                  <  3:138665/149‑1 (MQ=255)
                         aCCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGggtgg                                                                                                                           >  7:1815/1‑149 (MQ=255)
                           cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         >  1:158458/1‑149 (MQ=255)
                           cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         <  6:177254/149‑1 (MQ=255)
                           cgccgTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGgtggtg                                                                                                                         >  8:184409/1‑149 (MQ=255)
                                        gTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGc                                                                                                            <  8:49730/149‑1 (MQ=255)
                                         tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                              >  1:355962/1‑114 (MQ=255)
                                         tCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCCGCCAGCAGCCTgg                                                                                                                                              <  2:355962/114‑1 (MQ=255)
                                           gggggCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGAAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCg                                                                                                          <  7:184409/148‑1 (MQ=255)
                                                    gCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGt                                                                                                >  2:8555/1‑149 (MQ=255)
                                                     cAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTc                                                                                               <  2:49584/149‑1 (MQ=255)
                                                     cAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTc                                                                                               <  4:9717/149‑1 (MQ=255)
                                                      aGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcg                                                                                              >  2:118197/1‑149 (MQ=255)
                                                         tACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTACAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTcgcg                                                                                            >  2:67572/1‑148 (MQ=255)
                                                           cGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTa                                                                                          <  3:282230/148‑1 (MQ=255)
                                                           cGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTAGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTa                                                                                          <  1:310655/148‑1 (MQ=255)
                                                             gCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCATGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTc                                                                                                              >  5:110111/1‑126 (MQ=255)
                                                             gCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCATGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTc                                                                                                              <  6:110111/126‑1 (MQ=255)
                                                              cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       >  5:152013/1‑148 (MQ=255)
                                                              cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       <  8:68427/148‑1 (MQ=255)
                                                              cccAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTacac                                                                                       >  7:68427/1‑148 (MQ=255)
                                                                      gAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggc                                                                              >  2:298524/1‑149 (MQ=255)
                                                                      gAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTggc                                                                              >  2:8186/1‑149 (MQ=255)
                                                                             gggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGt                                                                        <  3:365536/148‑1 (MQ=255)
                                                                              ggCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTg                                                                       >  2:113083/1‑148 (MQ=255)
                                                                                ccTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAc                                                                    >  3:165354/1‑149 (MQ=255)
                                                                                 cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcg                                                                   <  1:45022/149‑1 (MQ=255)
                                                                                 cTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCAcg                                                                   <  3:196632/149‑1 (MQ=255)
                                                                                  ttGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGc                                                                      >  3:281084/1‑145 (MQ=255)
                                                                                  ttGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGc                                                                      <  4:281084/145‑1 (MQ=255)
                                                                                      aGCTACAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTa                                                              <  1:67572/149‑1 (MQ=255)
                                                                                       gCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAg                                                             >  3:355911/1‑149 (MQ=255)
                                                                                          gCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcg                                                          >  4:45704/1‑149 (MQ=255)
                                                                                           cAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcg                                                          <  1:298524/148‑1 (MQ=255)
                                                                                            aaCAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAgtcgtc                                                        >  4:292168/1‑149 (MQ=255)
                                                                                                          gCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCt                                          >  1:194144/1‑149 (MQ=255)
                                                                                                           cTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTg                                         <  8:1815/149‑1 (MQ=255)
                                                                                                              aaTTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTc                                      <  3:341549/149‑1 (MQ=255)
                                                                                                                ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCgg                                    >  3:313837/1‑149 (MQ=255)
                                                                                                                ttcttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCgg                                    >  7:187229/1‑149 (MQ=255)
                                                                                                                   ttcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAc                                 >  6:52245/1‑149 (MQ=255)
                                                                                                                    tcttcGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCAcc                                >  4:71535/1‑149 (MQ=255)
                                                                                                                         gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTg                            <  3:45704/148‑1 (MQ=255)
                                                                                                                         gCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTg                            <  7:245982/148‑1 (MQ=255)
                                                                                                                          cTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTTGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAg                          <  5:128201/149‑1 (MQ=255)
                                                                                                                             tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTGGTAGTCGTCCATGTCCAGCCCCTGCTCGGCACCtttta                           >  6:265007/1‑143 (MQ=255)
                                                                                                                             tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCtt                       >  7:37856/1‑149 (MQ=255)
                                                                                                                             tAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTGTGAGCtt                       >  4:231607/1‑149 (MQ=255)
                                                                                                                             tAGAAACCGCCCAGCCTGCTGGCGGCCTGGTGATCGAGCAGGCGGCTGGTGGCAACGAGTTCGCCGGTGCTCAGGACGCGTACACCAATGTGCCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGa                           <  5:265007/145‑1 (MQ=255)
                                                                                                                               gAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTg                      <  3:400981/148‑1 (MQ=255)
                                                                                                                                tcgatcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc                    >  2:203550/8‑149 (MQ=255)
                                                                                                                                aaaCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGc                    >  8:212443/1‑149 (MQ=255)
                                                                                                                                  aCGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCt                  <  2:194144/149‑1 (MQ=255)
                                                                                                                                   cGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 <  1:8186/149‑1 (MQ=255)
                                                                                                                                   atcaccAGGCCGCCAGCAGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCtt                 >  7:55355/5‑149 (MQ=255)
                                                                                                                                     gCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGa               <  4:165354/149‑1 (MQ=255)
                                                                                                                                      cccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGaa              <  6:152013/149‑1 (MQ=255)
                                                                                                                                       ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt             >  6:20476/1‑149 (MQ=255)
                                                                                                                                       ccAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGCTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAAt             >  1:266481/1‑149 (MQ=255)
                                                                                                                                          ggCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCt          >  5:1112/1‑149 (MQ=255)
                                                                                                                                             tgctgGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTgcg       >  3:328541/1‑149 (MQ=255)
                                                                                                                                                  gcggcCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGaa  >  2:158138/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTCGCCCCAAAGCACGGCGATGGTACCGCCGTTGCGGTAGTCGGGGGCGACGCAGGTACGGCCCAGTTCGAGGATCGGGCCTTGCAGCTGCAACAGGCCGTGCAGGCTGAATTCTTCTTCGCTGTAGAAACGGCCCAGGCTGCTGGCGGCCTGGTGATCGAGCAGGCGGGTGGTGGCAACCAGTTCGCCGGTGCTCAGGTCGCGTACACCAATGTGGCGGCAGTGCACGTCGTAGTCGTCCATGTCCAGGCCCTGCTCGGCACCTTTGAGCTTGGCCTTGAATTCTGCGCTGAA  >  NC_002947/1066194‑1066488

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