Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,081,045 T→A 12.5% intergenic (+36/+26) PP_5332 → / ← PP_5333 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,081,0450TA12.5% 189.2 / 18.7 72intergenic (+36/+26)PP_5332/PP_5333hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (5/4);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GGATGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGAT  >  NC_002947/6080896‑6081193
                                                                                                                                                     |                                                                                                                                                    
ggatggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCtt                                                                                                                                                      <  1:1038871/150‑1 (MQ=255)
 gatggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTg                                                                                                                                                     >  1:347410/1‑150 (MQ=255)
         cGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             <  1:1193102/150‑1 (MQ=255)
           aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           <  2:1936902/150‑1 (MQ=255)
           aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           <  1:2148412/150‑1 (MQ=255)
                tgtgCATCGAGACAGCGAATGTTTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  2:2379211/150‑1 (MQ=255)
                tgtgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      >  2:1824782/1‑150 (MQ=255)
                  tgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCggg                                                                                                                                    <  2:2041031/150‑1 (MQ=255)
                   gCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                              >  2:621967/1‑139 (MQ=255)
                   gCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                              <  1:621967/139‑1 (MQ=255)
                             aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCaa                                                                                                                         >  1:2158422/1‑150 (MQ=255)
                                gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAgg                                                                                                                      >  1:1478061/1‑150 (MQ=255)
                                 aaTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGGGTAGAGGCCGGGGCAATTATCAAAGGt                                                                                                                     >  1:1074/1‑150 (MQ=255)
                                     tgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                 <  2:1088350/150‑1 (MQ=255)
                                     tgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                 <  1:2138091/150‑1 (MQ=255)
                                       tgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGt                                                                                                               <  1:2109486/150‑1 (MQ=255)
                                                     ggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTATCGGACGCTTCGCGGGCGTGCCCGCGATGAGGCAGGAGATATTATCAACGCACTGCGACTTCCTCCACCCTc                                                                                                 >  1:1223633/1‑150 (MQ=255)
                                                         ttGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCt                                                                                             >  1:453527/1‑150 (MQ=255)
                                                             aGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTc                                                                                         <  2:347410/150‑1 (MQ=255)
                                                               ccTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTg                                                                                       >  1:2346254/1‑150 (MQ=255)
                                                                cTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGc                                                                                      <  1:1462800/150‑1 (MQ=255)
                                                                 tGGGGCGAGTCATTCGCGGGGCGTCCAGATCGGTTGTGGAAGTATCTGAGCCCTGCTCTTGTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCa                                                                                     <  2:2361899/150‑1 (MQ=255)
                                                                            aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                                    >  2:566835/1‑124 (MQ=255)
                                                                            aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                                    <  1:566835/124‑1 (MQ=255)
                                                                             ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAg                                                                                <  1:800475/143‑1 (MQ=255)
                                                                             ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAg                                                                                >  2:800475/1‑143 (MQ=255)
                                                                                gCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGcc                                                                      >  1:2444328/1‑150 (MQ=255)
                                                                                        tCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGc                                                              <  2:2346254/150‑1 (MQ=255)
                                                                                        tCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGc                                                              >  2:701580/1‑150 (MQ=255)
                                                                                          gaGATCGGTTGTGAAATTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            <  2:1223633/149‑1 (MQ=255)
                                                                                          cAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            <  2:1630731/150‑1 (MQ=255)
                                                                                          cAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            >  1:1665514/1‑150 (MQ=255)
                                                                                           aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:253197/71‑1 (MQ=255)
                                                                                           aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                          >  2:253197/1‑71 (MQ=255)
                                                                                           aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAg                                                           <  2:317226/150‑1 (MQ=255)
                                                                                           aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAg                                                           <  1:1274840/150‑1 (MQ=255)
                                                                                                     tgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCCGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTg                                                 >  1:2200904/1‑150 (MQ=255)
                                                                                                       aaaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGTGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATTCAGCGCATAGGCCGCATTGGCCAGGTCATTGCTgcc                                               >  2:1004499/1‑148 (MQ=255)
                                                                                                       aaaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               <  2:1478061/150‑1 (MQ=255)
                                                                                                        aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATGGGCCGCATCGGCCAGGTCATTGCTGAcc                                              <  2:2152032/150‑1 (MQ=255)
                                                                                                        aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAcc                                              <  2:1665514/150‑1 (MQ=255)
                                                                                                               tGAGCCCTGTTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCa                                       <  2:1074/150‑1 (MQ=255)
                                                                                                               tGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCa                                       >  2:1476712/1‑150 (MQ=255)
                                                                                                                     cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                                 <  1:1824782/150‑1 (MQ=255)
                                                                                                                     cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                                 >  2:174862/1‑150 (MQ=255)
                                                                                                                          ctTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCATCTGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTg                            <  1:1208231/150‑1 (MQ=255)
                                                                                                                           tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGt                                                                                          <  2:264874/87‑1 (MQ=255)
                                                                                                                           tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGt                                                                                          >  1:264874/1‑87 (MQ=255)
                                                                                                                           tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                           <  1:1476712/150‑1 (MQ=255)
                                                                                                                           tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                           <  1:174862/150‑1 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                           >  2:427445/1‑121 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                           <  1:427445/121‑1 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCTTAAACATAGGCCGCATCGGCCAGGGCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:643834/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  2:1441895/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:1170359/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  2:870636/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCGGAGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:1217445/1‑150 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               >  1:686128/1‑115 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               <  2:686128/115‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgg           >  1:1889107/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAcc                                                                                              <  2:915803/67‑1 (MQ=37)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAcc                                                                                              >  1:915803/1‑67 (MQ=37)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCa                                                       >  1:613374/1‑106 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCa                                                       <  2:613374/106‑1 (MQ=255)
                                                                                                                                            tCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg          <  1:870636/150‑1 (MQ=255)
                                                                                                                                            tCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg          >  2:2188302/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgggg         <  1:1441895/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgt        <  1:701580/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgt        <  2:453527/150‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               <  1:1101159/94‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               <  1:1102499/94‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               >  2:1101159/1‑94 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               >  2:1102499/1‑94 (MQ=255)
                                                                                                                                                   tttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGa   <  2:1856110/149‑1 (MQ=255)
                                                                                                                                                    ttGNCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  >  2:1506930/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  >  2:1352993/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  <  2:1889107/150‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGATTGACGCTGACCCACAGCGGGGCGTAGat  >  1:1697577/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGATGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGAT  >  NC_002947/6080896‑6081193

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