Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,389 C→T 6.0% intergenic (+417/‑103) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,3890CT 6.0% 280.8 / 4.7 94intergenic (+417/‑103)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base C (41/47);  new base T (3/3);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCGGCG  >  NC_002947/1899243‑1899534
                                                                                                                                                  |                                                                                                                                                 
ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCCCCGTCGAc                                                                                                                                                >  1:922920/1‑150 (MQ=255)
ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAc                                                                                                                                                <  1:150691/150‑1 (MQ=255)
 gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACa                                                                                                                                               >  2:701827/1‑150 (MQ=255)
    tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATGCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            >  1:1243524/1‑150 (MQ=255)
    tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            <  1:208737/150‑1 (MQ=255)
    tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            >  1:955013/1‑150 (MQ=255)
    tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            <  1:939116/150‑1 (MQ=255)
        cgcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTca                                                                                                                                        <  2:516764/150‑1 (MQ=255)
         gcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTGGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcac                                                                                                                                       <  1:1284444/150‑1 (MQ=255)
               aCGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGTCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTc                                                                                                                                 >  1:249182/1‑150 (MQ=255)
                       tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                                                                                         >  2:505496‑M2/55‑150 (MQ=255)
                       tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                                                                                         >  2:855685‑M2/55‑150 (MQ=255)
                       TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAgcccg                                                                                                                                          >  1:1147881‑M1/1‑132 (MQ=255)
                       TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAgcccg                                                                                                                                          <  2:1147881‑M1/133‑2 (MQ=255)
                       TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                                                                                         <  1:1019109‑M1/150‑1 (MQ=255)
                        cccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                                                                                        >  2:1421273‑M2/54‑150 (MQ=255)
                        CCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                                                                                        >  2:955013‑M1/1‑150 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       >  1:1153826‑M2/53‑150 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       >  2:308284‑M2/53‑150 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                                                                                                >  2:106841‑M2/53‑141 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                                                                                                <  1:106841‑M2/89‑1 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       >  2:1001149‑M2/53‑150 (MQ=255)
                         ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       >  2:1119965‑M2/53‑150 (MQ=255)
                         CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       <  1:57793‑M1/150‑1 (MQ=255)
                         CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                                                                                       <  1:412811‑M1/150‑1 (MQ=255)
                          cacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcg                                                                                                                     >  1:187775‑M2/52‑150 (MQ=255)
                           acagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTGCCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcg                                                                                                                     >  2:406547‑M2/51‑150 (MQ=255)
                                  cTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                                                                                              >  2:93056/1‑150 (MQ=255)
                                      gccgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGccc                                                                                                          >  1:1239510/1‑150 (MQ=255)
                                       caaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGcccc                                                                                                         >  1:286496‑M2/39‑150 (MQ=255)
                                        aaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCg                                                                                                        >  2:1160777‑M2/38‑150 (MQ=255)
                                        aaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCg                                                                                                        >  1:831982‑M2/38‑150 (MQ=255)
                                           ccTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCtt                                                                                                     <  2:781260/150‑1 (MQ=255)
                                            cTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTg                                                                                                    >  2:368548/1‑150 (MQ=255)
                                              tgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCa                                                                                                  <  2:1153826‑M2/119‑1 (MQ=255)
                                               aGGCGTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTCGACAGCCCTCACCACTTCGCGACAAAGGCGCCGTGTTGCCCCGCTTGCac                                                                                                 >  2:94844/1‑149 (MQ=255)
                                                 gggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCa                                                                                               <  2:286496‑M2/122‑1 (MQ=255)
                                                         atatCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCg                                                                                       <  1:1281238/150‑1 (MQ=255)
                                                            tacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgggg                                                                                    <  1:308284‑M2/133‑1 (MQ=255)
                                                             cccTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGtgttgccccggcatgctgatgcaagcggggc                                                                                   <  2:1353443‑M1/150‑32 (MQ=255)
                                                             acctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGc                                                                                   <  1:1001149‑M2/134‑1 (MQ=255)
                                                              cctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                                                                                   <  1:121934‑M2/102‑1 (MQ=255)
                                                              cctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                                                                                   >  2:121934‑M2/16‑117 (MQ=255)
                                                                  tgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCg                                                                              <  1:141410‑M2/139‑1 (MQ=255)
                                                                  tgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCg                                                                              >  1:1318580‑M2/12‑150 (MQ=255)
                                                                   tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                                                                                               >  1:364561‑M2/11‑100 (MQ=255)
                                                                   tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                                                                                               <  2:364561‑M2/90‑1 (MQ=255)
                                                                   tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                                                                                     <  2:1338091‑M2/132‑1 (MQ=255)
                                                                   tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                                                                                     >  1:1338091‑M2/11‑142 (MQ=255)
                                                                    ggAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGtgttgccctggcatgctgatgcaagcggggcaacacgg                                                                            <  1:970239‑M1/150‑39 (MQ=255)
                                                                    ggAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggggcaacacggc                                                                            <  2:355134‑M1/150‑39 (MQ=255)
                                                                          ggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttt                                                                      >  2:585951‑M2/4‑150 (MQ=255)
                                                                           gcgtGCCCGCGAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTGAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttttt                                                                    <  2:187775/147‑1 (MQ=255)
                                                                            cATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttttt                                                                    <  1:1300942‑M2/149‑1 (MQ=255)
                                                                                cccGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGt                                                                <  2:457690/150‑1 (MQ=255)
                                                                                        gaggCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCcaca                                                        <  2:831982/148‑1 (MQ=255)
                                                                                         gggCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAg                                                       <  1:428168/150‑1 (MQ=255)
                                                                                             cGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAgc                                                   <  2:1243524/150‑1 (MQ=255)
                                                                                                tataGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAgcccg                                                >  2:1103995/1‑150 (MQ=255)
                                                                                                       ggCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATgg                                         <  1:1421273/150‑1 (MQ=255)
                                                                                                          ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggggcaacacggc                                                                            <  2:1024373‑M1/112‑39 (MQ=255)
                                                                                                          ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggggcaacacggc                                                                            >  1:1024373‑M1/1‑74 (MQ=255)
                                                                                                           cATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCcac                                                         >  2:665639/1‑130 (MQ=255)
                                                                                                           cATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCcac                                                         <  1:665639/130‑1 (MQ=255)
                                                                                                            aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCtttt                                                                     >  1:120082/1‑117 (MQ=255)
                                                                                                            aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCtttt                                                                     <  2:120082/117‑1 (MQ=255)
                                                                                                             tAAATTGCTGATCTGCCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACa                                   <  1:406547/150‑1 (MQ=255)
                                                                                                             tAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACa                                   >  2:647132/1‑150 (MQ=255)
                                                                                                               aaTTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGc                                 <  2:1019109/150‑1 (MQ=255)
                                                                                                                  tGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGaa                              <  1:368548/150‑1 (MQ=255)
                                                                                                                  tGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGaa                              <  2:324194/150‑1 (MQ=255)
                                                                                                                  tGCTGATCTACCTCATCCCGGCCTGTCAACCTCGACAGACCTACCCACTTCACGACCAAGGCGCCGTGTTGCCGCGCTGGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGaa                              <  1:94844/150‑1 (MQ=255)
                                                                                                                       aTCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttttt                                                                    >  1:1145051/1‑107 (MQ=255)
                                                                                                                       aTCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttttt                                                                    <  2:1145051/107‑1 (MQ=255)
                                                                                                                       aTCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACg                         >  2:1352513/1‑150 (MQ=255)
                                                                                                                           aCCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGc                     >  1:1173490/1‑150 (MQ=255)
                                                                                                                              tCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCcaca                                                        <  2:232132/112‑1 (MQ=255)
                                                                                                                              tCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCcaca                                                        >  1:232132/1‑112 (MQ=255)
                                                                                                                              tCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTc                  <  2:1318580/150‑1 (MQ=255)
                                                                                                                                 tCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGc               <  1:855685/150‑1 (MQ=255)
                                                                                                                                  cccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGcc              <  1:93056/150‑1 (MQ=255)
                                                                                                                                   ccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGcccc                                                                                                         >  2:193506/1‑58 (MQ=255)
                                                                                                                                   ccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGcccc                                                                                                         <  1:193506/58‑1 (MQ=255)
                                                                                                                                   ccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGccc             <  1:647132/150‑1 (MQ=255)
                                                                                                                                      ggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgg                                                                                      >  1:1230537/1‑74 (MQ=39)
                                                                                                                                      ggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgg                                                                                      <  2:1230537/74‑1 (MQ=255)
                                                                                                                                      ggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCaa                                     >  1:1216496/1‑123 (MQ=255)
                                                                                                                                      ggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCaa                                     <  2:1216496/123‑1 (MQ=255)
                                                                                                                                          gTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgc      <  1:505496/150‑1 (MQ=255)
                                                                                                                                           tCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCcggcgcct                                                                        <  1:1333791‑M2/83‑9 (MQ=255)
                                                                                                                                           tCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCcggcgcct                                                                        >  2:1333791‑M2/1‑75 (MQ=255)
                                                                                                                                              aCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggg                                                                                      <  2:1289584‑M1/66‑29 (MQ=255)
                                                                                                                                              aCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggg                                                                                      >  1:1289584‑M1/1‑38 (MQ=255)
                                                                                                                                              aCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgcggcg  >  2:978458/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCGGCG  >  NC_002947/1899243‑1899534

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: