Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,458,128 C→A 6.6% intergenic (+21/+19) PP_1276 → / ← algA hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,458,1280CA 6.6% 280.3 / 5.1 89intergenic (+21/+19)PP_1276/algAhypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base A (3/3);  total (46/43)
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.58e-01

CTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGCG  >  NC_002947/1457979‑1458277
                                                                                                                                                     |                                                                                                                                                      
cTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGc                                                                                                                                                        <  1:1507218/150‑1 (MQ=255)
 tCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGcc                                                                                                                                                       <  2:356229/150‑1 (MQ=255)
   gAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCa                                                                                                                                                     >  1:1364874/1‑150 (MQ=255)
     aaTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAgg                                                                                                                                                   <  2:1980216/150‑1 (MQ=255)
     aaTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAgg                                                                                                                                                   <  1:643429/150‑1 (MQ=255)
      aTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAggg                                                                                                                                                  >  1:435206/1‑150 (MQ=255)
          cTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                              <  1:1770533/150‑1 (MQ=255)
          cTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                              >  2:1290805/1‑150 (MQ=255)
             ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTc                                                                                                                                           >  2:2138297/1‑150 (MQ=255)
                caAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGc                                                                                                                                        >  1:1161298/1‑150 (MQ=255)
                      tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTAc                                                                                                                                  >  2:394059/1‑150 (MQ=255)
                          ggCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTgcg                                                                                                                              >  2:319960/1‑150 (MQ=255)
                               aaTCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTCTCTTATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATg                                                                                                                         <  1:1845021/150‑1 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTTCGCGATGGTCtt                                                                                                                    >  2:645129/1‑150 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCtt                                                                                                                    <  2:504985/150‑1 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCtt                                                                                                                    >  1:922926/1‑150 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGGCAGCATTTACTGCGCGATGGTCtt                                                                                                                    >  1:443874/1‑150 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGAAAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCtt                                                                                                                    >  1:1450798/1‑150 (MQ=255)
                                      cGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGATGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCa                                                                                                                  <  2:1864641/150‑1 (MQ=255)
                                             gAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑Ac                                                                                                          <  1:2138297/150‑1 (MQ=255)
                                              aaGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCACTGGACAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACGGAA‑Acg                                                                                                         >  2:96666/1‑150 (MQ=255)
                                                ggggAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                              >  1:2352880/1‑112 (MQ=255)
                                                ggggAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                              <  2:2352880/112‑1 (MQ=255)
                                                         gatTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTAc                                                                                                                                  <  1:551125/115‑1 (MQ=255)
                                                         gatTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTAc                                                                                                                                  >  2:551125/1‑115 (MQ=255)
                                                         gatTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGAt                                                                                              >  1:320267/1‑150 (MQ=255)
                                                         gatTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGAt                                                                                              <  2:922926/150‑1 (MQ=255)
                                                           tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgg                                                                                            >  2:1426145/1‑150 (MQ=255)
                                                           tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgg                                                                                            <  2:1161298/150‑1 (MQ=255)
                                                             tGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          <  1:1220703/150‑1 (MQ=255)
                                                             tGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          <  1:784833/150‑1 (MQ=255)
                                                                 gCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCGCGCCCCTGGACAGCCCAGCGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTgg                                                                                      <  2:1040457/150‑1 (MQ=255)
                                                                  cGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGt                                                                                                                       >  2:1117346/1‑117 (MQ=255)
                                                                  cGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGt                                                                                                                       <  1:1117346/117‑1 (MQ=255)
                                                                   gACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGAc                                                                                    >  1:1229782/1‑150 (MQ=255)
                                                                    aCCGTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACg                                                                                   >  2:993059/1‑150 (MQ=255)
                                                                     ccTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGt                                                                                  <  2:168333/150‑1 (MQ=255)
                                                                     ccTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGt                                                                                  <  1:993059/150‑1 (MQ=255)
                                                                      cTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTg                                                                                 >  1:98123/1‑150 (MQ=255)
                                                                       tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGGCCCGAGCGTTGAATAAGCAGGCTCGCCCCGGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTgc                                                                                >  2:482680/1‑150 (MQ=255)
                                                                       tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATCAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTgc                                                                                >  2:1253931/1‑150 (MQ=255)
                                                                       tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACATgc                                                                                >  2:1404548/1‑150 (MQ=255)
                                                                            gggTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCATGGGCGAGTCAGCATAAACTCCGCGGCGG‑CATGATGTAAGATGCCGCGCTCGATGGGGGTGGACGTGCGCCCa                                                                           >  1:2377118/1‑150 (MQ=255)
                                                                              gtggtgGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATa                                                                         >  1:406489/1‑150 (MQ=255)
                                                                                ggtgggtgGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCAcc                                                                                                                <  1:1174336/110‑1 (MQ=255)
                                                                                ggtgggtgGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCAcc                                                                                                                >  2:1174336/1‑110 (MQ=255)
                                                                                 gtgggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATacac                                                                      <  1:2186599/150‑1 (MQ=255)
                                                                                 gtgggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATacac                                                                      >  2:2186599/1‑150 (MQ=255)
                                                                                 gtgggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAAGAAGAGCTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGCTGGGGGTGGACGTGCGCCCATacac                                                                      >  1:644111/1‑150 (MQ=255)
                                                                                    ggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATc                                                                   <  1:1253931/150‑1 (MQ=255)
                                                                                       ggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTc                                                                >  1:1079318/1‑150 (MQ=255)
                                                                                         cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          <  2:525555/122‑1 (MQ=255)
                                                                                         cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          >  1:525555/1‑122 (MQ=255)
                                                                                         cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGa                                                              >  2:413499/1‑150 (MQ=255)
                                                                                         cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGa                                                              >  2:2013526/1‑150 (MQ=255)
                                                                                          ccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  1:319960/150‑1 (MQ=255)
                                                                                          ccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  1:1290805/150‑1 (MQ=255)
                                                                                           cATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑Acgccg                                                                                                      >  2:1036659/1‑108 (MQ=255)
                                                                                           cATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑Acgccg                                                                                                      <  1:1036659/108‑1 (MQ=255)
                                                                                             ttCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTgc                                                                                <  1:1724353/128‑1 (MQ=255)
                                                                                             ttCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTgc                                                                                >  2:1724353/1‑128 (MQ=255)
                                                                                             gatgACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAgcg                                                          >  2:1559739/4‑150 (MQ=255)
                                                                                                 aCATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCg                                                      >  1:482652/1‑150 (MQ=255)
                                                                                                   aTCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCGCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                    >  1:1009562/1‑150 (MQ=255)
                                                                                                      tACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  2:751256/138‑1 (MQ=255)
                                                                                                      tACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             >  1:751256/1‑138 (MQ=255)
                                                                                                         gTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATc                                              <  1:1559739/150‑1 (MQ=255)
                                                                                                           cATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCt                                            <  1:394059/150‑1 (MQ=255)
                                                                                                                      gCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAg                                 <  2:406489/150‑1 (MQ=255)
                                                                                                                            aaTAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATg                                                   <  2:1612060/126‑1 (MQ=255)
                                                                                                                            aaTAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATg                                                   >  1:1612060/1‑126 (MQ=255)
                                                                                                                               aaGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGAc                        <  2:320267/150‑1 (MQ=255)
                                                                                                                                  cAGGCTCGCCCGTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGc                     <  2:435206/150‑1 (MQ=255)
                                                                                                                                   tgacTCGCCCCTGGGCTGACCAGGGGCGAGCCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             >  1:566962/4‑109 (MQ=38)
                                                                                                                                   tgacTCGCCCCTGGGCTGACCAGGGGCGAGCCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  2:566962/106‑1 (MQ=38)
                                                                                                                                   tgacTCGCCCCTGGGCTGACCAGGGGCAAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgg                                                                                            <  1:106973/75‑1 (MQ=25)
                                                                                                                                   tgacTCGCCCCTGGGCTGACCAGGGGCAAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgg                                                                                            >  2:106973/4‑78 (MQ=25)
                                                                                                                                   aGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCa                    >  2:1415939/1‑150 (MQ=255)
                                                                                                                                   aGGCGCGACCCGGGTCAGCCCAGCGGCGAGTCAGCATTTACTGCGTGAAGGACTTAGATAAG‑ATGTTATTTGTCATCCGGGTGGACTTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCa                    <  1:542930/150‑1 (MQ=255)
                                                                                                                                     gCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                    <  2:228740/116‑1 (MQ=255)
                                                                                                                                     gCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                    >  1:228740/1‑116 (MQ=255)
                                                                                                                                      cTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCt                 <  1:1415939/150‑1 (MQ=255)
                                                                                                                                       tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  2:1426376/105‑1 (MQ=255)
                                                                                                                                       tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             >  1:1426376/1‑105 (MQ=255)
                                                                                                                                           cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGa                         >  1:1229579/1‑137 (MQ=255)
                                                                                                                                           cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGa                         <  2:1229579/137‑1 (MQ=255)
                                                                                                                                           cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga              <  1:2218126/148‑1 (MQ=255)
                                                                                                                                           cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga              >  2:2218126/1‑148 (MQ=255)
                                                                                                                                            ccTGGTCAGCCCAGGGGCGAGTCGGCATTTACTGCGCGGTGGTCTTCACCGAA‑ACGCCGCGCTCGATATGACTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgatga           <  1:413499/150‑1 (MQ=255)
                                                                                                                                               ggTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGAtct        <  2:1364874/150‑1 (MQ=255)
                                                                                                                                                  cAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCa     <  2:2065014/150‑1 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          >  2:1892371/1‑62 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATgggg                                                                                          <  1:1892371/62‑1 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             <  1:149336/91‑1 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                             >  2:149336/1‑91 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat             >  1:579345/1‑139 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat             <  2:579345/139‑1 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGCg  >  2:392649/1‑150 (MQ=255)
                                                                                                                                                     cccAGGGGCGAGTCAGCATTTACTGCGCGATGGCTGATACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGCg  >  1:1145628/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAA‑ACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGCG  >  NC_002947/1457979‑1458277

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: ‑