Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,744,888 T→A 8.0% intergenic (‑14/+327) PP_1542 ← / ← PP_1543 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,744,8880TA 8.0% 218.8 / 6.8 74intergenic (‑14/+327)PP_1542/PP_1543hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base A (3/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GGCGTCGTCATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAGGATCA  >  NC_002947/1744746‑1745028
                                                                                                                                              |                                                                                                                                            
ggCGTCGTCATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCg                                                                                                                                         <  4:324644/148‑1 (MQ=255)
 gCGTCGTCATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGgt                                                                                                                                       >  5:158898/1‑149 (MQ=255)
    tcgtcATAGCGTGGTCTGGTGAGGAATTGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCTGTCCCAAGGTGTTCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACGGGCCAAGCCGgtgct                                                                                                                                    <  2:141025/149‑3 (MQ=255)
      gtcATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAATGTGTTCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGtgttgt                                                                                                                                  <  6:296846/149‑1 (MQ=255)
      gtcATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGtgttgt                                                                                                                                  >  2:259939/1‑149 (MQ=255)
            gcctggcctgGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGa                                                                                                                            <  3:232052/149‑1 (MQ=255)
             ccTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGACTTTCTCCATGCAAGCACCGGCTTgg                                                                                                                                            <  8:348268/132‑1 (MQ=255)
             ccTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGACTTTCTCCATGCAAGCACCGGCTTgg                                                                                                                                            >  7:348268/1‑132 (MQ=255)
                    tGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTg                                                                                                                                             <  1:122063/124‑1 (MQ=255)
                    tGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTg                                                                                                                                             >  2:122063/1‑124 (MQ=255)
                     ggTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGTCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTtata                                                                                                                          <  5:293199/142‑1 (MQ=255)
                     ggTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCACAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTATCCATGCAAGCACCGGCTTGGGCGGTGTTGTGCCTGATa                                                                                                                          >  6:293199/1‑142 (MQ=255)
                      gTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAAt                                                                                                                  >  8:117601/1‑149 (MQ=255)
                        gAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAAtt                                                                                                                 >  5:298553/1‑148 (MQ=255)
                          ggAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTaa                                                                                                              <  4:226347/149‑1 (MQ=255)
                          ggAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCCAAGCCGgtgct                                                                                                                                    <  6:188828/127‑3 (MQ=255)
                          ggAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCCAAGCCGgtgct                                                                                                                                    >  5:188828/1‑125 (MQ=255)
                           gAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTaa                                                                                                              <  8:334955/148‑1 (MQ=255)
                           gAAATCGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTaa                                                                                                              >  7:334955/1‑148 (MQ=255)
                                 gTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCAt                                                                                                        <  4:219108/148‑1 (MQ=255)
                                      tGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAGGCACCGGCTTGGCCGGtgttgt                                                                                                                                  >  2:48341/1‑117 (MQ=255)
                                      tGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCACGCACCGGCTTGGCCGGtgttgt                                                                                                                                  <  1:48341/117‑1 (MQ=255)
                                      tGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGt                                                                                                  >  5:63939/1‑149 (MQ=255)
                                        ttGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCTAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGGATTTGGGTTa                                                                                                >  2:249409/1‑149 (MQ=255)
                                          gCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTatt                                                                                              >  4:251685/1‑149 (MQ=255)
                                          gCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTat                                                                                               <  7:117601/148‑1 (MQ=255)
                                            cGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTattga                                                                                            <  1:17132/149‑1 (MQ=255)
                                              gTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTg                                                                                                                                             >  8:273943/1‑98 (MQ=255)
                                              gTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTg                                                                                                                                             <  7:273943/98‑1 (MQ=255)
                                              gTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTattgatt                                                                                          <  2:88393/149‑1 (MQ=255)
                                              gTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTattgatt                                                                                          >  2:159823/1‑149 (MQ=255)
                                                 cGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGt                                                                                        <  6:158898/148‑1 (MQ=255)
                                                 cGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGt                                                                                        >  3:192970/1‑148 (MQ=255)
                                                     aTCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATAc                                                                                   >  1:215326/1‑149 (MQ=255)
                                                         cGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTg                                                                                                     >  1:113417/1‑127 (MQ=255)
                                                         cGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTg                                                                                                     <  2:113417/127‑1 (MQ=255)
                                                         cGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGt                                                                               <  4:192970/149‑1 (MQ=255)
                                                                 gTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACgaga                                                                       <  2:215326/149‑1 (MQ=255)
                                                                     cAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTggg                                                                                                   <  7:275238/117‑1 (MQ=255)
                                                                     cAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTggg                                                                                                   >  8:275238/1‑117 (MQ=255)
                                                                         gtgtGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATTCCAGTCAACGAGATGAATAcc                                                               >  5:274442/1‑149 (MQ=255)
                                                                          tgtgCCAGAGCAAGTAGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATAccc                                                              <  6:63939/149‑1 (MQ=255)
                                                                            tgCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGa                                                            >  3:250260/1‑149 (MQ=255)
                                                                              ccAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATg                                                          <  1:159823/149‑1 (MQ=255)
                                                                              ccAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATg                                                          >  6:104287/1‑149 (MQ=255)
                                                                              ccAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGACAACGAGATGAATACCCGATg                                                          <  3:251685/149‑1 (MQ=255)
                                                                               cAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATAc                                                                                   >  2:90908/1‑123 (MQ=255)
                                                                                            gctgctTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGAt                                                                                                                           >  4:96520/1‑70 (MQ=255)
                                                                                                   ctcGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCa                                      <  4:250260/148‑1 (MQ=255)
                                                                                                      ggTGGCTTCGATCTTGGCTTTCTCCATGCAATCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGTCAATTTAATGGCATTTGTGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAg                                     <  1:249409/146‑1 (MQ=255)
                                                                                                        tGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAAt                                >  6:328194/1‑149 (MQ=255)
                                                                                                             tctaTCATGGCTTTCTCCATGCAAGCACCGGCTTGGACGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCgca                            <  5:277965/145‑1 (MQ=255)
                                                                                                              cGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCgcag                           <  1:354684/148‑1 (MQ=255)
                                                                                                                aTCATGGCTTTCTCCATGCAAGCACCGGCGTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGcagcag                        >  2:1757/1‑149 (MQ=255)
                                                                                                                 tCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCgca                            >  6:277965/1‑144 (MQ=255)
                                                                                                                      gCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCtctc                   <  5:104287/148‑1 (MQ=255)
                                                                                                                              cATGCAAGCACCGGCCTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATtg            <  1:1757/147‑1 (MQ=255)
                                                                                                                               aTGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCgcag                           >  8:3243/1‑131 (MQ=255)
                                                                                                                               aTGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATtga           >  7:203689/1‑147 (MQ=255)
                                                                                                                               aTGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATtga           <  8:203689/147‑1 (MQ=255)
                                                                                                                                tGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATtgatga        <  5:328194/149‑1 (MQ=255)
                                                                                                                                 gCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAg       >  4:239350/1‑149 (MQ=255)
                                                                                                                                  cAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAgg      <  7:49205/149‑1 (MQ=255)
                                                                                                                                  acaacACCGGCCAAGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAg       >  4:263667/5‑148 (MQ=255)
                                                                                                                                  acaacACCGGCCAAGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAg       >  8:288181/5‑148 (MQ=255)
                                                                                                                                   aaGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACcgcg                                                   <  3:325295/103‑1 (MQ=255)
                                                                                                                                   aaGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAGGa     >  1:346583/1‑149 (MQ=255)
                                                                                                                                   aaGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAAACAGGAAATCGCAGAAGCTCTCGGCAATGATGAGGa     >  8:105205/1‑149 (MQ=255)
                                                                                                                                   aaGCACCGGCTTCGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACcgcg                                                   >  4:325295/1‑103 (MQ=255)
                                                                                                                                      cACCGGCCAAGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAGGATCa  <  7:288181/149‑1 (MQ=255)
                                                                                                                                         cGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCa                                      <  7:92942/110‑1 (MQ=255)
                                                                                                                                         cGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCa                                      >  8:92942/1‑110 (MQ=255)
                                                                                                                                           gCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAgg                                    >  8:151452/1‑110 (MQ=255)
                                                                                                                                           gCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAgg                                    <  7:151452/110‑1 (MQ=255)
                                                                                                                                             ttGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCAttt                                                                                                      >  5:6641/1‑42 (MQ=255)
                                                                                                                                             ttGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCAttt                                                                                                      <  6:6641/42‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GGCGTCGTCATAGCCTGGGCTGGTGAGGAAATGGTTTCTGTTGCCGGTACGGCATCCCGGGTTCAGTCCCAAGGTGTGCCAGAGCAAGTCGAGCTGCTTCTCGGTGGCTTCGATCATGGCTTTCTCCATGCAAGCACCGGCTTGGCCGGTGTTGTGCCTGATAGTGGCAATTTAATGGCATTTGGGTTATTGATTGTAATACCAGTCAACGAGATGAATACCCGATGGACCGCGAGCCAACAAACCAGGAAATCGCAGCAGCTCTCGGCATTGATGAGGATCA  >  NC_002947/1744746‑1745028

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: