Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,800 A→G 5.2% 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 5.2% 239.1 / 3.2 77intergenic (‑30/+56)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

TTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCC  >  NC_000913/1315661‑1315948
                                                                                                                                           |                                                                                                                                                    
tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                             >  7:280921/1‑149 (MQ=255)
tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                             >  5:38803/1‑149 (MQ=255)
     tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                                                                        <  3:82340/149‑1 (MQ=255)
      gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                                                                        <  4:2233/148‑1 (MQ=255)
         gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                    >  1:59369/1‑149 (MQ=255)
          ctctTGCGAGTTTTGCCCGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCCGAACGCATGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAACTCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                                   <  6:110264/149‑1 (MQ=255)
               gCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATAGATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCTGATGCAATTATc                                                                                                                               <  8:73329/148‑1 (MQ=255)
                   gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGTATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                          <  8:58459/149‑1 (MQ=255)
                   gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAtt                                                                                                                                  <  1:186534/141‑1 (MQ=255)
                   gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAtt                                                                                                                                  >  2:186534/1‑141 (MQ=255)
                   gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                          <  7:62923/149‑1 (MQ=255)
                    ttttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                            >  8:105419/1‑146 (MQ=255)
                    ttttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                            <  7:105419/146‑1 (MQ=255)
                     gtttccAGTGCCCTGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGATTAAATACATCTTCAATGGTCTTCATATTCTTATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                                                                         <  2:109470/144‑1 (MQ=255)
                        gCCAGTGCCCGGGTTAATTGTTTTTCGGCGCTGTAGTTATCGGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGAAATTATCTTAATCAGCt                                                                                                                     >  7:327851/1‑149 (MQ=255)
                                   ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                          >  4:78844/1‑149 (MQ=255)
                                   ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                          <  6:38803/149‑1 (MQ=255)
                                     ttAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCAAGTGAAATCACTTCGATCCATACGTACCGCAGATGCAATTATATTTATCAGCTGATATTaataat                                                                                                         >  5:117642/1‑148 (MQ=255)
                                         ttGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACg                                                                                                    >  5:207417/1‑149 (MQ=255)
                                          tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCAAGTGAAATCACTGCTAGCCATACGTCCCTAAGATGAAATTATCTTTATCATATGATATTAATAATTAActg                                                                                                   >  7:349848/1‑147 (MQ=255)
                                            tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                  <  7:13469/148‑1 (MQ=255)
                                            tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                  <  6:167154/148‑1 (MQ=255)
                                                  tGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                           >  6:355229/1‑149 (MQ=255)
                                                    cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCa                                                                                            <  3:55962/146‑1 (MQ=255)
                                                              aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTAATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGAAATTATCTTTATCAGATGATATTAATAATTAACGGTGAATCAGGAttttcattat                                                                               >  7:336663/1‑142 (MQ=255)
                                                               tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGaa                                                                                               >  2:303912/1‑132 (MQ=255)
                                                               tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGaa                                                                                               <  1:303912/132‑1 (MQ=255)
                                                                   aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                          <  5:165916/149‑1 (MQ=255)
                                                                     aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTAGCCTGATTTATTac                                                                        >  2:168928/1‑149 (MQ=255)
                                                                      gCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                                                                        >  3:111524/1‑148 (MQ=255)
                                                                           tGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttat                                                                           <  2:153123/140‑1 (MQ=255)
                                                                           tGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttat                                                                           >  1:153123/1‑140 (MQ=255)
                                                                            gAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATATGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCt                                                                 >  7:23562/1‑149 (MQ=255)
                                                                             aatatataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTCAATCAGGATTTGCCTGATTTATTACCACCGCtt                                                                >  7:14853/1‑149 (MQ=255)
                                                                              ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGGGAAATCACTGGGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAAGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                               <  5:32004/149‑1 (MQ=255)
                                                                               taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                                                              >  3:197778/1‑149 (MQ=255)
                                                                               taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                                                              <  5:355229/149‑1 (MQ=255)
                                                                                 aataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTTAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                            <  3:258613/149‑1 (MQ=255)
                                                                                  ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                            >  6:179269/1‑148 (MQ=255)
                                                                                            aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                 >  7:22416/1‑149 (MQ=255)
                                                                                            aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                 >  5:88254/1‑149 (MQ=255)
                                                                                              tGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttttt                                               <  2:59369/149‑1 (MQ=255)
                                                                                                  cTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAttt                                                                                      <  1:227833/106‑1 (MQ=255)
                                                                                                  cTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAttt                                                                                      >  2:227833/1‑106 (MQ=255)
                                                                                                    taatattcatattTTTCTGCATTGAAATCACTGCGAGCCATACGGCACGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCg                                          <  8:336663/146‑1 (MQ=255)
                                                                                                     catattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGc                                         <  8:280921/148‑1 (MQ=255)
                                                                                                          tcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                       >  1:325175/1‑145 (MQ=255)
                                                                                                           catattTTTCTCCAGTGAAATCTATGCTAGCCATACGGGCGGCAGATGCATTTATCTTTATCTTATGATATTATTAATTAACTGTGTATCAGGTTTTGCCTGATTTATTTCCACCGCTTTGTTGAGCGGCTTTTTTTGCCGCTTCAGa                                   <  8:361233/148‑1 (MQ=255)
                                                                                                           catattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                   <  6:207417/148‑1 (MQ=255)
                                                                                                           catattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACAACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGAt                                  >  4:263813/1‑149 (MQ=255)
                                                                                                            atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                         <  1:284245/77‑1 (MQ=255)
                                                                                                            atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                                                                         >  2:284245/1‑77 (MQ=255)
                                                                                                            atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                                 >  3:113917/1‑149 (MQ=255)
                                                                                                            atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                                 >  3:218952/1‑149 (MQ=255)
                                                                                                             tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                >  2:78654/1‑149 (MQ=255)
                                                                                                             tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCAGCTTCAGATgc                                >  6:44009/1‑149 (MQ=255)
                                                                                                              attTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcg                               >  1:348826/1‑149 (MQ=255)
                                                                                                               tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                              <  4:218952/149‑1 (MQ=255)
                                                                                                               tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                              <  1:168928/149‑1 (MQ=255)
                                                                                                                   cgttctgTGAAATCTCTGCGTGCAATACTGCCCTCAGATTGATTTTTCTTTATCAGCTGATATTAATTATTATCGGTGTATGTGGTTTTGGCTGATTTATTTCCACTGGTTTTTTGAGCGCGTTTTTTACGTGCTTCTGATTCGCTTTg                          <  8:218535/143‑1 (MQ=255)
                                                                                                                       cAGTTAAATCACTGCGAGCCATACGTCCCGCAGATGCAATTATCTTTTTCAGCTTATATTAATAATTAACGGTGAATAAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                         <  6:88254/146‑1 (MQ=255)
                                                                                                                       cAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                      <  3:263813/149‑1 (MQ=255)
                                                                                                                       cAGTGAAATCACTGCGAGACATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                      <  6:117642/149‑1 (MQ=255)
                                                                                                                        aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                     <  8:22416/149‑1 (MQ=255)
                                                                                                                          tGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                          >  7:218535/1‑142 (MQ=255)
                                                                                                                            aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                 >  3:50/1‑149 (MQ=255)
                                                                                                                             tgcatcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                >  2:244013/6‑149 (MQ=255)
                                                                                                                               catcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt              >  2:349966/4‑148 (MQ=255)
                                                                                                                                cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                           <  7:286609/71‑1 (MQ=255)
                                                                                                                                cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                           >  8:286609/1‑71 (MQ=255)
                                                                                                                                 aCTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACACGCTTCAGATGCGCGTTGCGGATCATTTTTaa            <  3:78844/149‑1 (MQ=255)
                                                                                                                                     cgggCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa          <  1:244013/144‑1 (MQ=255)
                                                                                                                                     cgggCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaa           <  1:349966/142‑1 (MQ=255)
                                                                                                                                     cGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt         >  4:123033/1‑148 (MQ=255)
                                                                                                                                       aGCCATAAGGCCCTCATATGCAATTATCTTTATCATCTGATAATAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc      <  4:197778/149‑1 (MQ=255)
                                                                                                                                        gCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc      >  2:214577/1‑148 (MQ=255)
                                                                                                                                        gCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc      <  1:214577/148‑1 (MQ=255)
                                                                                                                                           aTACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgcc  >  3:162515/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
TTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCC  >  NC_000913/1315661‑1315948

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