Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,000,948 C→T 63.3% T11I (ACA→ATA)  USA300HOU_0955 → possible glutaredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,000,9480CT63.3% 44.2 / 60.3 67T11I (ACA→ATA) USA300HOU_0955possible glutaredoxin
Reads supporting (aligned to +/- strand):  ref base C (14/10);  new base T (22/20);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

AATAAATAGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACTGATGAAA  >  CP000730/1000812‑1001082
                                                                                                                                         |                                                                                                                                      
aataaataGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTTTTGTTCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGATTGAGATGTATGGTAACATCATTTAATTCACCAAATTGCACATCt                                                                                                                                    >  2:1298/1‑141 (MQ=255)
   aaataGTAACGATTAAGAATTGTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGc                                                                                                                                 <  1:773101/141‑1 (MQ=255)
    aataGTAACGATTAAGACTTTTAGAAAAAAATGGTTTTTCTTGATCTATAACTAAATTTTAAGCCATGAATTTTATTCATTCAATTAGGAGATTGAGATGTATGGTAACATAATTTAGTTTCCTATGTTGCatatttagct                                                                                                                                >  2:740783/1‑135 (MQ=255)
      taGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGt                                                                                                                              >  1:506049/1‑141 (MQ=255)
       aGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTa                                                                                                                             >  2:657943/1‑141 (MQ=255)
       aGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTa                                                                                                                             >  1:570877/1‑141 (MQ=255)
        gTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTaa                                                                                                                            >  1:746012/1‑141 (MQ=255)
        gTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTaa                                                                                                                            >  1:81620/1‑141 (MQ=255)
           aCGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGc                                                                                                                         <  1:53240/141‑1 (MQ=255)
           aCGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATAACAAAAAATTAAGGTTATTCATTCCCCTCATCAGCCACAAGGGATGGAAACAGGATGTACTTCACCAAGTTGCACATCTTGCCGTAAAGc                                                                                                                         >  2:359875/1‑141 (MQ=255)
               ttAAGAATTTTAGAAAAAAATGTTTTTTTTTTTTTTTTTACTAAATTTTAGGCAATGAATTTTGTTTATTCTGATAGGAGAGTGAGATGTATGTTTTTATAATTTAATTCACCATTTTGTATTTCTTGCAGTAAAgttaat                                                                                                                     >  2:727696/1‑136 (MQ=255)
                tAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAg                                                                                                                    >  1:414736/1‑141 (MQ=255)
                 aaGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGc                                                                                                                   >  2:286283/1‑141 (MQ=255)
                 aaGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGc                                                                                                                   >  2:602708/1‑141 (MQ=255)
                  aGAATTTTAGAGAAAAAGTTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGCGAGTGCGATGTATGGTAACATTAGTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCa                                                                                                                  <  2:70451/141‑1 (MQ=255)
                  aGAATTATAAAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCa                                                                                                                  <  2:674515/141‑1 (MQ=255)
                   gAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGGGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGTCGTAAAGCGAAAGCAt                                                                                                                 >  2:310037/1‑141 (MQ=255)
                   gAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCAt                                                                                                                 >  1:27570/1‑141 (MQ=255)
                      ttttAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGt                                                                                                              >  1:725071/1‑141 (MQ=255)
                      ttttAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGt                                                                                                              <  2:669535/141‑1 (MQ=255)
                           gAAAAAAATGTTTTGTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTAGTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACaa                                                                                                         >  2:328765/1‑141 (MQ=255)
                            aaaaaaaTGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAg                                                                                                        >  1:297650/1‑141 (MQ=255)
                                  aTGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATg                                                                                                  >  2:718352/1‑141 (MQ=255)
                                  aTGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTCCACCTTGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATg                                                                                                  <  2:759609/141‑1 (MQ=255)
                                  aTGTTTTTTCTTGATCTATAACTAAATCTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATg                                                                                                  <  1:666519/141‑1 (MQ=255)
                                    gTTTTTTCTTGATCTATAACTAAATTTTAAGCACGCACGTTCACTCTTGGAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGAc                                                                                                <  2:164042/141‑1 (MQ=255)
                                                      actgacttttAAGCAATGAAGTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGt                                                                              <  2:572577/135‑1 (MQ=255)
                                                      aCTAAATTTTAAGCAATGTATTTTATTCATCCCAATAGGAGAGTGAGATGTATGGTAACATTAGTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGt                                                                              <  2:744180/141‑1 (MQ=255)
                                                      aCTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGt                                                                              <  2:251169/141‑1 (MQ=255)
                                                               tAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATAttttt                                                                     >  1:312299/1‑141 (MQ=255)
                                                                aaGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATAtttttt                                                                    <  1:443895/141‑1 (MQ=255)
                                                                   cAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTCTGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGCAAGCATGGTTACAAGAACATGACATTCCGTATACCGAGCGTAATATTTTTTCTg                                                                 <  1:512439/141‑1 (MQ=255)
                                                                   cAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTg                                                                 <  1:695987/141‑1 (MQ=255)
                                                                   cAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTg                                                                 <  1:146998/141‑1 (MQ=255)
                                                                    aaTGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGa                                                                <  1:324292/141‑1 (MQ=255)
                                                                     aTGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGaa                                                               >  1:138896/1‑141 (MQ=255)
                                                                        gctttgagtcgttcaAATAGGAGAGTGAGTTGCGGGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACAt                                                            <  1:371952/130‑1 (MQ=255)
                                                                        aaTTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACAt                                                            <  1:359875/141‑1 (MQ=255)
                                                                          ttttATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCCACAAGTTGCAAAACTTGCCATAAAGCGAAAGCATGGTTACCAGAAAATGACATTCCGTATACGAAACGAGATATTGTTTCTTAACAttt                                                          >  2:599817/1‑141 (MQ=255)
                                                                          ttttAGTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATGTACTCCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACAttt                                                          <  1:84639/141‑1 (MQ=255)
                                                                               ttcattcaAAGAAGAGAGTGACTTGATGGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACaa                                                     <  2:197254/141‑1 (MQ=255)
                                                                                 aatttaAATCGGAGAGTTAGACTTATGGTAACATTATTTACTTAACACACACGCACATCTTCCCCTCACGCTCTCCACTTAATACTAAGACATGTCATACCGCATACGGAGCGTAATATGCCTTCTGAACATTTAACAAtt                                                   <  2:197728/140‑1 (MQ=255)
                                                                                  gttcaAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATtg                                                  >  2:172816/2‑141 (MQ=255)
                                                                                   ttcaAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATtga                                                 <  2:570877/141‑1 (MQ=255)
                                                                                   ttcaAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATtga                                                 >  1:667332/1‑141 (MQ=255)
                                                                                    tcaAATAGGATAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATtgat                                                <  2:381557/141‑1 (MQ=255)
                                                                                            ggAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACCATTGATGAAATTa                                         >  1:725978/1‑141 (MQ=255)
                                                                                            ggAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTa                                         >  2:82752/1‑141 (MQ=255)
                                                                                            ggAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTa                                         >  2:161927/1‑141 (MQ=255)
                                                                                             gagaGTGAGATGTATGATAACATTATTTAATTCACCATGCTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTaa                                        <  2:201253/141‑1 (MQ=255)
                                                                                               gagTGAGATGTATGGTAACATTATTAGCTTCACCAAGTTGCATATTTTGCCGGAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGc                                      <  2:320565/141‑1 (MQ=255)
                                                                                                agTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATCTTTTTTCTGAACATTTAACAATTGATGAAATTAAGCa                                     >  1:144661/1‑141 (MQ=255)
                                                                                                agTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCa                                     <  1:630165/141‑1 (MQ=255)
                                                                                                agTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCa                                     <  1:718352/141‑1 (MQ=255)
                                                                                                     gacgTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGGAATTAAGCAAatat                                >  2:597592/4‑141 (MQ=255)
                                                                                                     gaTGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAatat                                <  1:435588/141‑1 (MQ=255)
                                                                                                      aTGTATGGTAACATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATAtt                               >  1:606397/1‑141 (MQ=255)
                                                                                                                 cATTATTTACTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTg                    >  1:327888/1‑141 (MQ=255)
                                                                                                                     attTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGa                >  2:491693/1‑141 (MQ=255)
                                                                                                                          cTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTCAGCAAATATTAAAAATGAATGAAGACGGTa           >  2:503417/1‑141 (MQ=255)
                                                                                                                          cTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTa           >  2:306160/1‑141 (MQ=255)
                                                                                                                          cTTCACCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGAAGGTa           >  2:410275/1‑141 (MQ=255)
                                                                                                                            tCACCAAGTTGCACATTTTGCCGTAAAGCGAAAGCATGGTTACAAGAACAGGACATTCCGTATACGGAGCGTAATATTCTTTCTGAACATTTAACAATTGAGGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACt         <  1:544240/141‑1 (MQ=255)
                                                                                                                            tCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACt         >  2:251060/1‑141 (MQ=255)
                                                                                                                              aCCAAGTTGCATATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACtga       >  2:758811/1‑141 (MQ=255)
                                                                                                                              aCCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACtga       >  2:84023/1‑141 (MQ=255)
                                                                                                                                aaaGTTGCATATTTTGCCGAAAATCGAAAGCATGGTTACAACAACATCACATTCCGTATACGGAGCGTAATATTTTTCCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACtgatg     <  1:20963/140‑1 (MQ=255)
                                                                                                                                   gTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGACGACGGTACTGATGaaa  >  2:492160/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                      
AATAAATAGTAACGATTAAGAATTTTAGAAAAAAATGTTTTTTCTTGATCTATAACTAAATTTTAAGCAATGAATTTTATTCATTCAAATAGGAGAGTGAGATGTATGGTAACATTATTTACTTCACCAAGTTGCACATCTTGCCGTAAAGCGAAAGCATGGTTACAAGAACATGACATTCCGTATACGGAGCGTAATATTTTTTCTGAACATTTAACAATTGATGAAATTAAGCAAATATTAAAAATGACTGAAGACGGTACTGATGAAA  >  CP000730/1000812‑1001082

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: