Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 962,317 A→T 7.8% T81S (ACC→TCC)  msbA → fused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273962,3170AT 7.8% 670.0 / 37.9 217T81S (ACC→TCC) msbAfused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base A (89/111);  new base T (6/11);  total (95/122)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 217 total aligned reads displayed.

AGCGTTAATCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTATTACC  >  CP009273/962172‑962463
                                                                                                                                                 |                                                                                                                                                  
aGCGTTAATCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCAcc                                                                                                                                                  <  2:2063984/148‑1 (MQ=255)
    ttAATCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCt                                                                                                                                              <  1:1062970/148‑1 (MQ=255)
     tAATCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTa                                                                                                                                             >  2:2305118/1‑148 (MQ=255)
      aaTCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTAt                                                                                                                                            <  1:823243/148‑1 (MQ=255)
       aTCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATg                                                                                                                                           <  1:325811/148‑1 (MQ=255)
         ccTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGATGATGATTTTACGTGGTATCTCCAGCTATGtc                                                                                                                                         <  1:1301853/148‑1 (MQ=255)
            cAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTcc                                                                                                                                      >  2:1585994/1‑148 (MQ=255)
             aaCGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCa                                                                                                                                     >  1:1939047/1‑148 (MQ=255)
              aCGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAg                                                                                                                                    >  2:612690/1‑148 (MQ=255)
               cGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGc                                                                                                                                   >  1:1024196/1‑148 (MQ=255)
                 cagccagcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATTATGGCTTTGGTAAAACAGATCTCTCCGTGCTGGTTTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCATCTATGTCTCCAGCTa                                                                                                                                 <  2:2369255/148‑1 (MQ=255)
                 cagccagcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTa                                                                                                                                 >  1:2765618/1‑148 (MQ=255)
                 cagccagcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATAACCAGCTATGTCTCCAGCTa                                                                                                                                 <  1:1817294/148‑1 (MQ=255)
                    ccagcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTg                                                                                                                              >  2:570857/1‑148 (MQ=255)
                     cagcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGt                                                                                                                             <  2:2730270/148‑1 (MQ=255)
                      agcGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTa                                                                                                                            <  1:3061653/148‑1 (MQ=255)
                        cGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTAtc                                                                                                                          <  2:1547630/148‑1 (MQ=255)
                        cGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTAtc                                                                                                                          <  2:949825/148‑1 (MQ=255)
                        cGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTAtc                                                                                                                          >  1:1359757/1‑148 (MQ=255)
                          aTACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTAtctc                                                                                                                        <  2:802647/148‑1 (MQ=255)
                           tACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCTCCAGCTATGTCTCCAGCTACTGTATCTcc                                                                                                                       >  2:1626739/1‑148 (MQ=255)
                           tACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCTCCAGCTATGTCTCCAGCTACTGTATCTcc                                                                                                                       <  1:448744/148‑1 (MQ=255)
                           tACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTcc                                                                                                                       <  1:2149821/148‑1 (MQ=255)
                            aCCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCt                                                                                                                      <  2:114121/148‑1 (MQ=255)
                            aCCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCt                                                                                                                      <  1:916714/148‑1 (MQ=255)
                             ccTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTTCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTg                                                                                                                     >  1:753990/1‑148 (MQ=255)
                             ccTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTTCTGGGGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTg                                                                                                                     >  1:2413052/1‑148 (MQ=255)
                              cTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTgg                                                                                                                    <  1:2001176/148‑1 (MQ=255)
                               ttCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTggg                                                                                                                   >  2:101771/1‑148 (MQ=255)
                                tCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGt                                                                                                                  >  1:1746702/1‑148 (MQ=255)
                                tCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGt                                                                                                                  >  2:1296679/1‑148 (MQ=255)
                                  aTGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTAt                                                                                                                <  1:540306/148‑1 (MQ=255)
                                   tGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATc                                                                                                               >  1:1104678/1‑148 (MQ=255)
                                    gTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGCGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCa                                                                                                              <  1:3099666/148‑1 (MQ=255)
                                      tATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAgg                                                                                                            >  1:1576381/1‑148 (MQ=255)
                                       aTCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCTCCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGa                                                                                                           >  2:2006016/1‑148 (MQ=255)
                                        tCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGaa                                                                                                          >  2:7554/1‑148 (MQ=255)
                                             ccTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAggtg                                                                                                     >  2:2880872/1‑148 (MQ=255)
                                              cTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAggtgg                                                                                                    <  1:101579/148‑1 (MQ=255)
                                               ttAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAggtggt                                                                                                   >  1:3195671/1‑148 (MQ=255)
                                               ttAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAggtggt                                                                                                   >  1:794499/1‑148 (MQ=255)
                                               ttAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAggtggt                                                                                                   >  1:2325728/1‑148 (MQ=255)
                                                 aaGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTaa                                                                                                 <  2:2251853/148‑1 (MQ=255)
                                                  aGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAAt                                                                                                <  1:2725090/148‑1 (MQ=255)
                                                   gCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCTCCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATg                                                                                               <  2:280874/148‑1 (MQ=255)
                                                    ccACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGa                                                                                              <  2:865861/148‑1 (MQ=255)
                                                      aCTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGAcc                                                                                            >  2:940119/1‑148 (MQ=255)
                                                       cttcttGATGATGGATTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCa                                                                                           <  2:2432922/148‑1 (MQ=255)
                                                          cttGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGc                                                                                        >  2:2599221/1‑148 (MQ=255)
                                                                gatgGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTcgcc                                                                                  >  1:2711517/1‑148 (MQ=255)
                                                                  tgGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTcgccgc                                                                                <  1:3215235/148‑1 (MQ=255)
                                                                   ggCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCgccgcc                                                                               >  1:1091050/1‑148 (MQ=255)
                                                                   ggCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCgccgcc                                                                               >  2:1990365/1‑148 (MQ=255)
                                                                      tttGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGt                                                                            <  1:348786/148‑1 (MQ=255)
                                                                      tttGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGt                                                                            >  2:2227870/1‑148 (MQ=255)
                                                                           tAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTCATGACCATGCGTCGCCGCCTGTTTGGt                                                                       >  2:555713/1‑148 (MQ=255)
                                                                             aaaCAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTca                                                                     >  1:163595/1‑148 (MQ=255)
                                                                               aCAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTcaca                                                                   <  2:1625222/148‑1 (MQ=255)
                                                                                  gATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACAATGCGTCGCCGCCTGTTTGGTCACatga                                                                <  2:1467714/148‑1 (MQ=255)
                                                                                   aTCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACatgat                                                               <  1:291237/148‑1 (MQ=255)
                                                                                    tCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACatgatg                                                              <  2:426839/148‑1 (MQ=255)
                                                                                    tCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACatgatg                                                              <  2:217287/148‑1 (MQ=255)
                                                                                     cGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATgg                                                             <  1:1268318/148‑1 (MQ=255)
                                                                                         ccGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAAt                                                         <  2:3203670/148‑1 (MQ=255)
                                                                                          cGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATg                                                        <  2:2249015/148‑1 (MQ=255)
                                                                                            tGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGcc                                                      <  2:524105/148‑1 (MQ=255)
                                                                                              cgggTGTGGATGCCTCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAg                                                    <  2:2710845/146‑1 (MQ=255)
                                                                                              cTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAg                                                    <  2:1895646/148‑1 (MQ=255)
                                                                                              cTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAg                                                    <  1:2235098/148‑1 (MQ=255)
                                                                                                ggTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGtt                                                  >  2:1371823/1‑148 (MQ=255)
                                                                                                 gtgtGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGttt                                                 >  1:2842261/1‑148 (MQ=255)
                                                                                                  tgtgGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTc                                                <  1:1586178/148‑1 (MQ=255)
                                                                                                  tgtgGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTc                                                <  1:2474949/148‑1 (MQ=255)
                                                                                                  tgtgGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTc                                                <  1:1501959/148‑1 (MQ=255)
                                                                                                     ggATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCAtt                                             <  1:1296679/148‑1 (MQ=255)
                                                                                                     ggATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCAtt                                             <  2:2584240/148‑1 (MQ=255)
                                                                                                     ggATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCAtt                                             >  2:3069105/1‑148 (MQ=255)
                                                                                                      gATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTc                                            >  2:204074/1‑148 (MQ=255)
                                                                                                             cTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACa                                     <  1:1613329/148‑1 (MQ=255)
                                                                                                               ggtggtGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACaaa                                   <  2:2711517/148‑1 (MQ=255)
                                                                                                                gtggtgATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAAc                                  <  2:1746702/148‑1 (MQ=255)
                                                                                                                 tggtgATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACa                                 <  1:1254989/148‑1 (MQ=255)
                                                                                                                 tggtgATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACa                                 <  2:2413052/148‑1 (MQ=255)
                                                                                                                   gtgATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGt                               <  1:1594620/148‑1 (MQ=255)
                                                                                                                      aTCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCaa                            <  1:2234593/148‑1 (MQ=255)
                                                                                                                        cGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACg                          <  2:2381097/148‑1 (MQ=255)
                                                                                                                          ggCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACggg                        <  2:1527566/148‑1 (MQ=255)
                                                                                                                              gatgatTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACg                    >  2:3298751/1‑148 (MQ=255)
                                                                                                                                 gatTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCtg                 >  1:2278729/1‑148 (MQ=255)
                                                                                                                                  atTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCtgt                <  1:819587/148‑1 (MQ=255)
                                                                                                                                   ttttACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCtgtt               >  2:2552770/1‑148 (MQ=255)
                                                                                                                                      tACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTc            <  1:1154748/148‑1 (MQ=255)
                                                                                                                                       aCGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCa           >  2:401444/1‑148 (MQ=255)
                                                                                                                                           ggTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTa       <  2:279178/148‑1 (MQ=255)
                                                                                                                                            gTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTTACAAACAGTCAACGGGTACGCTGTTGTCACGTAt      <  2:3151650/148‑1 (MQ=255)
                                                                                                                                            gTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTAt      <  1:920559/148‑1 (MQ=255)
                                                                                                                                             tATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTAtt     <  2:1939047/148‑1 (MQ=255)
                                                                                                                                             tATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTAtt     >  2:439881/1‑148 (MQ=255)
                                                                                                                                              aTCTCCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTATTa    >  2:2298979/1‑148 (MQ=255)
                                                                                                                                                cACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTATTAcc  >  1:601630/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCGTTAATCCTCAACGCAGCCAGCGATACCTTCATGTTATCGCTCCTTAAGCCACTTCTTGATGATGGCTTTGGTAAAACAGATCGCTCCGTGCTGGTGTGGATGCCGCTGGTGGTGATCGGGCTGATGATTTTACGTGGTATCACCAGCTATGTCTCCAGCTACTGTATCTCCTGGGTATCAGGAAAGGTGGTAATGACCATGCGTCGCCGCCTGTTTGGTCACATGATGGGAATGCCAGTTTCATTCTTTGACAAACAGTCAACGGGTACGCTGTTGTCACGTATTACC  >  CP009273/962172‑962463

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