Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,674,690 G→T 8.8% intergenic (+316/‑95) yhjD → / → yhjE inner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,674,6900GT 8.8% 215.7 / 2.9 78intergenic (+316/‑95)yhjD/yhjEinner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base T (6/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCCAACTCATAACCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATACACGCCGATCAAC  >  NC_000913/3674549‑3674838
                                                                                                                                             |                                                                                                                                                    
tCCAACTCATAACCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTg                                                                                                                                               <  1:215066/149‑1 (MQ=255)
      tCATAACCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGCTTTTTGACCAtt                                                                                                                                         >  2:182313/1‑149 (MQ=255)
        aTAACCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTTACCATTTc                                                                                                                                       >  2:186879/1‑149 (MQ=255)
          aaCCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGTATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCATTTTTGCCTGATTTTTGACCATTTCcg                                                                                                                                     <  1:497229/149‑1 (MQ=255)
          aaCCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCcg                                                                                                                                     >  2:199748/1‑149 (MQ=255)
            ccccAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCAAAAAATTATCCCCTTTTTGCCTGATTTTTGCCCATTTCcgcg                                                                                                                                   >  2:13764/1‑149 (MQ=255)
             cccAGCATATAAATCCAGTTGGTAACTTGTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGa                                                                                                                                  <  2:428931/149‑1 (MQ=255)
                  cagataAATCCAGTTTGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCCTTTCATAGAAGAAAAATCACTGGCCTAAACATTCTCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGt                                                                                                                             <  1:687754/146‑1 (MQ=255)
                  cATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGt                                                                                                                             <  2:407783/149‑1 (MQ=255)
                    tcaaaaGCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTa                                                                                                                           <  2:657788/146‑1 (MQ=255)
                    tataCATCCAGTTCGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTa                                                                                                                           <  2:148526/149‑1 (MQ=255)
                    tataAATCCCGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTa                                                                                                                           <  2:161929/149‑1 (MQ=255)
                    tataAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTa                                                                                                                           <  2:471600/149‑1 (MQ=255)
                    tataAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTGCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTa                                                                                                                           <  2:417521/149‑1 (MQ=255)
                      ttaaTCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTCCTTTATTGTTGACCATTTCCGCGATTTGTTaca                                                                                                                         <  2:151291/147‑1 (MQ=255)
                      taAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTaca                                                                                                                         <  1:629202/149‑1 (MQ=255)
                       aaaTCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTacac                                                                                                                        >  2:668653/1‑149 (MQ=255)
                        aagccAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTAcaca                                                                                                                       <  1:268328/146‑1 (MQ=255)
                         aTCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAACATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGACTTGTTACACAt                                                                                                                      >  2:27843/1‑149 (MQ=255)
                              gTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGaaa                                                                                                                 <  1:274612/149‑1 (MQ=255)
                                 ggTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAatat                                                                                                              >  2:570502/1‑149 (MQ=255)
                                  gttatttttATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATc                                                                                                             <  2:107565/144‑1 (MQ=255)
                                   tAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGCCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCa                                                                                                            <  1:598334/149‑1 (MQ=255)
                                     aCTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAatattcat                                                                                                          >  2:116812/1‑145 (MQ=255)
                                     aCTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATACCACt                                                                                                          <  1:454056/149‑1 (MQ=255)
                                           atttaACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTAATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCt                                                                                                    <  2:264602/149‑1 (MQ=255)
                                            tttaACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGcga                                                                                                   >  1:215031/1‑147 (MQ=255)
                                            tttaACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCATAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGcgg                                                                                                   >  1:105111/1‑147 (MQ=255)
                                             ttaACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCtgt                                                                                                  <  1:112634/149‑1 (MQ=255)
                                               aaCCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGc                                                                                                <  1:138716/149‑1 (MQ=255)
                                                  cTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTCCTGTGCGTa                                                                                             >  1:53056/1‑149 (MQ=255)
                                                     aaaCCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAata                                                                                          <  2:207408/149‑1 (MQ=255)
                                                        ccAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTTTTTTTTTTTTAACCATTACCGCGCTTTCGTATACATTGAAAGTTCACTTTGGGTGTGCGTAGTATgg                                                                                       >  2:309454/1‑149 (MQ=255)
                                                         cAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGACTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGc                                                                                      <  2:79477/149‑1 (MQ=255)
                                                          aGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCt                                                                                     <  1:12097/149‑1 (MQ=255)
                                                            ttttATGCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTAt                                                                                   >  2:181915/1‑149 (MQ=255)
                                                            ttttATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTAt                                                                                   >  2:518171/1‑149 (MQ=255)
                                                                    aCTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTTTTTTTTTTTTGACCATTTACGCGATTTGTGAGACTTTGACATATCACTTTTGCTTAGAGTATTATGACTATTCGttctc                                                                           >  2:472869/1‑146 (MQ=255)
                                                                    aCTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGc                                                                           <  1:207702/149‑1 (MQ=255)
                                                                         ttATAAAATTATGTGAAGCATTTCATAGAAGAAACCTCACGGGCCTAAACATTATCCCCTGTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCaaaa                                                                      >  2:491088/1‑149 (MQ=255)
                                                                          tataCAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCaaaaa                                                                     >  2:119791/1‑149 (MQ=255)
                                                                          tataAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGGCaaaaa                                                                     >  1:341700/1‑149 (MQ=255)
                                                                           ataAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGcacagaaa                                                                    >  1:672914/1‑142 (MQ=255)
                                                                            taAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTTCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAAt                                                                   >  2:100581/1‑149 (MQ=255)
                                                                              aaaTTATGTGCAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACGATTTCCGCGATCAGCGACACACCGTAATAGCACTTTTGCTGTGCGTAATATGGCTAATCGTTAGCCAAAACATaa                                                                 >  1:200869/1‑149 (MQ=255)
                                                                                 ttATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACTCATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGaa                                                              <  1:686527/149‑1 (MQ=255)
                                                                                  tATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGaaa                                                             <  1:570683/149‑1 (MQ=255)
                                                                                      tgAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTCTACGGTAGGCAAAAAATAAGAAAAGag                                                         >  1:499471/1‑148 (MQ=255)
                                                                                       gAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTTTTTTTTTTTTTACCATTTCAGCGCTTTGTTACACTTTGAAATATCACTTTTGCTGTGCGTAATATTGCTCTTCGTTAGCCAAAAAATAAGAAAAGAtt                                                        >  2:342344/1‑149 (MQ=255)
                                                                                        aaGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTTTTTTTTTTTTCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCCATTCGTTACCCACACAATAAGAAAAGATTa                                                       >  2:293103/1‑149 (MQ=255)
                                                                                         aGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTGTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTAt                                                      <  2:241600/149‑1 (MQ=255)
                                                                                              ttCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATAGCTATTCGTTAGCCAAAAAATAAAAAAAGATTATgcaag                                                 >  1:574869/1‑149 (MQ=255)
                                                                                                 aTAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGCAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATgcagccaa                                              >  1:211116/1‑144 (MQ=255)
                                                                                                  tAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAAc                                             <  2:215031/149‑1 (MQ=255)
                                                                                                  tAGAAGAAAAATCACTGGCCTAAACATTACCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAACTAAGAAAAGATTATGCAAGCAAc                                             <  1:372690/149‑1 (MQ=255)
                                                                                                    gaagaaAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATCAGAAAAGATTATGCAAGCAACAg                                           >  1:181525/1‑149 (MQ=255)
                                                                                                     aagaaAAATCACTGGCCTAAACATGATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGc                                          <  1:46172/149‑1 (MQ=255)
                                                                                                      agaaAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGcc                                         <  1:668653/149‑1 (MQ=255)
                                                                                                      agaaAAAGCACTGGCCTAAACATTCTCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCGGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGcc                                         <  1:677044/149‑1 (MQ=255)
                                                                                                       gaaAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTTACCATTTTCGCGACTTGTTACACCTTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCCAAAAATAAGAAAAGATTATGCAAGCTACAGCca                                        >  2:347256/1‑149 (MQ=255)
                                                                                                       gaaAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATCGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCca                                        >  1:547795/1‑149 (MQ=255)
                                                                                                          aaaTCACTGGCCTAAACATTATCCCCTTTTTTTTTTTTTTTTGACCATTACGGCGATTTGTTACACATTGAAATATCCCTTTTGCTGTGCGTAATCTTGCTATTCGTTATCCAAAAAATCAGACACGATTATGCCAGCTACGGCCtcaa                                     >  2:353254/1‑149 (MQ=255)
                                                                                                          aaaTCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCacaa                                     <  1:181915/149‑1 (MQ=255)
                                                                                                           aaTCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGCCTATTCGTCATCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCacaac                                    <  1:614678/149‑1 (MQ=255)
                                                                                                                 tGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTCTGCAAGCAACAGCCACAACACTCGa                              <  1:413602/149‑1 (MQ=255)
                                                                                                                 tGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACAATCGa                              >  1:270499/1‑149 (MQ=255)
                                                                                                                  ggCCTAAACATTATCCCCTTTTTTTTTTTTTTTTTACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCATCACCCACAACACTCCAc                             >  2:260926/1‑149 (MQ=255)
                                                                                                                  ggCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGAc                             <  1:199748/149‑1 (MQ=255)
                                                                                                                   gCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACTATTTAAGCTATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGAcc                            >  1:337865/1‑149 (MQ=255)
                                                                                                                      tAAACATTATCCCCTTTTTTTTTTTTTTTTTAACATTACAGCGACTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTATCCAAAAAATAAGAATAGATTATTCTATCTACAGCCACAACACTCGACCACg                         >  2:168494/1‑149 (MQ=255)
                                                                                                                        aaCATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCCAGCAACAGCCACAACAGTCGACCACgcg                       >  1:181359/1‑147 (MQ=255)
                                                                                                                         aCATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGGCCACGGGc                      >  1:179317/1‑149 (MQ=255)
                                                                                                                                ccccTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATa               <  1:116812/149‑1 (MQ=255)
                                                                                                                                 tccTTTGTGCCTTATTTTTGACCATTTCCGCTCTTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGTCTATGCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACCCTCGACCACGAGCAAGAATac              <  1:345127/148‑1 (MQ=255)
                                                                                                                                 cccTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATac              >  1:424935/1‑149 (MQ=255)
                                                                                                                                    ttttttCCTGATTTTTTACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACCCGACCACGAGCAAGAATACACg           >  2:73333/1‑149 (MQ=255)
                                                                                                                                          ccTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATACACGCCGATc     <  1:100581/149‑1 (MQ=255)
                                                                                                                                            tGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATACACGCCGATCaa   <  2:424935/149‑1 (MQ=255)
                                                                                                                                             gATTTTTGACCATTTCCGCGATTTGTTAGACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATACTCCCCGGTCAAc  >  2:389913/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TCCAACTCATAACCCCAGCATATAAATCCAGTTGGTAACTTTTATTTAACCTGAAACCAGTTTTATCCACTATTTATAAAATTATGTGAAGCATTTCATAGAAGAAAAATCACTGGCCTAAACATTATCCCCTTTTTGCCTGATTTTTGACCATTTCCGCGATTTGTTACACATTGAAATATCACTTTTGCTGTGCGTAATATGGCTATTCGTTAGCCAAAAAATAAGAAAAGATTATGCAAGCAACAGCCACAACACTCGACCACGAGCAAGAATACACGCCGATCAAC  >  NC_000913/3674549‑3674838

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