Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,796 A→G 7.1% intergenic (‑26/+60) yciF ← / ← yciG putative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,7960AG 7.1% 350.8 / 12.9 126intergenic (‑26/+60)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (66/51);  new base G (4/5);  total (70/56)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

Only 100 of 126 total aligned reads displayed.

ACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT  >  NC_000913/1315648‑1315942
                                                                                                                                                          |                                                                                                                                                  
 cTTAATTTTTCAGTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGTTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAATCAGGTGAATAAAAACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTACAGTGAATTAACTTCGAGc                                                                                                                                                  >  2:1504974/1‑150 (MQ=255)
        tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACgg                                                                                                                                           >  2:1204087/1‑150 (MQ=255)
        tttCATTTGATGTTGCTCTTGCGAGTGTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACgg                                                                                                                                           <  1:1296088/150‑1 (MQ=255)
           cATTTGATGTTGCTCTTGCGAGTTTTGACAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGccc                                                                                                                                        >  2:57323/1‑150 (MQ=255)
            aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                       >  2:1645442/1‑150 (MQ=255)
            aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                       >  1:1012015/1‑150 (MQ=255)
             tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                                                                      >  1:871741/1‑150 (MQ=255)
              ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCa                                                                                                                                     >  2:750193/1‑150 (MQ=255)
               tGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAg                                                                                                                                    >  1:1257425/1‑150 (MQ=255)
               tGATGTAGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TCCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAg                                                                                                                                    <  1:275895/150‑1 (MQ=255)
                gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                   <  2:538012/150‑2 (MQ=255)
                 aTGTTGCTCTTGCGAGTTTTGCCAGTTCCCGTTTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTTAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                                                                  <  1:766475/150‑3 (MQ=255)
                 aTGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                                                                  <  2:1186439/150‑3 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATATATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                 >  1:1498180/1‑150 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtga                                                                                                                                 <  1:1668966/150‑4 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                 >  2:386884/1‑150 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                 >  1:451614/1‑150 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                 >  2:1155895/1‑150 (MQ=255)
                  tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                 >  1:1430330/1‑150 (MQ=255)
                    ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCa                                                                                                                               >  1:810208/1‑150 (MQ=255)
                    ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCa                                                                                                                               <  2:696351/150‑1 (MQ=255)
                      gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgatttc                                                                                                                             <  1:750193/150‑8 (MQ=255)
                      gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                             <  1:968504/150‑1 (MQ=255)
                      gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                             >  1:351676/1‑150 (MQ=255)
                      gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATAATTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                             <  1:1008847/150‑1 (MQ=255)
                        tctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTa                                                                                                                           >  1:1232163/1‑150 (MQ=255)
                          ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCACTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATc                                                                                                                         <  2:1550901/150‑1 (MQ=255)
                           tGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCGGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCt                                                                                                                        <  2:1232988/150‑1 (MQ=255)
                                gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                                                                   <  2:1403552/150‑1 (MQ=255)
                                     gCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCATCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGATg                                                                                                              >  2:623771/1‑150 (MQ=255)
                                       cAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAATTACAGCTTCAAAGGTC‑‑‑‑‑‑TTCATATACATATTTTTCTCCAGTGAAATCACTCCGAGCCATACGGCCAGCTGATCCAATAAGCTAAAACAGCTGat                                                                                                            >  2:1560833/1‑150 (MQ=255)
                                       aaGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                                                            <  2:1059465/149‑1 (MQ=255)
                                           gCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                                                                  <  1:1485586/124‑1 (MQ=255)
                                             ccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                                                            >  1:308737/1‑150 (MQ=255)
                                             ccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                      <  1:1140437/150‑1 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                   >  2:1690834/1‑150 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                   >  1:13897/1‑150 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                   >  1:315619/1‑150 (MQ=255)
                                                   tAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTaa                                                                                                >  2:597405/1‑150 (MQ=255)
                                                    aaTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                               >  1:168111/1‑150 (MQ=255)
                                                        ttttttCTGCGCTGTAGGTATCTTAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATTTTTATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                           <  2:1092049/149‑1 (MQ=255)
                                                        gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATGTTTATCAGCTGATATTAATAATTAACGGTg                                                                                           <  2:1107492/150‑1 (MQ=255)
                                                        gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                           <  2:451614/150‑1 (MQ=255)
                                                          ttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGaa                                                                                         <  1:394133/150‑1 (MQ=255)
                                                              cTGCGCTGTAGGTATCTGAAAGCTGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                     >  1:1153854/1‑150 (MQ=255)
                                                                gcgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGa                                                                                   >  1:671681/1‑150 (MQ=255)
                                                                 cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                  <  1:1329209/150‑1 (MQ=255)
                                                                  gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                                                                 <  2:1012015/150‑1 (MQ=255)
                                                                  gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCATCAATGGTC‑‑‑‑‑‑TACATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGAAGATGCAAAAATCTTTATAAGCTGATATTATTAAATAACGGTGAATCATGAta                                                                                 >  2:1696133/1‑149 (MQ=255)
                                                                   cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAttt                                                                                >  1:997016/1‑150 (MQ=255)
                                                                        ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTTCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                           >  2:1513979/1‑150 (MQ=255)
                                                                        ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                           >  1:797407/1‑150 (MQ=255)
                                                                        ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                           >  2:1520600/1‑150 (MQ=255)
                                                                        ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCACGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                           >  1:168388/1‑150 (MQ=255)
                                                                           aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                        <  1:386884/150‑1 (MQ=255)
                                                                           aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                        <  1:1565247/150‑1 (MQ=255)
                                                                            tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                       >  1:1567445/1‑150 (MQ=255)
                                                                            tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                       >  2:115053/1‑150 (MQ=255)
                                                                             cTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TACATATTCATATTTTTCTCAAGTGAAAACACTCCGAGCCAAACGGACCGCAGAAGAAATTAACTATATCAGCTGATAATAATAATTAACGGTGAATCAGGATTTGCCTGATtta                                                                      >  2:1541520/1‑150 (MQ=255)
                                                                              tGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttat                                                                     >  1:619030/1‑150 (MQ=255)
                                                                               gAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                    >  2:334079/1‑150 (MQ=255)
                                                                                aaaGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                   <  2:1037660/150‑1 (MQ=255)
                                                                                   gCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                <  1:767985/150‑1 (MQ=255)
                                                                                   gCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                <  2:1430330/150‑1 (MQ=255)
                                                                                    cAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccac                                                               >  1:1586651/1‑150 (MQ=255)
                                                                                     aGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccacc                                                              >  2:773446/1‑150 (MQ=255)
                                                                                     aGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccacc                                                              <  2:449968/150‑1 (MQ=255)
                                                                                             aaataCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                      >  1:502022/1‑150 (MQ=255)
                                                                                             aaataCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                      >  2:117538/1‑150 (MQ=255)
                                                                                                taCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAc                                                   <  2:1202128/150‑1 (MQ=255)
                                                                                                 aCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAcc                                                  <  1:1204087/150‑1 (MQ=255)
                                                                                                  cATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg                                                 <  2:1257425/150‑1 (MQ=255)
                                                                                                    tCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGcc                                               >  2:489164/1‑150 (MQ=255)
                                                                                                     cTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt                                              >  1:1493131/1‑150 (MQ=255)
                                                                                                       tCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttt                                            <  1:115053/150‑1 (MQ=255)
                                                                                                       tCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttt                                            >  1:123347/1‑150 (MQ=255)
                                                                                                        cAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                           >  1:1367236/1‑150 (MQ=255)
                                                                                                         aaTGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttttt                                          >  1:1009666/1‑150 (MQ=255)
                                                                                                         aaTGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACCGCCCGCAGATGCAATTATCTTTATCAGAAGATATTAATAATTAACGGAGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCatttt                                          >  1:1003154/1‑150 (MQ=255)
                                                                                                          aTGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTGCAGTGAAATCACTGCGAGCCATACGGCGCGCATATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGTATCAGTATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttttt                                         <  1:623771/150‑1 (MQ=255)
                                                                                                          aTGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTACAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGGGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttttt                                         <  1:57323/150‑1 (MQ=255)
                                                                                                              tC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                     <  1:973554/150‑1 (MQ=255)
                                                                                                               c‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCg                                    <  1:669403/150‑1 (MQ=255)
                                                                                                                           attcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                              <  2:1389768/150‑1 (MQ=255)
                                                                                                                           attcatattTTTATCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                              <  2:724195/150‑1 (MQ=255)
                                                                                                                               atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                          >  2:773674/1‑150 (MQ=255)
                                                                                                                               atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                          >  2:655418/1‑150 (MQ=255)
                                                                                                                                  tttttCTCCAGTGAAATCACTGCGAGCCATTCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                       <  2:149048/150‑1 (MQ=255)
                                                                                                                                  tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                       <  2:1336931/150‑1 (MQ=255)
                                                                                                                                     ttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                    <  1:1155895/150‑1 (MQ=255)
                                                                                                                                     ttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                    <  2:810208/150‑1 (MQ=255)
                                                                                                                                        tcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCgg                 <  2:871741/150‑1 (MQ=255)
                                                                                                                                           tgTGAAATCACTGCGAGCCATACGGCCCGCAGATTCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTTTTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATc              <  1:1595305/149‑1 (MQ=255)
                                                                                                                                                aaTCACTTCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtttt         <  2:1068050/150‑1 (MQ=255)
                                                                                                                                                aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtttt         >  1:21105/1‑150 (MQ=255)
                                                                                                                                                 gcatcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa       >  2:845579/5‑150 (MQ=255)
                                                                                                                                                 ataaCTGCGAGCAATACGGCCCGAAGATGCAATTATCTTTATAAGCTGATATTAATAATTAACGGTGAATCAGGATTTGACTGATTTATTACAACCGCTTTGTTGACCTCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt        >  2:490872/4‑150 (MQ=255)
                                                                                                                                                 aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt        <  2:419685/150‑1 (MQ=255)
                                                                                                                                                  catcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa       >  1:377763/4‑150 (MQ=255)
                                                                                                                                                    aCTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaa     <  1:489164/150‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                  
ACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT  >  NC_000913/1315648‑1315942

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