Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 34,112 Δ1 bp 100% intergenic (+74/‑188) carB → / → caiF carbamoyl‑phosphate synthase large subunit/cai operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,1120T.81.6% 2313.4 / inf 817intergenic (+74/‑188)carB/caiFcarbamoyl‑phosphate synthase large subunit/cai operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (62/85);  new base . (289/378);  total (352/465)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 898 total aligned reads displayed.

ACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAACGA  >  NC_000913/33964‑34260
                                                                                                                                                    |                                                                                                                                                    
aCCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt                                                                                                                                                      >  4:110971/1‑149 (MQ=255)
aCCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt                                                                                                                                                      >  1:3204686/1‑149 (MQ=255)
  cgcgATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttt                                                                                                                                                    >  4:2774743/1‑149 (MQ=255)
     gATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttt                                                                                                                                                 <  1:562948/149‑1 (MQ=255)
     gATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttt                                                                                                                                                 >  3:1236984/1‑149 (MQ=255)
        ggCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttt                                                                                                                                              >  1:7046663/1‑149 (MQ=255)
          cgcTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGc                                                                                                                                            <  1:6020039/149‑1 (MQ=255)
          caacctaTGCCGATCCGAATGAAAAGTTATTTTCGTGGCAGGAAATGCACGCACAAATCAAAAAATACGGTGTCTTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGc                                                                                                                                            <  1:6641536/143‑1 (MQ=255)
             tGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCa                                                                                                                                         <  4:2322002/149‑1 (MQ=255)
              gAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCAtt                                                                                                                                       <  3:3932879/149‑1 (MQ=255)
              gAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCAtt                                                                                                                                       >  1:6887751/1‑149 (MQ=255)
                aTGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTaa                                                                                                                                     <  2:4076203/149‑1 (MQ=255)
                 tGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAg                                                                                                                                    >  4:896058/1‑149 (MQ=255)
                             tGAAAAATTAATTGCTGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCGGCTAATTTGATCGAATGACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTGATTAAATCTTTTgg                                                                                                                        >  4:3792740/1‑149 (MQ=255)
                              gAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAGTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                                                                                        <  3:3687007/149‑1 (MQ=255)
                                aaaaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGa                                                                                                                     <  2:1247401/149‑1 (MQ=255)
                                  aaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAAt                                                                                                                   >  1:2788084/1‑149 (MQ=255)
                                   aGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                                                                                  >  1:5128305/1‑149 (MQ=255)
                                    gTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                                                                                  <  1:3500644/149‑1 (MQ=255)
                                          tCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttt                                                                                                           <  1:6307558/149‑1 (MQ=255)
                                            ggTGCAGGAAATGCAGGCATAGATCAAAGAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATATCTAATACGGTGCACTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                                                                                         <  3:5875232/149‑1 (MQ=255)
                                            ggTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                                                                                         <  1:5780251/149‑1 (MQ=255)
                                             gTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttttt                                                                                                                                               >  1:5375029/1‑111 (MQ=255)
                                             gTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttttt                                                                                                                                               <  2:5375029/111‑1 (MQ=255)
                                              tGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                                                                                       <  3:6376299/149‑1 (MQ=255)
                                              tGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                                                                                       <  2:4896758/149‑1 (MQ=255)
                                              tGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                                                                                       <  4:6096280/149‑1 (MQ=255)
                                              tGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                                                                                       <  2:5901762/149‑1 (MQ=255)
                                               gCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                                                                                       <  4:67177/149‑1 (MQ=255)
                                               gCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgat                                                                                                      >  4:7003594/1‑149 (MQ=255)
                                                 aGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatga                                                                                                    <  2:938011/149‑1 (MQ=255)
                                                  ggAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatga                                                                                                    >  2:2098343/1‑149 (MQ=255)
                                                   gAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGAc                                                                                                   >  2:386887/1‑149 (MQ=255)
                                                          acgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                                                                                           <  4:1701335/149‑1 (MQ=255)
                                                          acgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                                                                                           >  1:1081274/1‑149 (MQ=255)
                                                            gcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGa                                                                                         <  2:4296659/149‑1 (MQ=255)
                                                             cacaGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGa                                                                                         <  3:819114/149‑1 (MQ=255)
                                                                  aTCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAg                                                                                    <  1:1240306/149‑1 (MQ=255)
                                                                    cAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAg                                                                                    >  2:7100499/1‑146 (MQ=255)
                                                                     aaaTAATAGCGTGTCATGGCAGATATTTTTCATCCGCTTATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATTTTTTGGGGAAACGATATTTTTGTTGACATAAGCAGGATTTAGCTc                                                                                 >  4:5871794/1‑149 (MQ=255)
                                                                     aaaTAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTca                                                                                <  2:365232/149‑1 (MQ=255)
                                                                      aataatAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAt                                                                                        >  2:4420291/1‑140 (MQ=255)
                                                                      aataatAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcac                                                                               >  1:4729411/1‑149 (MQ=255)
                                                                       ataataGCGTGTCATGGCAGATATTTTTCCTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCACTAAGTAAATCTTTTGGGCAATCGATATTTTTGCTGACCTAAGCAGCATTTAGCTcaca                                                                              <  4:254658/149‑1 (MQ=255)
                                                                       ataataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcaca                                                                              <  1:1544336/149‑1 (MQ=255)
                                                                       ataataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcaca                                                                              >  1:2424525/1‑149 (MQ=255)
                                                                       ataataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcaca                                                                              <  3:2704475/149‑1 (MQ=255)
                                                                       ataataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcaca                                                                              <  2:6624941/149‑1 (MQ=255)
                                                                        taataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttt                                                                                                                                                 >  2:3457574/1‑82 (MQ=255)
                                                                        taataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                                                                                       >  2:5574802/1‑139 (MQ=255)
                                                                         aataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCacac                                                                             <  3:417270/149‑1 (MQ=255)
                                                                          ataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACtt                                                                           >  3:302749/1‑149 (MQ=255)
                                                                               gtgtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGa                                                                      <  2:6281134/149‑1 (MQ=255)
                                                                                       gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCAtt                                                                                                                                       >  4:4047696/1‑76 (MQ=255)
                                                                                       gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGaa                                                               <  1:7261022/148‑1 (MQ=255)
                                                                                       gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGaa                                                               >  2:7261022/1‑148 (MQ=255)
                                                                                       gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAg                                                              <  1:229886/149‑1 (MQ=255)
                                                                                            tatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                                         <  2:6018532/149‑1 (MQ=255)
                                                                                              tttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATa                                                       >  2:1945866/1‑149 (MQ=255)
                                                                                               ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGc                                                          >  3:4762931/1‑145 (MQ=255)
                                                                                                tttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACt                                                     >  3:5382510/1‑149 (MQ=255)
                                                                                                     tCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGa                                                <  2:5058930/149‑1 (MQ=255)
                                                                                                      ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGa                                                <  1:4900385/149‑1 (MQ=255)
                                                                                                      ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATa                                                                                                <  4:3507151/100‑1 (MQ=255)
                                                                                                      ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGat                                               <  3:2024816/149‑1 (MQ=255)
                                                                                                       cGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGata                                              <  3:2407033/149‑1 (MQ=255)
                                                                                                       cGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGata                                              >  2:2745816/1‑149 (MQ=255)
                                                                                                        gCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                                                                                        >  2:1750027/1‑74 (MQ=255)
                                                                                                        gCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGatat                                             >  1:1892720/1‑149 (MQ=255)
                                                                                                        gCTAATTTGATCGAATAACTAATACGGTTCTCAGATGAGGACCGTTTTTTTTTGCCCTTTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGata                                              >  2:1055619/1‑149 (MQ=255)
                                                                                                            aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATc                                          <  3:6752484/149‑1 (MQ=255)
                                                                                                               tGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCcac                                       <  3:6587924/149‑1 (MQ=255)
                                                                                                                gATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACaa                                     >  1:4284386/1‑149 (MQ=255)
                                                                                                                gATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACaa                                     >  4:3270245/1‑149 (MQ=255)
                                                                                                                 aTCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGatat                                                                                                             <  1:3740132/77‑1 (MQ=255)
                                                                                                                 aTCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATc                                                  <  4:5583506/135‑1 (MQ=255)
                                                                                                                    gAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAtttt                                 >  3:5056250/1‑149 (MQ=255)
                                                                                                                     aaTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                                         >  4:4766091/1‑124 (MQ=255)
                                                                                                                     aaTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTa                                >  2:6942808/1‑149 (MQ=255)
                                                                                                                      aTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTaa                               <  3:5650119/149‑1 (MQ=255)
                                                                                                                        aaCTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAata                             >  4:649746/1‑149 (MQ=255)
                                                                                                                          cTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatga                                                                                                    <  4:1875845/76‑1 (MQ=255)
                                                                                                                          cTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATg                           <  4:5273313/149‑1 (MQ=255)
                                                                                                                             aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTggg                                                                                                                       <  1:4564430/54‑1 (MQ=255)
                                                                                                                             aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGcc                        >  2:4700571/1‑149 (MQ=255)
                                                                                                                              tACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCt                       >  3:4467124/1‑149 (MQ=255)
                                                                                                                                cGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCtt                      <  4:1341835/149‑1 (MQ=255)
                                                                                                                                 ggTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGt                    <  4:4988221/149‑1 (MQ=255)
                                                                                                                                  gTTCTCTGATGAGGTCCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGtt                   <  4:4526105/149‑1 (MQ=255)
                                                                                                                                  gTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATcc                                         <  3:3399109/128‑1 (MQ=255)
                                                                                                                                     ctctGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATc                                                                        >  3:988315/1‑93 (MQ=255)
                                                                                                                                      tctGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAAt               <  4:302749/149‑1 (MQ=255)
                                                                                                                                      tctGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAAt               <  2:1444709/149‑1 (MQ=255)
                                                                                                                                       ctGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAAt               <  1:2098343/149‑1 (MQ=255)
                                                                                                                                        tgatgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTg             >  2:3253382/1‑149 (MQ=255)
                                                                                                                                        tgatgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTg             >  2:6783969/1‑149 (MQ=255)
                                                                                                                                          atgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCt           <  3:6160928/149‑1 (MQ=255)
                                                                                                                                           tgaGGTCCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGGCGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCtt          <  1:398952/149‑1 (MQ=255)
                                                                                                                                            gaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTAGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCACACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTc         <  2:2308889/149‑1 (MQ=255)
                                                                                                                                                    tttttttttGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAACGa  <  3:1774791/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAACGA  >  NC_000913/33964‑34260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑