Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 396,540 G→A A208T (GCA→ACA)  USA300HOU_0364 → possible lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730396,5400GA98.4% 198.5 / NA 64A208T (GCA→ACA) USA300HOU_0364possible lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/24);  minor base C (1/0);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATG  >  CP000730/396395‑396672
                                                                                                                                                 |                                                                                                                                    
aGCTGATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTCGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGcaca                                                                                                                                    <  2:322588/148‑1 (MQ=255)
    gATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATATATTCACATACAGATTTATATGATTATAAAGCGAACGTTGAAGGCACACaaa                                                                                                                                >  2:572362/1‑148 (MQ=255)
         aTTAGATATCATACCAAAATTAATGTTACAAGGTTCTGCTGACCTATTAAATGCAGTTGCAACTTCTAAAATCAAAGGTGACGAAGAAATTTATTTATATACAGATTTATATGATTTGATATTTGACGTTGAAGGCACACAAAAAAtt                                                                                                                           >  2:622704/1‑148 (MQ=255)
          ttAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAAATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAttt                                                                                                                          >  2:502641/1‑148 (MQ=255)
           tAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTa                                                                                                                         >  1:504382/1‑148 (MQ=255)
              atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                                                                      <  1:561847/148‑1 (MQ=255)
              atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                                                                      >  2:495407/1‑148 (MQ=255)
              atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAAAGTTGAAGGCACCAAAAAAATTTATGa                                                                                                                      >  2:16823/1‑148 (MQ=255)
              atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATAAAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                                                                      <  2:35989/148‑1 (MQ=255)
                    tACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATAAAGACTTATATGCTCTCACAGTTAACGTAGTAATCACAATAAAATTTTATGACTTAtt                                                                                                                <  2:368058/148‑1 (MQ=255)
                       cAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                                                                             <  1:559664/148‑1 (MQ=255)
                       aaaatTTAATGTTAAAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGATGGCACATAAAAAATTTATGACTTATTTaa                                                                                                             <  2:432133/147‑1 (MQ=255)
                        aataTTAATGTTACCAGGTTCTGTTGACCTATTAAATGACGTTCCCACTTCTAGCATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATTTGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                                                            <  2:443254/145‑1 (MQ=255)
                        aaaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                                                            <  1:451871/148‑1 (MQ=255)
                                 gTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                                                                   >  2:294952/1‑148 (MQ=255)
                                 gTTAAAAGGTGCTGTTGACCTATTAAATGAAGTTGCAACTTCGAAATTCACAGGTGAAGAAGAAATTTATTCACATACATATTTACATGATTTTAAAGCGAACGTTGAATGCACACAAAAAATTTATGACTTATTTAAGCCTATTTTa                                                                                                   <  2:70131/148‑1 (MQ=255)
                                    aCAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTagag                                                                                                >  2:414651/1‑148 (MQ=255)
                                     cAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAGTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTGAACCCTATTTTAaaga                                                                                               >  2:252919/1‑148 (MQ=255)
                                     cAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAgaga                                                                                               >  2:153786/1‑148 (MQ=255)
                                      aaGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                                                                              >  1:179151/1‑148 (MQ=255)
                                       aGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaa                                                                                             >  2:484343/1‑148 (MQ=255)
                                        ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                                                            >  1:254036/1‑148 (MQ=255)
                                                   cccaTTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaaa                                                                                 <  1:506774/145‑1 (MQ=255)
                                                    cTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                                                                <  1:332845/148‑1 (MQ=255)
                                                    atATTAAATGAAGTTGCAAGTTCTAAAATCACAGGTCAATAAGAAATTAATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAGGTATGACTTATTTATACCTATTTTAGAGAAAAAAGATAAAAAAt                                                                                <  2:145724/147‑1 (MQ=255)
                                                      aTTAAATGAAGTTGCAACTTCTAAAATGACAGGGGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACAAAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                                                                              <  1:602811/148‑1 (MQ=255)
                                                      aTTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                                                                              >  2:564080/1‑148 (MQ=255)
                                                      aTTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                                                                              <  2:389589/148‑1 (MQ=255)
                                                          aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                                                                          <  1:414651/148‑1 (MQ=255)
                                                          aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATCTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTCTTTTAGAGAAAAAAGATAAAAAATTAAGtg                                                                          >  2:8018/1‑148 (MQ=255)
                                                           aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtga                                                                         >  1:320837/1‑148 (MQ=255)
                                                              aaGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatga                                                                      >  2:459395/1‑148 (MQ=255)
                                                                gTTGCAACTTCTAAAATCACAGGTGAAGAAGAAAGTTATTCACATACAGATTTATATGATTTTAAAGCGAAAGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGata                                                                    >  2:377622/1‑148 (MQ=255)
                                                                          cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGaa                                                          >  2:159888/1‑148 (MQ=255)
                                                                            aaaaaaaaaaaagaagaagaaATTTATTCAAATACAGATTTATATGAAAATAAAGCGAACGTTAAAGCCACACAAAAAATTTATGACTTATTTCAACCTATTTTAGAAAAAAAAAAAAAAAAATTAAGTGATGATATCCAAATGAACt                                                        <  2:266110/136‑1 (MQ=255)
                                                                                 cacaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGCAGGCACACATAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAGAAAATTACGGGCTGATATCCAAATGAACTTAGAt                                                   >  1:162150/1‑148 (MQ=255)
                                                                                  acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTAGATa                                                  >  2:183141/1‑148 (MQ=255)
                                                                                     ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAg                                               <  1:588180/148‑1 (MQ=255)
                                                                                      gTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGGTGATATCCAAATGAACTTCGATAAAGt                                              >  1:38029/1‑148 (MQ=255)
                                                                                        gaagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCCAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCGAAATGAACTTCGATAAAGTGa                                            >  1:436753/1‑148 (MQ=255)
                                                                                          agaagaaATTTATTCACATACAGAATTATATGATTTTAAAGCGAACATTGAAGGCACACAAAAAATTTATGACGTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAAt                                          <  1:142814/148‑1 (MQ=255)
                                                                                             agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATGTCCAAATGAACTTCGATAAAGTGAATCaa                                       >  1:519305/1‑148 (MQ=255)
                                                                                             agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCaa                                       <  1:319855/148‑1 (MQ=255)
                                                                                             agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCaa                                       >  2:581315/1‑148 (MQ=255)
                                                                                              gaatataatgcagaTACAGACTTACATGATTTCAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAat                                      <  1:615853/135‑1 (MQ=255)
                                                                                              gaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAat                                      >  2:138588/1‑148 (MQ=255)
                                                                                               aaaaTTATTCACAAACAGATTTATATGATGTTAAAGCGAACGTTGAAAGAAAAAAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAatt                                     <  2:399474/148‑1 (MQ=255)
                                                                                               aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAGGTGTTGCTTTACAAATGAACTTCGATAAAGTGAATCAatt                                     >  1:238587/1‑148 (MQ=255)
                                                                                                        cacaTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATaa                            <  1:212860/148‑1 (MQ=255)
                                                                                                            tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAatat                        >  2:249032/1‑148 (MQ=255)
                                                                                                            tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAatat                        >  2:249301/1‑148 (MQ=255)
                                                                                                            tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAatat                        <  1:502641/148‑1 (MQ=255)
                                                                                                             aCAGATTTATATGATTTTAAAACGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAAtata                       <  1:217047/148‑1 (MQ=255)
                                                                                                                ggatcaaatGATTTTAAAGGCAAGAGCGAAGGCCCTCAAACACTTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGACGATATCCAAATGAACTTCGATGAAGTGAATCAATTATTGGATAAAAATAAAg                    >  1:465347/8‑148 (MQ=255)
                                                                                                                       atGATTTTAAAGCGAACGTTGAAGGCACGCAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTAGAAATAGTAAGACAATTATTGGATAAATATAAAGATaacaa             >  1:389275/1‑148 (MQ=255)
                                                                                                                       atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATaacaa             <  1:89929/148‑1 (MQ=255)
                                                                                                                         gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAACAAGATGCAAAATTAAGTGATGATATCCAAATGAACTTCGGTAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg           >  1:33845/1‑148 (MQ=255)
                                                                                                                         gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg           >  2:168347/1‑148 (MQ=255)
                                                                                                                          aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGCTGATATCCAAATGAACTTCGATCAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg          >  2:456156/1‑148 (MQ=255)
                                                                                                                          aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAGGATAACAAcgg          >  1:126438/1‑148 (MQ=255)
                                                                                                                          aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg          >  1:41895/1‑148 (MQ=255)
                                                                                                                          aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg          >  1:245582/1‑148 (MQ=255)
                                                                                                                                 aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTAt   >  1:218789/1‑148 (MQ=255)
                                                                                                                                 aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTAt   >  1:344557/1‑148 (MQ=255)
                                                                                                                                  gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATAAATTATTGGATAAATATAAAGATAACAACGGCGGTTATg  >  1:572955/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
AGCTGATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATG  >  CP000730/396395‑396672

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: