Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,807 C→T 7.4% intergenic (‑37/+49) 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,8070CT 7.4% 375.7 / 11.3 124intergenic (‑37/+49)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base C (61/52);  new base T (4/5);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATG  >  NC_000913/1315659‑1315954
                                                                                                                                                          |                                                                                                                                                   
cATTTGATGTTGCTCTTGCGAGTTTTGACAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGccc                                                                                                                                                    >  2:57323/1‑150 (MQ=255)
 aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                   >  1:1012015/1‑150 (MQ=255)
 aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                   >  2:1645442/1‑150 (MQ=255)
   ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCa                                                                                                                                                 >  2:750193/1‑150 (MQ=255)
    tGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAg                                                                                                                                                >  1:1257425/1‑150 (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:1430330/1‑150 (MQ=255)
       tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                             >  1:451614/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‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                                         <  1:968504/150‑1 (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)
                  ggAGTTATGCCAGTGCCAGGTTTAATTGTTTTTCTGCGCTGTATGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGAGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                  <  2:1057074/149‑1 (MQ=255)
                     gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                                                                               <  2:1403552/150‑1 (MQ=255)
                          gCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCATCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGATg                                                                                                                          >  2:623771/1‑150 (MQ=255)
                            aaGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                                                                        <  2:1059465/149‑1 (MQ=255)
                                gCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                                                                              >  2:1485586/1‑124 (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)
                                     ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                               >  1:315619/1‑150 (MQ=255)
                                     ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                               >  2:1690834/1‑150 (MQ=255)
                                        tAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTaa                                                                                                            >  2:597405/1‑150 (MQ=255)
                                             gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATGTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                       <  2:1107492/150‑1 (MQ=255)
                                             gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                       <  2:1195411/150‑1 (MQ=255)
                                             gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                       <  2:451614/150‑1 (MQ=255)
                                              tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGa                                                                                                      <  1:1694908/150‑1 (MQ=255)
                                              tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGa                                                                                                      >  1:581511/1‑150 (MQ=255)
                                               ttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGaa                                                                                                     <  1:394133/150‑1 (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)
                                                       gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                                                                             <  1:1645442/150‑1 (MQ=255)
                                                        cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAttt                                                                                            >  1:997016/1‑150 (MQ=255)
                                                           tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGcc                                                                                         <  2:1118779/150‑1 (MQ=255)
                                                             ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                                       >  2:1520600/1‑150 (MQ=255)
                                                             ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                                       >  1:797407/1‑150 (MQ=255)
                                                             ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCACGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                                                                       >  1:168388/1‑150 (MQ=255)
                                                                aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                    <  1:1565247/150‑1 (MQ=255)
                                                                aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                    <  1:386884/150‑1 (MQ=255)
                                                                 tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                                   >  2:115053/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                                                                                >  2:334079/1‑150 (MQ=255)
                                                                    gAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                                >  1:969543/1‑150 (MQ=255)
                                                                     aaaGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                               <  2:1037660/150‑1 (MQ=255)
                                                                        gCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                            <  2:1430330/150‑1 (MQ=255)
                                                                        gCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                            <  1:767985/150‑1 (MQ=255)
                                                                         cAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccac                                                                           >  1:1586651/1‑150 (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)
                                                                                       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:123347/1‑150 (MQ=255)
                                                                                            tCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttt                                                        <  1:115053/150‑1 (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‑‑‑‑‑‑TTCATATTCATATTTTTCTACAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGGGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttttt                                                     <  1:57323/150‑1 (MQ=255)
                                                                                                   tC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                                 <  1:973554/150‑1 (MQ=255)
                                                                                                    c‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCg                                                <  1:669403/150‑1 (MQ=255)
                                                                                                                attcatattTTTATCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                                          <  2:724195/150‑1 (MQ=255)
                                                                                                                   catattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                                       >  1:1394993/1‑150 (MQ=255)
                                                                                                                    atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                      >  2:773674/1‑150 (MQ=255)
                                                                                                                      attTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                    <  1:433686/150‑1 (MQ=255)
                                                                                                                       tttttCTCCAGTGAAATCACTGCGAGCCATTCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                                   <  2:149048/150‑1 (MQ=255)
                                                                                                                       tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                                   <  2:1336931/150‑1 (MQ=255)
                                                                                                                          ttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                <  2:810208/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)
                                                                                                                                    ttgcatcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                      >  1:1554284/7‑150 (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:484404/1‑150 (MQ=255)
                                                                                                                                      aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt                    <  2:419685/150‑1 (MQ=255)
                                                                                                                                            gCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt              >  2:703937/1‑150 (MQ=255)
                                                                                                                                            gCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt              >  1:806874/1‑150 (MQ=255)
                                                                                                                                                 ccATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc         <  1:1695178/150‑1 (MQ=255)
                                                                                                                                                  cATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgcc        >  2:1388874/1‑150 (MQ=255)
                                                                                                                                                     aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCt     <  2:1567445/150‑1 (MQ=255)
                                                                                                                                                        gCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg  <  2:1152523/150‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                   
CATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTC‑‑‑‑‑‑TTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATG  >  NC_000913/1315659‑1315954

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