Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 977,528 Δ1 bp coding (1177/1197 nt) lplA → sugar ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM406671977,5260C.98.6% 4429.7 / 44.5 947coding (1175/1197 nt)lplAsugar ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (8/5);  new base . (501/431);  total (510/437)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 972 total aligned reads displayed.

CCTGCAGGGGAAAATGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTAATGGTGTCT  >  AM406671/977378‑977670
                                                                                                                                                    |                                                                                                                                                
   gCAGGGGAAAATGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAAcctt                                                                                                                                               >  1:9822437/1‑147 (MQ=255)
    cAGGGGAAAATGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCttt                                                                                                                                              >  1:1477705/1‑149 (MQ=255)
       gggAAAATGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTgcg                                                                                                                                          <  2:4275797/149‑1 (MQ=255)
          aaaaTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTgcgcgc                                                                                                                                       >  2:6553112/1‑149 (MQ=255)
          aaaaTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTgcgcgc                                                                                                                                       <  1:9171108/149‑1 (MQ=255)
           aaaTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCt                                                                                                                                      <  2:5546735/149‑1 (MQ=255)
            aaTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGATAACTAGTGACCATTTCCTTGCCACCCTCCATTTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACGTTTGCGCGCTg                                                                                                                                     >  2:2306551/1‑149 (MQ=255)
             aTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACGAACCTTTGCGCGCTgg                                                                                                                                    >  1:7027856/1‑149 (MQ=255)
             aTGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTgg                                                                                                                                    >  1:7175795/1‑149 (MQ=255)
              tGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTggg                                                                                                                                   >  2:2497325/1‑149 (MQ=255)
                     aTGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGa                                                                                                                            >  1:303508/1‑149 (MQ=255)
                     aTGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGa                                                                                                                            <  1:3558775/149‑1 (MQ=255)
                      tGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGaa                                                                                                                           <  1:3431155/149‑1 (MQ=255)
                      tGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGaa                                                                                                                           <  1:410234/149‑1 (MQ=255)
                           cAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTAc                                                                                                                      >  2:5630371/1‑149 (MQ=255)
                           cAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTAc                                                                                                                      <  1:3266766/149‑1 (MQ=255)
                            aaGATAACCGGATTGCAGCCGTTGCTGGNATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACa                                                                                                                     >  1:5786530/1‑149 (MQ=255)
                             aGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACaa                                                                                                                    <  1:8289011/149‑1 (MQ=255)
                              gATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAAc                                                                                                                   >  2:8050783/1‑149 (MQ=255)
                              gATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAAc                                                                                                                   >  1:4756349/1‑149 (MQ=255)
                                   ccGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTtata                                                                                                              <  1:2314711/149‑1 (MQ=255)
                                   ccGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTtata                                                                                                              >  1:171301/1‑149 (MQ=255)
                                   ccGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTtata                                                                                                              <  2:8595936/149‑1 (MQ=255)
                                      gATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTa                                                                                                           >  2:104210/1‑149 (MQ=255)
                                      gATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTa                                                                                                           >  1:3153314/1‑149 (MQ=255)
                                       aTTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTaa                                                                                                          >  2:5433693/1‑149 (MQ=255)
                                          gCAGCCCTGGCTGGTATGATGTTATCTTCTAAATAAAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATc                                                                                                       <  2:10269180/149‑1 (MQ=255)
                                             gCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTgg                                                                                                    >  1:10019519/1‑149 (MQ=255)
                                             gCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTgg                                                                                                    <  1:5438173/149‑1 (MQ=255)
                                              ccGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGga                                                                                                   >  1:9495825/1‑149 (MQ=255)
                                               cGTTGCTGGTATGATGTTATCTTCTACATTGAACAATAGTGACCATTTCCTTGCCACCCTCCAACTCATTGCCTGGTTACTTTATTCACACTCTCGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTCCAACTTATACAAAATCtacag                                                                                                  >  2:9332117/1‑145 (MQ=255)
                                               cGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGgag                                                                                                  >  1:3857106/1‑149 (MQ=255)
                                                  tGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGaca                                                                                               <  1:7802493/149‑1 (MQ=255)
                                                   gCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacac                                                                                              <  1:789650/149‑1 (MQ=255)
                                                     tGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacactt                                                                                            <  1:2500388/149‑1 (MQ=255)
                                                      ggTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacactta                                                                                           <  2:6088758/149‑1 (MQ=255)
                                                       gTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacacttac                                                                                          >  1:3940258/1‑149 (MQ=255)
                                                        tATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacacttaca                                                                                         <  2:9867801/149‑1 (MQ=255)
                                                        tATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacacttaca                                                                                         <  1:101960/149‑1 (MQ=255)
                                                          tgatgTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacacttacact                                                                                       >  2:4424526/1‑149 (MQ=255)
                                                           gatgTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGacacttacactt                                                                                      <  2:9519597/149‑1 (MQ=255)
                                                             tgTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGATGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTg                                                                                    >  2:6850608/1‑149 (MQ=255)
                                                               ttATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGc                                                                                  <  1:7670645/149‑1 (MQ=255)
                                                                tATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCa                                                                                 >  2:3050878/1‑149 (MQ=255)
                                                                tATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCa                                                                                 >  1:4229364/1‑149 (MQ=255)
                                                                  tcttctAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCaaa                                                                               >  2:8797195/1‑149 (MQ=255)
                                                                   cttctAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCaaaa                                                                              <  2:2653167/149‑1 (MQ=255)
                                                                      ctAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGAt                                                                           >  2:7758883/1‑149 (MQ=255)
                                                                        aaaTTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGt                                                                         <  1:2819889/149‑1 (MQ=255)
                                                                        aaaTTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGt                                                                         <  1:9297194/149‑1 (MQ=255)
                                                                        aaaTTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGt                                                                         <  2:6733258/149‑1 (MQ=255)
                                                                          aTTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTa                                                                       >  1:9683601/1‑149 (MQ=255)
                                                                          aTTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTa                                                                       <  1:5559611/149‑1 (MQ=255)
                                                                           ttGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTaa                                                                      >  1:5137138/1‑149 (MQ=255)
                                                                             gAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATa                                                                    >  1:600404/1‑149 (MQ=255)
                                                                              aaaaaTAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATAc                                                                   >  2:8960370/1‑149 (MQ=255)
                                                                               aaaaTAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATAcc                                                                  <  1:900407/149‑1 (MQ=255)
                                                                                aaaTAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCt                                                                 >  1:6656090/1‑149 (MQ=255)
                                                                                 aaTAGTGACCATTTCCTTGCCACCCTCCAATATATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGAGAGTTGAAAGTACAACTTATACTAAATCTGGAGACACTAAAAAAAAAAAAAAAGATGTTAATACCtt                                                                <  2:1977878/149‑1 (MQ=255)
                                                                                  aTAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTa                                                               >  1:529421/1‑149 (MQ=255)
                                                                                    aGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGTCTCTGGAGAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACATTAgg                                                             >  2:7130701/1‑149 (MQ=255)
                                                                                    aGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAgg                                                             >  2:5837519/1‑149 (MQ=255)
                                                                                      tGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCt                                                           <  2:83253/149‑1 (MQ=255)
                                                                                        aCCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTg                                                         >  2:9188424/1‑149 (MQ=255)
                                                                                        aCCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTg                                                         <  1:6875739/149‑1 (MQ=255)
                                                                                         ccATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGa                                                        >  1:6451816/1‑149 (MQ=255)
                                                                                         ccATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGa                                                        <  1:10389019/149‑1 (MQ=255)
                                                                                           aTTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAAc                                                      >  2:1783605/1‑149 (MQ=255)
                                                                                            tttCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAAcc                                                     >  2:1320223/1‑149 (MQ=255)
                                                                                              tCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCa                                                   >  1:6223082/1‑149 (MQ=255)
                                                                                               ccTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCaa                                                  >  2:3435837/1‑149 (MQ=255)
                                                                                                  tGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTg                                               >  2:6328215/1‑149 (MQ=255)
                                                                                                   gCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTgg                                              >  2:7645847/1‑149 (MQ=255)
                                                                                                    ccACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGt                                             <  2:2732337/149‑1 (MQ=255)
                                                                                                    ccACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGt                                             >  2:1649801/1‑149 (MQ=255)
                                                                                                     cACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTa                                            <  1:5033265/149‑1 (MQ=255)
                                                                                                       cccTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACTCTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACa                                          <  2:5351985/149‑1 (MQ=255)
                                                                                                       cccTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACa                                          >  1:747487/1‑149 (MQ=255)
                                                                                                        ccTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACaa                                         >  1:1243555/1‑149 (MQ=255)
                                                                                                        ccTCCAATTTATTGACTGGTTACGTTATTCGGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACGGCGGATACTAAATCGGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACaa                                         <  2:4627371/149‑1 (MQ=255)
                                                                                                             aaTTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGAc                                    <  1:3751386/149‑1 (MQ=255)
                                                                                                             aaTTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGAc                                    <  1:4504862/149‑1 (MQ=255)
                                                                                                              aTTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGAcc                                   >  2:8203419/1‑149 (MQ=255)
                                                                                                               tttATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCt                                  >  1:8042925/1‑149 (MQ=255)
                                                                                                                 tATTGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGc                                >  1:8869066/1‑149 (MQ=255)
                                                                                                                   ttGACTGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCaa                              >  1:9990676/1‑149 (MQ=255)
                                                                                                                        tGGTTACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGa                         >  1:2161369/1‑149 (MQ=255)
                                                                                                                        tGGTAACTTTATTCAGACCCTGGTGAAAGCTATGAGCGCAAGAAAGTTAAAGGTATAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTATGCTTGAACCCAAGTGGTACAAAAGACCTTCAAAAAGa                         >  2:3985856/1‑149 (MQ=255)
                                                                                                                            tACTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGAtttt                     <  1:1198771/149‑1 (MQ=255)
                                                                                                                            tACTTTATTCAGACTCTGGACAAACCTTAGCGCGCTAAAGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTATGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTAATAAAAAAGAACTGCAAAAAGAtttt                     <  2:3756256/149‑1 (MQ=255)
                                                                                                                             aCTTTATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTg                    >  2:7223632/1‑149 (MQ=255)
                                                                                                                               tttATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGCTACAACTGATACTAAATCGGGAGACACTTGCACTTTGGCAAGAGATGTTAATACCTTAGGCTAGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTgct                  >  2:7571372/1‑147 (MQ=255)
                                                                                                                                 tATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGttt                >  2:4860281/1‑149 (MQ=255)
                                                                                                                                 tATTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTAAAACTTATAATAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGttt                >  2:3200644/1‑149 (MQ=255)
                                                                                                                                  aTTCAGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTc               >  2:5065385/1‑149 (MQ=255)
                                                                                                                                      aGACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTa           <  1:1962475/149‑1 (MQ=255)
                                                                                                                                       gACTCTGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTaa          <  1:8137449/149‑1 (MQ=255)
                                                                                                                                           ctGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTAATGgt      <  1:2656696/149‑1 (MQ=255)
                                                                                                                                            tGGACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTAATGgtg     <  1:8579077/149‑1 (MQ=255)
                                                                                                                                             ggACAAACCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTAATGgtgt    >  2:5141068/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCTGCAGGGGAAAATGGTTTGATGGTTCAAGATAACCGGATTGCAGCCGTTGCTGGTATGATGTTATCTTCTAAATTGAAAAATAGTGACCATTTCCTTGCCACCCTCCAATTTATTGACTGGTTACTTTATTCAGACTCTGGACAAACCCTTTGCGCGCTGGGGAGTTGAAGGTACAACTTATACTAAATCTGGAGACACTTACACTTTGGCAAAAGATGTTAATACCTTAGGCTTGAACCCAAGTGGTACAAAAGACCTGCAAAAAGATTTTGGTTTCCGTAATGGTGTCT  >  AM406671/977378‑977670

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: