Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,166,794 T→C 6.0% F139F (TTT→TTC USA300HOU_1093 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,166,7940TC 6.0% 207.6 / 2.9 67F139F (TTT→TTCUSA300HOU_1093hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ATGTCGATTTGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTCAAATGGTGTCTTT  >  CP000730/1166658‑1166932
                                                                                                                                        |                                                                                                                                          
aTGTCGATTTGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACg                                                                                                                                        >  1:235627/1‑141 (MQ=255)
   tCGATTTGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACgtgt                                                                                                                                     >  1:139223/1‑141 (MQ=255)
   tCGATTTGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACgtgt                                                                                                                                     <  2:504337/141‑1 (MQ=255)
        ttGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAAt                                                                                                                                <  2:283740/141‑1 (MQ=255)
         tGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAAAGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAAtt                                                                                                                               <  2:183374/141‑1 (MQ=255)
            cAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGa                                                                                                                            >  2:413813/1‑141 (MQ=255)
            cAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGa                                                                                                                            >  2:177578/1‑141 (MQ=255)
             agagGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGaa                                                                                                                           >  1:241882/1‑141 (MQ=255)
                   tGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGt                                                                                                                     <  2:538076/141‑1 (MQ=255)
                        ttttAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCt                                                                                                               <  1:332873/141‑1 (MQ=255)
                          ttAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTCGACGTGTCTAATTTGAATTTAGTAATTCTa                                                                                                              >  1:702821/1‑141 (MQ=255)
                          ttAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTCGACGTGTCTAATTTGAATTTAGTAATTCTa                                                                                                              >  1:702795/1‑141 (MQ=255)
                            agagCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTAcc                                                                                                            <  2:563828/141‑1 (MQ=255)
                             gagCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCa                                                                                                           <  2:185625/141‑1 (MQ=255)
                             gagCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCa                                                                                                           >  1:123689/1‑141 (MQ=255)
                              agCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCGATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCaa                                                                                                          >  2:98323/1‑141 (MQ=255)
                                 aTTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGt                                                                                                       >  1:80582/1‑141 (MQ=255)
                                      aTTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATc                                                                                                  <  1:470243/141‑1 (MQ=255)
                                           aTGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAg                                                                                             >  2:425603/1‑141 (MQ=255)
                                                  ttAGTTTACTTTGCGATCTATCACGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTc                                                                                      <  2:354539/141‑1 (MQ=255)
                                                    aGTTTACTTTGCGATCTATCAAGCAGTCAATGACAAAATCGACATTACTGATATAACTCTCGACACTTTCCAACCCGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcac                                                                                    <  1:664512/141‑1 (MQ=255)
                                                    aGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcac                                                                                    <  1:754868/141‑1 (MQ=255)
                                                    aGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcac                                                                                    <  1:656497/141‑1 (MQ=255)
                                                    aGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcac                                                                                    <  2:490811/141‑1 (MQ=255)
                                                    aGTTTACTTTGCGATCTATAAAAAAAAAAAAAAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcac                                                                                    <  2:757407/141‑1 (MQ=255)
                                                     gTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTcaca                                                                                   >  2:192601/1‑141 (MQ=255)
                                                              gCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAAc                                                                          <  1:413813/141‑1 (MQ=255)
                                                               cGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACt                                                                         >  1:232076/1‑141 (MQ=255)
                                                               cGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACt                                                                         >  1:350635/1‑141 (MQ=255)
                                                               cGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACt                                                                         >  1:109836/1‑141 (MQ=255)
                                                                 aTCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACtat                                                                       <  1:97932/141‑1 (MQ=255)
                                                                   cTATCAAGCAGAAAAAAAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAg                                                                     <  2:76663/141‑1 (MQ=255)
                                                                     aTCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGaa                                                                   <  1:249646/141‑1 (MQ=255)
                                                                               aaaTGCAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTAAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTa                                                         >  2:307226/1‑141 (MQ=255)
                                                                               aaaTGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTa                                                         >  2:51639/1‑141 (MQ=255)
                                                                                  tGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACaa                                                      <  1:37506/141‑1 (MQ=255)
                                                                                  tGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACaa                                                      >  1:302727/1‑141 (MQ=255)
                                                                                  tGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACaa                                                      <  1:92799/141‑1 (MQ=255)
                                                                                   gAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACaaa                                                     <  2:39878/141‑1 (MQ=255)
                                                                                    aaaaaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTACTCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAAt                                                    <  1:303979/141‑1 (MQ=255)
                                                                                    aaaaaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAAt                                                    >  1:659345/1‑141 (MQ=255)
                                                                                     aaaaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACTTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATATAATGAGAAGTTACAAATa                                                   >  2:190353/1‑141 (MQ=255)
                                                                                     aaaaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATa                                                   <  2:721193/141‑1 (MQ=255)
                                                                                     aaaaaTGGACATTAATGATATTAATAGCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTGAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTAGAAATa                                                   >  1:509585/1‑141 (MQ=255)
                                                                                      aaaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATaa                                                  >  1:476153/1‑141 (MQ=255)
                                                                                       aaaTGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATaaa                                                 <  1:549900/141‑1 (MQ=255)
                                                                                           ggACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCt                                             <  2:401308/141‑1 (MQ=255)
                                                                                           ggACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCt                                             <  1:84343/141‑1 (MQ=255)
                                                                                            gACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTa                                            <  2:581893/141‑1 (MQ=255)
                                                                                            gACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTa                                            >  1:591131/1‑141 (MQ=255)
                                                                                              cATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATg                                          <  1:420020/141‑1 (MQ=255)
                                                                                              cATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATg                                          >  2:186123/1‑141 (MQ=255)
                                                                                               aTTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATga                                         <  2:147432/141‑1 (MQ=255)
                                                                                                     gATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATtat                                   >  2:352167/1‑141 (MQ=255)
                                                                                                     gATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATACTACCAACGTTAATCATTAGCTGGCTCACAACATTTCACGCTAGAATGAGAAGTTACAAATAAAATCTATGAGATtat                                   >  1:272170/1‑141 (MQ=255)
                                                                                                     atcatTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATtat                                   >  2:177494/4‑141 (MQ=255)
                                                                                                      atatTAATATCGACAATTTCCAATCTGTCTTTTTCGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATtata                                  <  2:702795/141‑1 (MQ=255)
                                                                                                      atatTAATATCGACAATTTCCAATCTGTCTTTTTCGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATtata                                  <  2:702821/141‑1 (MQ=255)
                                                                                                        atTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATAcc                                >  1:368279/1‑141 (MQ=255)
                                                                                                         ttAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCt                               <  1:450832/141‑1 (MQ=255)
                                                                                                             tatCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAg                           >  1:485155/1‑141 (MQ=255)
                                                                                                                     ggctaccttgttacTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAAc                   <  1:321734/128‑1 (MQ=255)
                                                                                                                        tCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACAtt                >  2:601069/1‑141 (MQ=255)
                                                                                                                              cTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCatttgaatgttg          <  1:371400/141‑12 (MQ=255)
                                                                                                                              cTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCatttgaatgttg          <  1:371216/141‑12 (MQ=255)
                                                                                                                                gTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAAGGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTAAAATGgt        >  2:305573/1‑141 (MQ=255)
                                                                                                                                 tCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTCAAATGgtg       <  1:620950/141‑1 (MQ=255)
                                                                                                                                   ttttttGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTCAAATGGTGTc     <  1:138243/141‑1 (MQ=255)
                                                                                                                                      tttGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTCAAATGGTGTCttt  <  1:195426/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
ATGTCGATTTGACAGAGGTTGTCGTTTTAGAGCATTTAATTCAATGGTTCTTAGTTTACTTTGCGATCTATCAAGCAGTAAATGAAAAAATGGACATTAATGATATTAATATCGACAATTTCCAATCTGTCTTTTTTGACGTGTCTAATTTGAATTTAGTAATTCTACCAACGTTAATCATTAGCTGGGTCACAATATTTAACTATAGAATGAGAAGTTACAAATAAAATCTATGAGATTATACCTTCAGACACCAACATTCAAATGGTGTCTTT  >  CP000730/1166658‑1166932

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