Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 6.7% intergenic (+57/+27) chiQ → / ← fur chitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,1730AG 6.7% 161.5 / 2.7 61intergenic (+57/+27)chiQ/furchitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base G (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGC  >  NC_000913/710027‑710314
                                                                                                                                                  |                                                                                                                                             
tgtCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAgg                                                                                                                                             <  1:233127/149‑1 (MQ=255)
 gtCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGt                                                                                                                                            <  1:837341/149‑1 (MQ=255)
     gCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                        >  2:1024039/1‑149 (MQ=255)
          cAACGCAAACCGGAAATGATCACTCGGTGAAAGCCGATTGTGATAACGTCAGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttt                                                                                                                                   <  1:195302/149‑2 (MQ=255)
             cGCTAAGCGTAAATGAGCACGCGGCTAAAGCCGAAGGTGATAAAGCCTGGCATGAAATACGCAGTAATAACAAATATGTGAGAGCTGTAACTCGCGCTTTTCTTATTTCCCTTGCATTAAAAAGCCAACCTGCGGGTTGGCttttttat                                                                                                                                <  2:600270/149‑5 (MQ=255)
             cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCg                                                                                                                                >  1:820047/1‑149 (MQ=255)
                aaaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGGCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTGTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGGTc                                                                                                                             >  1:644193/1‑149 (MQ=255)
                   ccGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACTCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAgg                                                                                                                          <  1:218817/149‑1 (MQ=255)
                      gAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctg                                                                                                                       <  2:332070/149‑1 (MQ=255)
                      gAAATGATCAGGCGGAGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTATCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctg                                                                                                                       <  1:674533/149‑1 (MQ=255)
                       aaaTGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctgg                                                                                                                      <  1:499922/149‑1 (MQ=255)
                            aTCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAt                                                                                                                 >  2:943400/1‑149 (MQ=255)
                            aTCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAt                                                                                                                 <  1:562732/149‑1 (MQ=255)
                              cAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAttt                                                                                                               <  1:3982/149‑1 (MQ=255)
                              cAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCGGTAACTCTCGCGTTTCTTCTTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGCCTGGCTTAttt                                                                                                               <  2:819675/149‑1 (MQ=255)
                                ggCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTCGCGTTTCTCGTACAGGATGGCTTATTTGc                                                                                                             >  2:189451/1‑149 (MQ=255)
                                    gTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTc                                                                                                         >  2:214237/1‑149 (MQ=255)
                                    gTGAAAGCCGATTGTGATAAAGTCTGGAATGAAAAACGAACGAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATGTCCCTTGCATAAACAAGCCAACCCGCAGGTTGGCGTTTCTCGATCAGGCGGGCTTATTTGCCTTc                                                                                                         <  2:673533/149‑1 (MQ=255)
                                               ttGTGATAAAGTCTGGCAGGAAATACGCAGTTATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGtt                                                                                              <  2:93589/149‑1 (MQ=255)
                                               tggtgATAAATTCTTGCAGGAAATACTCAGTACTAAAAAATAAGGGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCAAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGtt                                                                                              <  2:900204/147‑1 (MQ=255)
                                                tgtgATAAAGTCTGTCAGGAAATACGCAGTAATAACAAATTAGTGAGAGCTGTAACTGTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCTACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTc                                                                                             <  1:703167/149‑1 (MQ=255)
                                                   gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAta                                                                                          >  1:488311/1‑148 (MQ=255)
                                                   gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGATTATTTGCCTTCGTGCGCATGTACATc                                                                                          >  2:704953/1‑149 (MQ=255)
                                                          tCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                                                                   <  2:754009/149‑1 (MQ=255)
                                                              gCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAAt                                                                               >  1:636016/1‑149 (MQ=255)
                                                                aGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCg                                                                             >  1:970643/1‑149 (MQ=255)
                                                                     aTACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCTTGCGCATTTTCATCTTCGCGGCAATCGCCTTc                                                                        >  1:116634/1‑149 (MQ=255)
                                                                      tACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTACg                                                                       >  1:998820/1‑149 (MQ=255)
                                                                       aagCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGAATTTCTTATTACCCTAGCATAAAAACGCCAACCCGCAGGTTGGCTTTTCTCTATCATGCTGGCTTAATAGACTTCGTGAGCATGTTCATCTTCGCTGCACTCTCCtgagg                                                                      >  1:559922/3‑145 (MQ=255)
                                                                        cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGc                                                                     <  1:596476/149‑1 (MQ=255)
                                                                        cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGCTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGc                                                                     >  1:967546/1‑149 (MQ=255)
                                                                         gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTGCCTTGCATAAAAAACCCAACCCGCAGGTTGGCTTTTATCGTTCAGGCTGGCTTATTTGCCTACGTCCGCAAGTTCATCTTCGCGGCAATCGCCTTCGGca                                                                    >  2:872470/1‑149 (MQ=255)
                                                                         gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcc                                                                    >  2:613958/1‑148 (MQ=255)
                                                                         gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTATTCTTGTTTCGCTTGCATAAAAAATCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGACGTCGTGCGCATCTTCATCTTCGCGGCAATCGCCTTCGGca                                                                    >  2:830399/1‑149 (MQ=255)
                                                                                   cAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGt                                                          >  1:895221/1‑149 (MQ=255)
                                                                                   cAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTAATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGt                                                          >  1:477711/1‑149 (MQ=255)
                                                                                    aaaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTa                                                         >  2:847606/1‑149 (MQ=255)
                                                                                    aaaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTa                                                         >  1:490678/1‑149 (MQ=255)
                                                                                    aaaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTAGCATAAAAAACCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTGGTGCGCATGTGCATCTTCGCGGCAATCGCCTTCGCCACAGTGACCGTa                                                         >  2:66459/1‑149 (MQ=255)
                                                                                     aaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTaa                                                        <  1:898698/149‑1 (MQ=255)
                                                                                     aaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTaa                                                        <  2:216650/149‑1 (MQ=255)
                                                                                       tAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAg                                                      <  2:146059/149‑1 (MQ=255)
                                                                                        aaGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGa                                                     >  1:149501/1‑149 (MQ=255)
                                                                                           tGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATag                                                  <  2:493387/149‑1 (MQ=255)
                                                                                           tGAGAGATGTAACTCTCGCTTTTCTTATTTACCTTGCATAAAAACGCCAACCCGCAGGTTTGCTTTTCTCGTTCAGGCGGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCACTCGCCTTCGGCACAGTGACCGTAAAGATag                                                  <  2:342406/149‑1 (MQ=255)
                                                                                               agCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGAc                                              <  2:652486/149‑1 (MQ=255)
                                                                                                 cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGAGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtg                                            >  1:783618/1‑149 (MQ=255)
                                                                                                     acctctCGCTTTTCTTATTTCGCTTGCATTAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGt                                        <  2:164341/147‑1 (MQ=255)
                                                                                                      aCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGtt                                       <  1:449351/149‑1 (MQ=255)
                                                                                                                    tATTTCCCTTGCATACAATAGCCAACACGCTGGTTGGCTTTTCTAGTTAAGGCTTGCTTATTTGCATTCTTGCGGATGTTCATCTTCGCGGCAATCGCATTCGGCACAGTGTCCGTAATGTTATAGACTGTGGTTAGTCAGGCGAATGc                         <  2:40231/149‑1 (MQ=255)
                                                                                                                    tATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGc                         <  1:547293/149‑1 (MQ=255)
                                                                                                                     aTTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTTAAGATAGAGACTGTGGTTAGTCAGGCGAATGcc                        >  1:305935/1‑149 (MQ=255)
                                                                                                                           cTTGCATAAAAAAGCCAACCCGCAGGTTGGATTTTCTCGTTCAGGCTTGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGTCACAGTGACCGTAAAGATAGAGACGGTGGTTAGTCAGGCGAATGCCATGttt                  >  2:17553/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCgg             >  1:298574/1‑149 (MQ=255)
                                                                                                                                   aaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACGGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCaa          >  2:775701/1‑149 (MQ=255)
                                                                                                                                    aaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAt         >  1:376690/1‑149 (MQ=255)
                                                                                                                                      aaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCATTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAttt       >  1:672135/1‑149 (MQ=255)
                                                                                                                                      aaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAttt       >  1:277347/1‑149 (MQ=255)
                                                                                                                                        gCCAACCCGCAGTGTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa     <  2:820047/149‑1 (MQ=255)
                                                                                                                                        cccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa     <  1:850332/148‑1 (MQ=255)
                                                                                                                                           aaCCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGc  >  1:653549/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGC  >  NC_000913/710027‑710314

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: ‑