Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,152 G→A 7.1% intergenic (‑52/+25) PP_3663 ← / ← pssA GGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,162,1520GA 7.1% 200.7 / 4.6 70intergenic (‑52/+25)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCGGCG  >  NC_002947/4162003‑4162298
                                                                                                                                                     |                                                                                                                                                  
gACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                                    >  1:1324099/1‑150 (MQ=255)
      aTGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                              <  2:1783006/150‑1 (MQ=255)
         cGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAggg                                                                                                                                           <  2:469705/150‑1 (MQ=255)
          gATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                          <  2:1252286/150‑1 (MQ=255)
          gATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                          <  2:950936/150‑1 (MQ=255)
               tGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTa                                                                                                                                     >  1:242131/1‑150 (MQ=255)
                  tCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAg                                                                                                                                  >  2:1910187/1‑150 (MQ=255)
                      gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGa                                                                                                                              >  2:1497834/1‑150 (MQ=255)
                          agagGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTc                                                                                                                          <  2:1324099/150‑1 (MQ=255)
                           gagGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCa                                                                                                                         <  2:1737398/150‑1 (MQ=255)
                             ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                                       >  1:339011/1‑150 (MQ=255)
                             ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                                       <  1:808464/150‑1 (MQ=255)
                             ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                                       <  2:725579/150‑1 (MQ=255)
                              gggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCagag                                                                                                                      >  1:485168/1‑150 (MQ=255)
                              gggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCagag                                                                                                                      >  2:1235131/1‑150 (MQ=255)
                                     gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAgg                                                                                                                                 >  2:1294493/1‑132 (MQ=255)
                                     gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAgg                                                                                                                                 <  1:1294493/132‑1 (MQ=255)
                                       aCTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCtt                                                                                                             >  1:362866/1‑150 (MQ=255)
                                         tGTTGCGTACACATTGTAATCCAGCGACTGACAAAGGGCCGTGTGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACGACTCTGGCCCTGTCGCTGGCTTGCCGGCTACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGa                                                                                                           <  2:68207/150‑1 (MQ=255)
                                         tGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGa                                                                                                           <  1:1235131/150‑1 (MQ=255)
                                          gTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACCGGCGATGTCAGAGCACCCGCTTGAt                                                                                                          <  2:896102/150‑1 (MQ=255)
                                                   acaTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTc                                                                                                 >  2:1743828/1‑150 (MQ=255)
                                                            tCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCt                                                                                        >  1:408808/1‑150 (MQ=255)
                                                            tCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCt                                                                                        >  2:431813/1‑150 (MQ=255)
                                                             ccAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCta                                                                                       >  1:1686117/1‑150 (MQ=255)
                                                              cAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCa                                                                                                                                        >  1:704901/1‑100 (MQ=255)
                                                              cAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCa                                                                                                                                        <  2:704901/100‑1 (MQ=255)
                                                              cAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCtat                                                                                      >  1:577048/1‑150 (MQ=255)
                                                               aGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCtata                                                                                     <  2:1686117/150‑1 (MQ=255)
                                                                 ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACa                                                                                                                                              <  1:407076/91‑1 (MQ=255)
                                                                 ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACa                                                                                                                                              >  2:407076/1‑91 (MQ=255)
                                                                 ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgc                                                                                   <  2:339011/150‑1 (MQ=255)
                                                                 ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgc                                                                                   <  1:1420632/150‑1 (MQ=255)
                                                                   ccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgcgc                                                                                 >  1:1601578/1‑150 (MQ=255)
                                                                    cTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCa                                                                                >  2:1731879/1‑150 (MQ=255)
                                                                          aaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGAc                                                                          <  2:242131/150‑1 (MQ=255)
                                                                            gggCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                                        >  1:357653/1‑150 (MQ=255)
                                                                            gggCCGTGGGCAGCAATATCATAATGACAGCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                                        >  1:1423449/1‑150 (MQ=255)
                                                                               ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                                    <  1:635873/71‑1 (MQ=255)
                                                                               ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                                    >  2:635873/1‑71 (MQ=255)
                                                                                 gTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGaag                                                                   <  2:485168/150‑2 (MQ=255)
                                                                                       agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                             <  2:88475/150‑1 (MQ=255)
                                                                                       agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                             >  1:88475/1‑150 (MQ=255)
                                                                                         caATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAAGCGGCTATAGCGCAACCGACGCAGAAACTTGCGCa                                                           <  1:1910187/150‑1 (MQ=255)
                                                                                              tCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgc                                                      >  1:169615/1‑150 (MQ=255)
                                                                                                aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                      >  1:1746110/1‑100 (MQ=255)
                                                                                                aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                      <  2:1746110/100‑1 (MQ=255)
                                                                                                 tAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCtata                                                                                     <  2:518682/116‑1 (MQ=255)
                                                                                                 tAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCtata                                                                                     >  1:518682/1‑116 (MQ=255)
                                                                                                    tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                                >  1:1757733/1‑150 (MQ=255)
                                                                                                        aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGtt                                            >  2:81621/1‑150 (MQ=255)
                                                                                                            aTGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGc                                        <  1:388920/150‑1 (MQ=255)
                                                                                                            aTGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGc                                        <  2:357653/150‑1 (MQ=255)
                                                                                                             tGCGCGGTGCTCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGcc                                       <  1:81621/150‑1 (MQ=255)
                                                                                                             tGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGcc                                       >  2:33458/1‑150 (MQ=255)
                                                                                                                gcgGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTtc                                    <  1:431813/150‑1 (MQ=255)
                                                                                                                  ggTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgctt                                                    >  1:1472071/1‑132 (MQ=255)
                                                                                                                  ggTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgctt                                                    <  2:1472071/132‑1 (MQ=255)
                                                                                                                       gcACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCa                             <  2:577048/150‑1 (MQ=255)
                                                                                                                                tgacccTGTCGCCGGCAAGCCAGCGACAGGGCCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCctct                     <  1:595974/146‑1 (MQ=39)
                                                                                                                                tgacccTGTCGCCGGCAAGCCAGCGACAGGGCCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCctct                     >  2:595974/4‑149 (MQ=39)
                                                                                                                                  gCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGcc                  >  1:1352651/1‑150 (MQ=255)
                                                                                                                                    ccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGa                >  2:551252/1‑150 (MQ=255)
                                                                                                                                    ccTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGa                >  1:1013574/1‑150 (MQ=255)
                                                                                                                                      tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGa                                                                                                           <  2:1476802/57‑1 (MQ=255)
                                                                                                                                      tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGa                                                                                                           >  1:1476802/1‑57 (MQ=255)
                                                                                                                                      tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAGCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAAc              <  2:1423449/150‑1 (MQ=255)
                                                                                                                                          gCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGAt          >  2:1376923/1‑150 (MQ=255)
                                                                                                                                           cTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGAtt         >  1:607398/1‑150 (MQ=255)
                                                                                                                                            tGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTc                                           <  2:215411/115‑1 (MQ=255)
                                                                                                                                            tGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTc                                           >  1:215411/1‑115 (MQ=255)
                                                                                                                                                ttGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCgg    >  2:1705800/1‑150 (MQ=255)
                                                                                                                                                  gCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCGgcg  <  1:70136/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCGGCG  >  NC_002947/4162003‑4162298

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: