Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 258,642 A→G 7.8% K59R (AAG→AGG)  PP_0207 → nitrate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947258,6420AG 7.8% 229.0 / 7.3 77K59R (AAG→AGG) PP_0207nitrate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (3/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CTCGCGCTGGCGCTCACTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCTGGAT  >  NC_002947/258494‑258784
                                                                                                                                                     |                                                                                                                                              
cTCGCGCTGGCGCTCACTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAg                                                                                                                                               >  2:243519‑M2/1‑149 (MQ=255)
     gcTGGCGCTCACTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACaa                                                                                                                                          >  2:1355230/1‑150 (MQ=255)
         gcgcTCACTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGAc                                                                                                                                      <  1:1429492/150‑1 (MQ=255)
             tCACCGGCCTGAACGCCAGCGCCGAAACCATACGCAGCGCCTTCGGCACGCAGGGGACGACCATCAAGTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCa                                                                                                                                  <  2:513598/150‑1 (MQ=255)
               aCTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGTTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCAcc                                                                                                                                <  1:355669/150‑1 (MQ=255)
                     cTGCACGCACGCGCAGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAg                                                                                                                          <  2:93449/150‑1 (MQ=255)
                     cTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                >  2:441455/1‑128 (MQ=255)
                     cTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                <  1:441455/128‑1 (MQ=255)
                     cTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAg                                                                                                                          >  1:1374623/1‑150 (MQ=255)
                      tGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAgg                                                                                                                         <  2:1374623/150‑1 (MQ=255)
                            cccgcgcCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACCGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGt                                                                                                                   >  2:1574449/4‑150 (MQ=255)
                            ccAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTCCCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGt                                                                                                                   <  2:640355/150‑1 (MQ=255)
                                 gcCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGgaag                                                                                                              >  1:113764/1‑150 (MQ=255)
                                    gAAACCATACGCATCGCCATCGGCACCCAGGACACAACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAAc                                                                                                           <  1:555714/150‑1 (MQ=255)
                                           tACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                >  2:5055/1‑106 (MQ=255)
                                           tACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                <  1:5055/106‑1 (MQ=255)
                                               cATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCgg                                                                                                <  2:1447684/150‑1 (MQ=255)
                                                aTCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                >  1:50342/1‑101 (MQ=255)
                                                aTCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCaa                                                                                                                                                <  2:50342/101‑1 (MQ=255)
                                                aTCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGt                                                                                               >  1:93015/1‑150 (MQ=255)
                                                  cGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTgc                                                                                             <  2:93015/150‑1 (MQ=255)
                                                   gCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGgaag                                                                                                              >  2:586533/1‑132 (MQ=255)
                                                   gCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGgaag                                                                                                              <  1:586533/132‑1 (MQ=255)
                                                          gCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACTTGCCGTCGTGCGGCAGgtacttgtccagca                                                                                                                                  >  2:80197‑M2/1‑91 (MQ=255)
                                                          gCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACTTGCCGTCGTGCGGCAGgtacttgtccagca                                                                                                                                  <  1:80197‑M2/105‑15 (MQ=255)
                                                              ccAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGgaagaa                                                                                                            <  1:396383/123‑1 (MQ=255)
                                                              ccAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGgaagaa                                                                                                            >  2:396383/1‑123 (MQ=255)
                                                              ccAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCaa                                                                                 <  1:1294837/150‑1 (MQ=255)
                                                               cAGGACACCACCATCAACTGCGCCACCGGGGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAAc                                                                                >  2:1256164/1‑150 (MQ=255)
                                                                     accaccATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATg                                                                          >  1:804523/1‑150 (MQ=255)
                                                                                      cACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGa                                                         >  1:1213541/1‑150 (MQ=255)
                                                                                      cACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGa                                                         <  1:1355230/150‑1 (MQ=255)
                                                                                             ggCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGgc                                                  <  2:1597953/150‑1 (MQ=255)
                                                                                                ctgctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGcc                                               >  1:741810/1‑150 (MQ=255)
                                                                                                 tgctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTc                                                           <  1:708629/137‑1 (MQ=255)
                                                                                                 tgctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTc                                                           >  2:708629/1‑137 (MQ=255)
                                                                                                 tgctgATCCGTGAGCTGGGCCTGCTGAACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAg                                                                                                                    <  1:781069/80‑1 (MQ=255)
                                                                                                 tgctgATCCGTGAGCTGGGCCTGCTGAACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAg                                                                                                                    >  2:781069/1‑80 (MQ=255)
                                                                                                  gctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTAc                                                                                                                             <  2:362538/70‑1 (MQ=255)
                                                                                                  gctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTAc                                                                                                                             >  1:362538/1‑70 (MQ=255)
                                                                                                   ctgATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGACATGGTCGCCGGCAAGCTCGACTCCGGCGCCATg                                            <  2:1257289/150‑1 (MQ=255)
                                                                                                     gATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGc                                          <  2:1101157/150‑1 (MQ=255)
                                                                                                        ccGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGa                                       <  2:967702/150‑1 (MQ=255)
                                                                                                          gTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACt                                     >  1:790677/1‑150 (MQ=255)
                                                                                                           tGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACtt                                    >  1:1158580/1‑150 (MQ=255)
                                                                                                                 gggCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTATCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAgg                              >  1:1342773/1‑150 (MQ=255)
                                                                                                                         tgGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATg                                                                          <  1:1257998/98‑1 (MQ=255)
                                                                                                                         tgGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATg                                                                          >  2:1257998/1‑98 (MQ=255)
                                                                                                                           gACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGATTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCGCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    <  1:1574449/150‑1 (MQ=255)
                                                                                                                           gACAAGTACCTGCCGCACGACGGCAAACATAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    <  1:1236330/150‑1 (MQ=255)
                                                                                                                            aCAAGTACCTGCCGCACGACGGCAAgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatggtggtgtcct                                                                            <  2:695001‑M2/93‑69 (MQ=255)
                                                                                                                            aCAAGTACCTGCCGCACGACGGCAAgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatggtggtgtcct                                                                            >  1:695001‑M2/1‑25 (MQ=255)
                                                                                                                              aaGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCAt                         >  2:1356172/1‑142 (MQ=255)
                                                                                                                              aaGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCAt                         <  1:1356172/142‑1 (MQ=255)
                                                                                                                                gTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCAc                                                                                                      <  2:1510010/63‑1 (MQ=255)
                                                                                                                                gTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCAc                                                                                                      >  1:1510010/1‑63 (MQ=255)
                                                                                                                                  aCCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    <  2:1128681/143‑1 (MQ=255)
                                                                                                                                  aCCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    <  2:1128501/143‑1 (MQ=255)
                                                                                                                                  aCCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    >  1:1128681/1‑143 (MQ=255)
                                                                                                                                  aCCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAAc                    >  1:1128501/1‑143 (MQ=255)
                                                                                                                                   ccTGCCGCACGACGGCAGgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatg                                                                                      >  2:587460‑M2/1‑18 (MQ=255)
                                                                                                                                   ccTGCCGCACGACGGCAGgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatg                                                                                      <  1:587460‑M2/76‑59 (MQ=255)
                                                                                                                                   ccTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGa                                                                                     <  1:740705/77‑1 (MQ=255)
                                                                                                                                   ccTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGa                                                                                     >  2:740705/1‑77 (MQ=255)
                                                                                                                                    cTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTc                                                     <  2:968298/108‑1 (MQ=255)
                                                                                                                                    cTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTc                                                     >  1:968298/1‑108 (MQ=255)
                                                                                                                                    cTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGcc           >  1:614251/1‑150 (MQ=255)
                                                                                                                                    cTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCATCGGGATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGcc           >  1:1251428/1‑150 (MQ=255)
                                                                                                                                     tGCCGCACGACGGCAGgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatggtggtgtcctgggtgccgatggcgatgc                                                          <  2:462133‑M2/102‑87 (MQ=255)
                                                                                                                                     tGCCGCACGACGGCAGgtacttgtccagcaggcccagctcacggatcagcaggccaccggtggcgcagttgatggtggtgtcctgggtgccgatggcgatgc                                                          >  1:462133‑M2/1‑16 (MQ=255)
                                                                                                                                     tGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGcgcc                                                                                          >  2:1468156/1‑70 (MQ=255)
                                                                                                                                     tGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGcgcc                                                                                          <  1:1468156/70‑1 (MQ=255)
                                                                                                                                      gCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGaaa                                                                            >  2:1545512/1‑83 (MQ=255)
                                                                                                                                      gCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGaaa                                                                            <  1:1545512/83‑1 (MQ=255)
                                                                                                                                      gCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATg                                                                          <  1:855798/85‑1 (MQ=255)
                                                                                                                                      gCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATg                                                                          >  2:855798/1‑85 (MQ=255)
                                                                                                                                         gCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCt      >  1:656554/1‑150 (MQ=255)
                                                                                                                                           acgacgGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGtt              >  1:1116540/1‑140 (MQ=255)
                                                                                                                                           acgacgGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGtt              <  2:1116540/140‑1 (MQ=255)
                                                                                                                                            cgacgGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGACATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCTGGa   >  1:608802/1‑150 (MQ=255)
                                                                                                                                             gacgGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCGTTCAACGGCGTTGCCCATCTGGAt  <  2:790677/150‑1 (MQ=255)
                                                                                                                                                 ggCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCt      <  1:267017/143‑1 (MQ=255)
                                                                                                                                                 ggCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCt      >  2:267017/1‑143 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CTCGCGCTGGCGCTCACTGGCCTGAACGCCAGCGCCGAAACCATACGCATCGCCATCGGCACCCAGGACACCACCATCAACTGCGCCACCGGTGGCCTGCTGATCCGTGAGCTGGGCCTGCTGGACAAGTACCTGCCGCACGACGGCAAGTACAAGGACGCCACCTACCAGGTCGAGTGGAAGAACTTCACCAGCGGTGCGCCGCTGACCAACGAAATGGTCGCCGGCAAGCTCGACTTCGGCGCCATGGCCGACTTCCCAGGGTCATTCAACGGCGTTGCCCATCTGGAT  >  NC_002947/258494‑258784

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: