Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,457,380 A→G 5.8% A136A (GCA→GCG)  PP_2150 → thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,457,3800AG 5.8% 310.1 / 6.2 106A136A (GCA→GCG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (50/47);  new base G (3/3);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 109 total aligned reads displayed.

CTGCCGCCTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGCCATGG  >  NC_002947/2457231‑2457528
                                                                                                                                                     |                                                                                                                                                     
ctGCCGCCTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCa                                                                                                                                                       <  1:2075440/150‑1 (MQ=255)
 tGCCGCCTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGTTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCAc                                                                                                                                                      >  1:70766/1‑150 (MQ=255)
     gccTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGc                                                                                                                                                  >  2:2100698/1‑150 (MQ=255)
       cTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGccc                                                                                                                                                <  1:444777/150‑1 (MQ=255)
               cTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTgcgc                                                                                                                                        >  1:2299070/1‑150 (MQ=255)
                tCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTgg                                                                                                                                                   >  1:924045/1‑138 (MQ=255)
                tCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTgg                                                                                                                                                   <  2:924045/138‑1 (MQ=255)
                  gcgTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAg                                                                                                                                     >  1:457095/1‑150 (MQ=255)
                   cccatggcctgcaggcgcgcagcgatcaagtcgacagcatgaccggcgacgatgctgttgaagtcgagctgcaccggggcatccttgcacagggctcgacctttgacgcgcaagtgctggtaacccacaCACTGGCCCAGGTGCGCCAGc                                                                                                                                    >  2:192162‑M1/130‑150 (MQ=255)
                    gTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCg                                                                                                                                   >  1:724184/1‑150 (MQ=255)
                         tCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTgg                                                                                                                              <  2:568239/150‑1 (MQ=255)
                         tCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTgg                                                                                                                              <  1:191624/150‑1 (MQ=255)
                         tCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTgg                                                                                                                              <  2:568573/150‑1 (MQ=255)
                          cAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTggg                                                                                                                             >  1:266322/1‑150 (MQ=255)
                                  gTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCCGc                                                                                                                     >  1:2155033/1‑150 (MQ=255)
                                    cAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTgc                                                                                                                                          >  2:2257986/1‑127 (MQ=255)
                                      gCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGt                                                                                                                                            <  2:899698/123‑1 (MQ=255)
                                      gCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGt                                                                                                                                            >  1:899698/1‑123 (MQ=255)
                                        aCCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTgc                                                                                                               <  2:1556035/150‑1 (MQ=255)
                                        aCCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACGTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTgc                                                                                                               >  2:1564105/1‑150 (MQ=255)
                                         ccTTGCCGAGTAGATCGATGGCGCTATCGATCTCACCGTAAACACATTGCTCGACATGTGGGGCTTCGGCCCTCCGGCCCGTCACAAGCAAGTGCCCGAACCACAGGCCCAGGCCCAGGTGCTCCCGCGTTTGCTTTCCCAGCACCTgcg                                                                                                              >  2:1159893/1‑150 (MQ=255)
                                         ccTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCa                                                                                                                                                       >  1:780995/1‑109 (MQ=255)
                                         ccTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCa                                                                                                                                                       <  2:780995/109‑1 (MQ=255)
                                         ccTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTgcg                                                                                                              >  1:1301822/1‑150 (MQ=255)
                                          cTTGCCGAGCAGAGGGGTGGCGCTTTCGAGCGCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTgcgc                                                                                                             <  2:756350/150‑1 (MQ=255)
                                          cTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTgcgc                                                                                                             >  2:2198204/1‑150 (MQ=255)
                                           ttGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTGACCAGCACTTGcgcg                                                                                                            <  1:2373810/150‑1 (MQ=255)
                                             gCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTc                                                                                                          <  2:1723490/150‑1 (MQ=255)
                                               cGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGTCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCaa                                                                                                        >  1:759099/1‑150 (MQ=255)
                                                gAGCAGAGCGATGGCGTTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGAGCGTCaaa                                                                                                       <  2:70766/150‑1 (MQ=255)
                                                gAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCaaa                                                                                                       <  2:1988274/150‑1 (MQ=255)
                                                  gCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAgg                                                                                                     >  1:2406188/1‑150 (MQ=255)
                                                  gCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCACGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAATTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCAGTTGCGCGTCAAAgg                                                                                                     >  2:226445/1‑150 (MQ=255)
                                                     gagCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCgtgt                                                                                                                                <  2:832926/120‑1 (MQ=255)
                                                                 tttCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCaa                                                                                      >  2:1505283/1‑150 (MQ=255)
                                                                  ttCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAg                                                                                     <  1:1260362/150‑1 (MQ=255)
                                                                  ttCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAg                                                                                     >  2:1095806/1‑150 (MQ=255)
                                                                      aTCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCa                                                                                                                    <  1:434956/115‑1 (MQ=255)
                                                                      aTCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCa                                                                                                                    >  2:434956/1‑115 (MQ=255)
                                                                       tctcACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGCCGAGCCCTGTGCAAGGATGc                                                                                >  2:1444880/1‑150 (MQ=255)
                                                                         tcACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGccc                                                                              >  1:2215508/1‑150 (MQ=255)
                                                                                 gagCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAg                                                                      <  2:639851/150‑1 (MQ=255)
                                                                                      atttcgcGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGa                                                                 >  2:836975/7‑150 (MQ=255)
                                                                                          cTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCa                                                                                                                       >  2:1201297/1‑92 (MQ=255)
                                                                                          cTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCa                                                                                                                       <  1:1201297/92‑1 (MQ=255)
                                                                                           tCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCt                                                                                            <  1:2292254/118‑1 (MQ=255)
                                                                                           tCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCt                                                                                            >  2:2292254/1‑118 (MQ=255)
                                                                                           tCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCa                                                           >  1:2370480/1‑150 (MQ=255)
                                                                                           tCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCa                                                           >  2:451072/1‑150 (MQ=255)
                                                                                           accggggcatccttgcacagggctcgacctttgacgcgcaagtgctggtaacccacaCGCTGGCCCAGGTGCGCCAg                                                                                                                                     <  2:1901372‑M1/20‑1 (MQ=255)
                                                                                           accggggcatccttgcacagggctcgacctttgacgcgcaagtgctggtaacccacaCGCTGGCCCAGGTGCGCCAg                                                                                                                                     >  1:1901372‑M1/58‑77 (MQ=255)
                                                                                             gACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAAc                                                         >  2:2129793/1‑150 (MQ=255)
                                                                                              aCTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACa                                                        <  1:451072/150‑1 (MQ=255)
                                                                                                ttGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGc                                                      >  2:1543585/1‑150 (MQ=255)
                                                                                                   tgGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTc                                                                                                   >  2:408951/1‑103 (MQ=255)
                                                                                                   tgGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTc                                                                                                   <  1:408951/103‑1 (MQ=255)
                                                                                                     ggtcTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCAtcgt                                                 >  1:85406/4‑150 (MQ=255)
                                                                                                     gggCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGTTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCAtcgt                                                 <  1:291751/150‑1 (MQ=255)
                                                                                                     gggCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCAtcgt                                                 >  1:2284425/1‑150 (MQ=255)
                                                                                                         ttCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TAAACAGCATCGTcccc                                             >  2:808198/1‑150 (MQ=255)
                                                                                                             gCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTc                                         <  1:1564105/150‑1 (MQ=255)
                                                                                                             gCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTc                                         >  1:98032/1‑150 (MQ=255)
                                                                                                             gCCCTCAGGCCCGTCACGAGCAAGCGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTc                                         <  2:2096455/150‑1 (MQ=255)
                                                                                                                cTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATg                                      <  2:2155033/150‑1 (MQ=255)
                                                                                                                   aGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCa                                                                                                                       <  2:2353646/67‑1 (MQ=255)
                                                                                                                   aGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCa                                                                                                                       >  1:2353646/1‑67 (MQ=255)
                                                                                                                     gCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGccc                                                                              <  2:1668150/106‑1 (MQ=255)
                                                                                                                     gCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGccc                                                                              >  1:1668150/1‑106 (MQ=255)
                                                                                                                       ccGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCgc                               >  1:2390947/1‑149 (MQ=255)
                                                                                                                       ccGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGa                               <  2:559238/150‑1 (MQ=255)
                                                                                                                        cGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGAc                              <  1:37023/150‑1 (MQ=255)
                                                                                                                          tCACGAGCAAGTGCCCGAACCACAGGCCCTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGGGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACtt                            <  1:1582735/150‑1 (MQ=255)
                                                                                                                          tCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACtt                            <  1:1543585/150‑1 (MQ=255)
                                                                                                                          tCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACtt                            <  2:311674/150‑1 (MQ=255)
                                                                                                                          tCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACtt                            <  2:105694/150‑1 (MQ=255)
                                                                                                                            aCGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATCCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGa                          <  2:759099/150‑1 (MQ=255)
                                                                                                                              gAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGa                                                                 >  2:2077003/1‑110 (MQ=255)
                                                                                                                              gAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGa                                                                 <  1:2077003/110‑1 (MQ=255)
                                                                                                                              gAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATc                        <  2:925973/150‑1 (MQ=255)
                                                                                                                              gAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATc                        >  1:1828921/1‑150 (MQ=255)
                                                                                                                               aGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCg                       <  2:2236153/150‑1 (MQ=255)
                                                                                                                                gcaagtgctggtaacccacaCGCTGGCGCACCTGGGCCAGTGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCaa                                                                                      >  2:1492851‑M1/21‑87 (MQ=255)
                                                                                                                                gcaagtgctggtaacccacaCGCTGGCGCACCTGGGCCAGTGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCaa                                                                                      <  1:1492851‑M1/67‑1 (MQ=255)
                                                                                                                                    gTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTgcg                  <  2:85406/150‑1 (MQ=255)
                                                                                                                                       cccGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTgcgcgc               >  2:1637227/1‑150 (MQ=255)
                                                                                                                                        ccGAACCACAGGCACTGGCAAAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGcc              <  1:226445/150‑1 (MQ=255)
                                                                                                                                          gAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCgg                                           >  1:995409/1‑119 (MQ=255)
                                                                                                                                            aCCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCa          <  2:2406188/150‑1 (MQ=255)
                                                                                                                                             ccACAGGCAATGCCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCACTGTGCAAGCATGACCCGGTGCAGCTCGTCT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAg         <  1:1099715/150‑1 (MQ=255)
                                                                                                                                              cacaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAgg        <  1:1095806/150‑1 (MQ=255)
                                                                                                                                              cacaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAgg        >  1:1778816/1‑150 (MQ=255)
                                                                                                                                              cacaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAg         <  1:1985748/149‑1 (MQ=255)
                                                                                                                                               ccacaCGCTGGCGCACCTGGGCCAGTGTGTGGGTTACCAGTACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAg                                                                      <  2:343610‑M1/83‑1 (MQ=255)
                                                                                                                                               ccacaCGCTGGCGCACCTGGGCCAGTGTGTGGGTTACCAGTACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAg                                                                      >  1:343610‑M1/6‑88 (MQ=255)
                                                                                                                                               acaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCCAGGTTGTCACGGTGCAAC‑CGTTTATCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGc       <  1:2100698/150‑1 (MQ=255)
                                                                                                                                               acaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGc       >  1:2160853/1‑150 (MQ=255)
                                                                                                                                               acaGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGc       >  1:1701055/1‑150 (MQ=255)
                                                                                                                                                  ggCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGCCAt    <  1:836975/150‑1 (MQ=255)
                                                                                                                                                  ggCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGCCAt    >  1:2081610/1‑150 (MQ=255)
                                                                                                                                                    CACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGCCATgg  >  2:1301822‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCCGCCTGACATGCTCGCGTTGGTCAGCCTGGGTCAGCACCTTGCCGAGCAGAGCGATGGCGCTTTCGATCTCACCGTAGAGCCATTGCTCGACTTGTGGGGCTTCGGCCCTCAGGCCCGTCACGAGCAAGTGCCCGAACCACAGGCACTGGCCCAGGTGCGCCAGCGTGTGGGTTACCAGCACTTGCGCGTCAAAGGTCGAGCCCTGTGCAAGGATGCCCCGGTGCAGCTCGACT‑TCAACAGCATCGTCGCCGGTCATGCTGTCGACTTGATCGCTGCGCGCCTGCAGGCCATGG  >  NC_002947/2457231‑2457528

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