Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,629 G→A 6.5% D181D (GAC→GAT)  PP_2941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,344,6290GA 6.5% 223.9 / 3.3 78D181D (GAC→GAT) PP_2941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base A (3/2);  total (41/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 = 9.93e-01

CCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATCGGCGTTATGCCTCG  >  NC_002947/3344481‑3344769
                                                                                                                                                    |                                                                                                                                            
ccAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGt                                                                                                                                             >  2:224016/1‑150 (MQ=255)
ccAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGt                                                                                                                                             >  2:567637/1‑150 (MQ=255)
 cAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTc                                                                                                                                            >  1:1193038/1‑150 (MQ=255)
  aGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTTCGCCCTGGACAGCGGCCGTAATGATGACGAAGCCCTCa                                                                                                                                           >  1:1007023/1‑150 (MQ=255)
    gTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACGATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATg                                                                                                                                         >  2:1028698/1‑150 (MQ=255)
                  gCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCAt                                                                                                                                          >  2:188668/1‑135 (MQ=255)
                  gCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCAt                                                                                                                                          <  1:188668/135‑1 (MQ=255)
                    tCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAAc                                                                                                                         >  2:340043/1‑150 (MQ=255)
                      ggCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGaa                                                                                                                          >  2:282791/1‑147 (MQ=255)
                      ggCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGaa                                                                                                                          <  1:282791/147‑1 (MQ=255)
                       gCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCATGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTgg                                                                                                                      <  1:813858/150‑1 (MQ=255)
                        cTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGc                                                                                                                     <  1:224016/150‑1 (MQ=255)
                               gCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTGATGAGTACGAACCCGTCATGCTTGCTGCCAGCGAACTGGCCGATGCt                                                                                                              >  1:806540/1‑150 (MQ=255)
                                       tCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCaaa                                                                                                      <  1:171917/150‑1 (MQ=255)
                                       tCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCaaa                                                                                                      >  2:1092627/1‑150 (MQ=255)
                                            ggTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaa                                                                                                 <  2:1193038/150‑1 (MQ=255)
                                            ggTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaa                                                                                                 <  2:1391455/150‑1 (MQ=255)
                                             gTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaag                                                                                                >  2:921574/1‑150 (MQ=255)
                                              ttGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaagg                                                                                               >  1:650285/1‑150 (MQ=255)
                                                  aTNCACTCACGCGCAAGGTCTGGGGCCNAAACCACTGGCCNGCGGTCNTNCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGgcgc                                                                                           >  1:982891/1‑150 (MQ=255)
                                                  aTCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaagg                                                                                               <  1:1236923/146‑1 (MQ=255)
                                                  aTCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaagg                                                                                               >  2:1236923/1‑146 (MQ=255)
                                                   cccACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGgcgcn                                                                                          <  2:982891/149‑2 (MQ=255)
                                                      aCTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGt                                                                                       >  2:64931/1‑150 (MQ=255)
                                                       cTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTc                                                                                      >  1:853934/1‑150 (MQ=255)
                                                         cACGCGCAAGGTCTGCGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCACTGCACATCGGCCGTAAAGATGACGAACCAGTCATGCTCGATGCTAGAGAACTTGCTGATTCTGGCGCAAAGGGAAGGGCGCCCGTccc                                                                                    <  1:324504/150‑1 (MQ=255)
                                                          aCGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaagg                                                                                               >  2:1212713/1‑138 (MQ=255)
                                                          aCGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaagg                                                                                               <  1:1212713/138‑1 (MQ=255)
                                                           cgcgCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCgg                                                                                  <  1:845182/150‑1 (MQ=255)
                                                                  ggTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGa                                                                           <  1:64931/150‑1 (MQ=255)
                                                                    tCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATg                                                                         >  2:1095425/1‑150 (MQ=255)
                                                                      tGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATgta                                                                       >  2:1412884/1‑150 (MQ=255)
                                                                          gCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAg                                                                                                                              >  1:1278381/1‑91 (MQ=255)
                                                                          gCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAg                                                                                                                              <  2:1278381/91‑1 (MQ=255)
                                                                          gCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTAc                                                                   >  2:994224/1‑150 (MQ=255)
                                                                           ccAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  >  1:949283/1‑150 (MQ=255)
                                                                           ccAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  <  1:921574/150‑1 (MQ=255)
                                                                           ccAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  <  2:650285/150‑1 (MQ=255)
                                                                           ccAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  >  1:189211/1‑150 (MQ=255)
                                                                              aaaCCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCt                                                               <  1:86542/150‑1 (MQ=255)
                                                                              aaaCCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCt                                                               <  1:567637/150‑1 (MQ=255)
                                                                                  cACTGTCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGGCGCGGTCGTCATGATGACGAGCGCATACTGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGACACGGGACGGGCGCCCGTCCCGGCGTCGGGTGTAGTACGGCTGtgt                                                           >  1:443211/1‑148 (MQ=255)
                                                                                  cACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAGCCCATCATGCTCGCTGCCAGCGAACTGGc                                                                                                                     <  1:1308044/92‑1 (MQ=255)
                                                                                  cACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAGCCCATCATGCTCGCTGCCAGCGAACTGGc                                                                                                                     >  2:1308044/1‑92 (MQ=255)
                                                                                  cACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGttt                                                           <  1:1412884/150‑1 (MQ=255)
                                                                                  cACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCATCATGCTCGCTGCCAGCGAACTGGc                                                                                                                     >  2:594023/1‑92 (MQ=255)
                                                                                  cACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCATCATGCTCGCTGCCAGCGAACTGGc                                                                                                                     <  1:594023/92‑1 (MQ=255)
                                                                                    cTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTttct                                                         >  1:756670/1‑150 (MQ=255)
                                                                                     tGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTttctt                                                        <  1:706224/150‑1 (MQ=255)
                                                                                             ggTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGc                                                >  1:399943/1‑150 (MQ=255)
                                                                                                 gTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCGGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTGTTCGGTCCGCCTTc                                            >  1:482916/1‑150 (MQ=255)
                                                                                                   gCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCAt                                          >  1:870917/1‑150 (MQ=255)
                                                                                                       aTCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTTGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGAc                                      <  1:1095425/150‑1 (MQ=255)
                                                                                                       aTCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaa                                                                                                 <  1:1069135/91‑1 (MQ=255)
                                                                                                       aTCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAaagggaa                                                                                                 >  2:1069135/1‑91 (MQ=255)
                                                                                                          tACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGCTGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCa                                   >  2:165284/1‑150 (MQ=255)
                                                                                                            ccATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTcc                                                  >  2:659129/1‑133 (MQ=255)
                                                                                                            ccATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTcc                                                  <  1:659129/133‑1 (MQ=255)
                                                                                                             cATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTgg                                                                                                            >  1:33052/1‑74 (MQ=255)
                                                                                                             cATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTgg                                                                                                            <  2:33052/74‑1 (MQ=255)
                                                                                                             cATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTc                                >  2:933464/1‑150 (MQ=255)
                                                                                                               ttCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCAt                              >  2:586177/1‑150 (MQ=255)
                                                                                                                tCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATa                             <  1:586177/150‑1 (MQ=255)
                                                                                                                 ccGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATAc                            <  1:994224/150‑1 (MQ=255)
                                                                                                                 ccGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTAATGCACGCTGCCAGCGAACTGGCCGATGGTGGCGCAAAGTGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATAc                            <  2:1228901/150‑1 (MQ=255)
                                                                                                                      cTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCTAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGcc                       >  1:1038944/1‑150 (MQ=255)
                                                                                                                      cTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGcc                       >  1:976497/1‑150 (MQ=255)
                                                                                                                       tGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGGTGTAGTACGGCTGTTTTTTCGGTCCGCCTTCATTGAACCATTAATACCAGCCg                      >  2:186206/1‑150 (MQ=255)
                                                                                                                         ggCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTc                    <  2:1075402/150‑1 (MQ=255)
                                                                                                                           cATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGa                  <  2:399943/150‑1 (MQ=255)
                                                                                                                             tCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATc                <  2:756670/150‑1 (MQ=255)
                                                                                                                                      atgatgACGAACCCGTCATTCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCAGATGTAGTACGGCTTTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATCGGCGTTATg       <  1:186206/150‑1 (MQ=255)
                                                                                                                                      atgatgACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATCGGCGTTATg       <  1:470889/150‑1 (MQ=255)
                                                                                                                                       tgatgaCGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  <  2:429592/90‑1 (MQ=255)
                                                                                                                                       tgatgaCGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACg                                                                  >  1:429592/1‑90 (MQ=255)
                                                                                                                                           gaCGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATCGGCGTTATGCCTCg  <  1:1011840/150‑1 (MQ=255)
                                                                                                                                             cGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTttctt                                                        <  1:203746/94‑1 (MQ=255)
                                                                                                                                             cGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTttctt                                                        >  2:203746/1‑94 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTCACGCGCAAGGTCTGGGGCCAAAACCACTGGCCTGCGGTCGTGCACATCTACCATTCCGCCCTGGGCATCGGCCGTAATGATGACGAACCCGTCATGCTCGCTGCCAGCGAACTGGCCGATGCTGGCGCAAAGGGAAGGGCGCCCGTCCCGGCGTCGGATGTAGTACGGCTGTTTCTTCGGTCCGCCTTCATCGACCCATTCATACCAGCCGTCGATCGGCGTTATGCCTCG  >  NC_002947/3344481‑3344769

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