Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,839 T→G 9.1% intergenic (‑70/+10) recA ← / ← pncC DNA recombination/repair protein RecA/NMN aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,8390TG 9.1% 228.2 / 14.6 77intergenic (‑70/+10)recA/pncCDNA recombination/repair protein RecA/NMN aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (38/32);  new base G (3/4);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

GGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAG  >  NC_000913/2823692‑2823987
                                                                                                                                                   |                                                                                                                                                     
ggAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCa                                                                                                                                                     >  1:1172807/1‑150 (MQ=255)
   gCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                  >  2:799562/1‑150 (MQ=255)
   gCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATAc                                                                                                                                                  <  1:841853/150‑1 (MQ=255)
       ttACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTa                                                                                                                                              >  1:111659/1‑150 (MQ=255)
       ttACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTa                                                                                                                                              >  1:1347024/1‑150 (MQ=255)
        tACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTAt                                                                                                                                             >  1:1045514/1‑150 (MQ=255)
         aCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAta                                                                                                                                            <  1:684019/150‑2 (MQ=255)
         aCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGGCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCTTACAGTAta                                                                                                                                            <  2:1423094/150‑2 (MQ=255)
          ccAAATTGTTTCTCAAGCTGGCCCAGTGCTGCCGCCCACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGGTCTGGTGAAGCAATTATACTGTATGCTCATACAGTATCa                                                                                                                                           >  1:1409041/1‑150 (MQ=255)
            aaaTTTTTTTTCCATCTGGCCCCTTTGCGGCGCCCAAGCGTTTTTTTTTTTTTCGTCGTTTCCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAg                                                                                                                                         <  1:1243607/150‑1 (MQ=255)
            aaaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGAGCATACAGTAtaatt                                                                                                                                         <  2:1347024/150‑5 (MQ=255)
             aaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                        >  2:1361044/1‑150 (MQ=255)
                tGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtt                                                                                                                                     <  2:39290/150‑1 (MQ=255)
                   ttCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTg                                                                                                                                  <  1:927743/150‑1 (MQ=255)
                    tctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGt                                                                                                                                 >  1:438722/1‑150 (MQ=255)
                     ctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                                                                                >  2:951521/1‑150 (MQ=255)
                     ctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTa                                                                                                                                >  2:591744/1‑150 (MQ=255)
                        aaTCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGaa                                                                                                                             <  1:1361044/150‑1 (MQ=255)
                             ggCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                        <  1:591744/150‑1 (MQ=255)
                               cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                      >  1:1196471/1‑150 (MQ=255)
                               cccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgttg                                                                                                                      <  2:1045514/150‑1 (MQ=255)
                                 cAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                                                                    <  1:781478/150‑1 (MQ=255)
                                  aGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCca                                                                                                                   >  1:589399/1‑150 (MQ=255)
                                   gTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcac                                                                                                                  <  1:1154493/150‑1 (MQ=255)
                                     gctgcCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACaa                                                                                                                >  2:1332469/1‑150 (MQ=255)
                                                   tttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCAt                                                                                                  <  2:728230/150‑1 (MQ=255)
                                                       tgtttgttTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGc                                                                                              <  1:320232/150‑1 (MQ=255)
                                                        gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCa                                                                                             >  1:1075048/1‑150 (MQ=255)
                                                          ttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGt                                                                                           >  2:531011/1‑150 (MQ=255)
                                                              tttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGcc                                                                                       >  1:1111138/1‑150 (MQ=255)
                                                               ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      <  1:110383/150‑1 (MQ=255)
                                                               ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCt                                                                                      >  2:131706/1‑150 (MQ=255)
                                                                 cgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTga                                                                                    <  1:686200/150‑1 (MQ=255)
                                                                     gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacgac                                                                                <  1:1065471/150‑1 (MQ=255)
                                                                      aTAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGacgacg                                                                               <  2:1172807/150‑1 (MQ=255)
                                                                      aTAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGacgacg                                                                               <  2:1120242/150‑1 (MQ=255)
                                                                      aTAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGacgacg                                                                               <  1:1307534/150‑1 (MQ=255)
                                                                       tAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGc                                                                              <  1:131706/150‑1 (MQ=255)
                                                                         gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                                                                            >  1:417480/1‑150 (MQ=255)
                                                                         gCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAc                                                                            <  2:927050/150‑1 (MQ=255)
                                                                          ccATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCAcc                                                                           >  2:861135/1‑150 (MQ=255)
                                                                             tttttACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCa                                                                        <  2:257073/150‑1 (MQ=255)
                                                                                   cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                                  <  1:1097817/150‑1 (MQ=255)
                                                                                   cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACgg                                                                  <  1:1294576/150‑1 (MQ=255)
                                                                                     ccTGTCATGCCGGGTAATATCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                                >  2:261339/1‑150 (MQ=255)
                                                                                     ccTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                                >  1:170821/1‑150 (MQ=255)
                                                                                      cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcg                                                               >  1:228475/1‑150 (MQ=255)
                                                                                      cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcg                                                               <  2:1241194/150‑1 (MQ=255)
                                                                                           aTGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCt                                                          <  1:824837/150‑1 (MQ=255)
                                                                                            tGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTg                                                         >  1:1166714/1‑150 (MQ=255)
                                                                                               cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAg                                                      >  2:964703/1‑150 (MQ=255)
                                                                                               cGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAg                                                      >  2:988768/1‑150 (MQ=255)
                                                                                                gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGc                                                     <  1:531011/150‑1 (MQ=255)
                                                                                                               gTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTaa                                     >  1:535909/1‑150 (MQ=255)
                                                                                                                    tatGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCt                                <  1:1013937/150‑1 (MQ=255)
                                                                                                                       gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTACCCGGCGGGAATGcttcag                             <  2:797416/150‑6 (MQ=255)
                                                                                                                         tCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCAcc                           <  2:497345/150‑1 (MQ=255)
                                                                                                                            gTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcgcg                        >  2:95969/1‑150 (MQ=255)
                                                                                                                              tGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCggg                      >  2:1023544/1‑150 (MQ=255)
                                                                                                                               gAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGc                     >  2:251026/1‑150 (MQ=255)
                                                                                                                                 aGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAg                   >  2:236347/1‑150 (MQ=255)
                                                                                                                                 aGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAg                   >  2:931174/1‑150 (MQ=255)
                                                                                                                                   cAATTATTCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 >  2:317041/1‑150 (MQ=255)
                                                                                                                                   acttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 >  2:86967/6‑150 (MQ=255)
                                                                                                                                    cttgataCTGTATGAGCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTgg                >  1:528576/5‑150 (MQ=255)
                                                                                                                                    aaTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTgg                >  1:293258/1‑150 (MQ=255)
                                                                                                                                     aTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGc               >  2:279287/1‑150 (MQ=255)
                                                                                                                                         taCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaaa           <  1:951521/150‑1 (MQ=255)
                                                                                                                                          aCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAg          >  1:949060/1‑150 (MQ=255)
                                                                                                                                           cTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGc         <  2:26177/150‑1 (MQ=255)
                                                                                                                                             gTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATg                                   >  2:76990/1‑122 (MQ=255)
                                                                                                                                             gTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATg                                   <  1:76990/122‑1 (MQ=255)
                                                                                                                                             gTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCaa       >  2:910573/1‑150 (MQ=255)
                                                                                                                                                tGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAAcc    <  1:988768/150‑1 (MQ=255)
                                                                                                                                                  cTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 >  2:775462/1‑135 (MQ=255)
                                                                                                                                                  cTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                 <  1:775462/135‑1 (MQ=255)
                                                                                                                                                  cTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAg  >  1:743569/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAG  >  NC_000913/2823692‑2823987

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