Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,800 A→G 7.1% intergenic (‑30/+56) 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,8000AG 7.1% 374.8 / 14.2 127intergenic (‑30/+56)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (66/52);  new base G (4/5);  total (70/57)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 127 total aligned reads displayed.

ATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCC  >  NC_000913/1315653‑1315948
                                                                                                                                                         |                                                                                                                                                    
aTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                    >  2:799134/1‑150 (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                                                                                                                                       >  1:451614/1‑150 (MQ=255)
             tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                       >  1:1430330/1‑150 (MQ=255)
             tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                       >  2:1155895/1‑150 (MQ=255)
             tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                       >  2:386884/1‑150 (MQ=255)
              gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                      >  1:1202128/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:351676/1‑150 (MQ=255)
                 gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATAATTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                                   <  1:1008847/150‑1 (MQ=255)
                   tctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTa                                                                                                                                 >  1:1232163/1‑150 (MQ=255)
                      tGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCGGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCt                                                                                                                              <  2:1232988/150‑1 (MQ=255)
                        ggAGTTATGCCAGTGCCAGGTTTAATTGTTTTTCTGCGCTGTATGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGAGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                            <  2:1057074/149‑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                                                                                                                                        >  2:1485586/1‑124 (MQ=255)
                                      gCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                                                                        <  1:1485586/124‑1 (MQ=255)
                                        ccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                                                                  >  1:308737/1‑150 (MQ=255)
                                        ccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                            <  1:707873/150‑1 (MQ=255)
                                           ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                         >  1:13897/1‑150 (MQ=255)
                                               aaTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                     >  1:168111/1‑150 (MQ=255)
                                                   gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATGTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                 <  2:1107492/150‑1 (MQ=255)
                                                   gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                 <  2:451614/150‑1 (MQ=255)
                                                   gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                 <  2:1195411/150‑1 (MQ=255)
                                                    tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGa                                                                                                >  1:581511/1‑150 (MQ=255)
                                                         cTGCGCTGTAGGTATCTGAAAGCTGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                           >  1:1153854/1‑150 (MQ=255)
                                                          tGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAgg                                                                                          >  1:1262565/1‑150 (MQ=255)
                                                           gcgcTGTAGGTATCTGAAAGCAGGTGAATATATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGa                                                                                         <  2:1598498/150‑1 (MQ=255)
                                                           gcgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGa                                                                                         >  1:671681/1‑150 (MQ=255)
                                                            cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                        <  1:1329209/150‑1 (MQ=255)
                                                             gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCATCAATGGTC‑‑‑‑‑‑TACATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGAAGATGCAAAAATCTTTATAAGCTGATATTATTAAATAACGGTGAATCATGAta                                                                                       >  2:1696133/1‑149 (MQ=255)
                                                              cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAttt                                                                                      >  1:997016/1‑150 (MQ=255)
                                                                 tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGcc                                                                                   <  2:1118779/150‑1 (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                                                                             >  2:115053/1‑150 (MQ=255)
                                                                       tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                             >  1:1567445/1‑150 (MQ=255)
                                                                        cTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATtta                                                                            >  2:335018/1‑150 (MQ=255)
                                                                        cTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TACATATTCATATTTTTCTCAAGTGAAAACACTCCGAGCCAAACGGACCGCAGAAGAAATTAACTATATCAGCTGATAATAATAATTAACGGTGAATCAGGATTTGCCTGATtta                                                                            >  2:1541520/1‑150 (MQ=255)
                                                                         tGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttat                                                                           >  1:619030/1‑150 (MQ=255)
                                                                          gAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                          >  1:969543/1‑150 (MQ=255)
                                                                          gAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                          >  2:334079/1‑150 (MQ=255)
                                                                           aaaGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                         <  2:1037660/150‑1 (MQ=255)
                                                                              gCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                      <  2:1430330/150‑1 (MQ=255)
                                                                               cAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccac                                                                     >  1:1586651/1‑150 (MQ=255)
                                                                                aGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccacc                                                                    <  2:449968/150‑1 (MQ=255)
                                                                                aGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccacc                                                                    >  2:773446/1‑150 (MQ=255)
                                                                                        aaataCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                            >  1:502022/1‑150 (MQ=255)
                                                                                           taCATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAc                                                         <  2:1202128/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)
                                                                                                         tC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                           <  1:973554/150‑1 (MQ=255)
                                                                                                                      attcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                                    <  2:1389768/150‑1 (MQ=255)
                                                                                                                          atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                >  2:773674/1‑150 (MQ=255)
                                                                                                                          atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                >  2:655418/1‑150 (MQ=255)
                                                                                                                            attTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                              <  1:433686/150‑1 (MQ=255)
                                                                                                                             tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                             <  2:1336931/150‑1 (MQ=255)
                                                                                                                                ttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                          <  1:1155895/150‑1 (MQ=255)
                                                                                                                                  ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCTCCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg                        >  1:337716/1‑150 (MQ=255)
                                                                                                                                   tcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCgg                       <  2:871741/150‑1 (MQ=255)
                                                                                                                                    ccAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                      >  2:589975/1‑150 (MQ=255)
                                                                                                                                     cAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                     >  2:1304910/1‑150 (MQ=255)
                                                                                                                                      tgTGAAATCACTGCGAGCCATACGGCCCGCAGATTCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTTTTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATc                    <  1:1595305/149‑1 (MQ=255)
                                                                                                                                      aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATc                    >  2:1528573/1‑150 (MQ=255)
                                                                                                                                           aaTCACTTCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtttt               <  2:1068050/150‑1 (MQ=255)
                                                                                                                                           aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtttt               >  1:21105/1‑150 (MQ=255)
                                                                                                                                            gcatcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa             >  2:845579/5‑150 (MQ=255)
                                                                                                                                            aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt              >  2:484404/1‑150 (MQ=255)
                                                                                                                                             catcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa             >  1:377763/4‑150 (MQ=255)
                                                                                                                                               aCTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaa           <  1:489164/150‑1 (MQ=255)
                                                                                                                                                  gCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt        >  1:806874/1‑150 (MQ=255)
                                                                                                                                                  gCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt        >  2:703937/1‑150 (MQ=255)
                                                                                                                                                      gCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg    <  2:964484/150‑1 (MQ=255)
                                                                                                                                                       ccATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc   <  1:1695178/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                    
ATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCC  >  NC_000913/1315653‑1315948

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: