Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,785,915 G→T 85.3% intergenic (‑4/+271) kbl ← / ← yibB glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,785,9150GT85.3% 591.4 / 106.1 232intergenic (‑4/+271)kbl/yibBglycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (13/21);  new base T (89/109);  total (102/130)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 233 total aligned reads displayed.

ACAAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGT  >  CP009273/3785768‑3786064
                                                                                                                                                   |                                                                                                                                                     
acaAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTAt                                                                                                                                                     >  2:1022260/1‑150 (MQ=255)
 caAAAGTTAATGACGTGGTTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTAtt                                                                                                                                                    <  1:784836/150‑1 (MQ=255)
  aaaaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATtc                                                                                                                                                   <  1:3584536/150‑1 (MQ=255)
  aaaaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATtc                                                                                                                                                   <  1:3584600/150‑1 (MQ=255)
     aGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTcc                                                                                                                                                <  1:652746/150‑1 (MQ=255)
     aGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTcc                                                                                                                                                <  2:3937009/150‑1 (MQ=255)
     aGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTcc                                                                                                                                                <  2:437267/150‑1 (MQ=255)
     aGTTAAGGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTcc                                                                                                                                                <  2:1655289/150‑1 (MQ=255)
      gTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCa                                                                                                                                               >  2:1157055/1‑150 (MQ=255)
      gTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCa                                                                                                                                               >  2:2027105/1‑150 (MQ=255)
      gTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCa                                                                                                                                               >  2:3828223/1‑150 (MQ=255)
         aaTGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGAc                                                                                                                                            >  2:3193590/1‑150 (MQ=255)
             aCGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAggg                                                                                                                                        <  1:4260815/150‑1 (MQ=255)
             aCGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGAGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCATATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAggg                                                                                                                                        >  1:9652/1‑150 (MQ=255)
               gTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCa                                                                                                                                      <  2:945614/150‑1 (MQ=255)
                    ttCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTc                                                                                                                                 >  1:959357/1‑150 (MQ=255)
                    ttCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTc                                                                                                                                 <  1:455372/150‑1 (MQ=255)
                         tCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGca                                                                                                                            >  1:538529/1‑150 (MQ=255)
                         tCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGca                                                                                                                            >  2:752272/1‑150 (MQ=255)
                              cacaGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACatat                                                                                                                       <  2:2210980/150‑1 (MQ=255)
                                caGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTa                                                                                                                     <  1:996974/150‑1 (MQ=255)
                                caGTGATATCTGCTTGCTGCGCAGACGTAATAAGGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTa                                                                                                                     <  1:1000326/150‑1 (MQ=255)
                                       atCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAAc                                                                                                              <  2:207518/150‑1 (MQ=255)
                                         cTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACtt                                                                                                            <  2:1109745/150‑1 (MQ=255)
                                         cTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACtt                                                                                                            <  1:3147268/150‑1 (MQ=255)
                                          tgcttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTta                                                                                                           >  1:3257160/1‑150 (MQ=255)
                                            cttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTtata                                                                                                         >  1:2894419/1‑150 (MQ=255)
                                             ttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAg                                                                                                        >  1:918885/1‑150 (MQ=255)
                                              tgctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAgg                                                                                                       >  1:165120/1‑150 (MQ=255)
                                              tgctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAgg                                                                                                       <  1:2485549/150‑1 (MQ=255)
                                              tgctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAgg                                                                                                       <  1:1173652/150‑1 (MQ=255)
                                                ctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTGACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTc                                                                                                     <  2:1574297/150‑1 (MQ=255)
                                                ctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACTGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTc                                                                                                     <  2:1459690/150‑1 (MQ=255)
                                                  gcgcAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTcgc                                                                                                   <  1:314072/150‑1 (MQ=255)
                                                   cgcAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTcgcg                                                                                                  >  2:4476670/1‑150 (MQ=255)
                                                   cgcAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTcgcg                                                                                                  <  1:3790940/150‑1 (MQ=255)
                                                   cgcAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTcgcg                                                                                                  <  1:3791045/150‑1 (MQ=255)
                                                     cAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGcgac                                                                                                <  2:347739/150‑1 (MQ=255)
                                                        aCGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACgag                                                                                             >  1:3040236/1‑150 (MQ=255)
                                                         cGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACgaga                                                                                            <  2:1694086/150‑1 (MQ=255)
                                                         cGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACgaga                                                                                            <  1:4344477/150‑1 (MQ=255)
                                                            aataatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATaa                                                                                         >  1:3852835/1‑150 (MQ=255)
                                                              taatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcg                                                                                       <  2:491029/150‑1 (MQ=255)
                                                              taatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcg                                                                                       <  1:2645038/150‑1 (MQ=255)
                                                                   cgcTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTa                                                                                  >  2:269843/1‑150 (MQ=255)
                                                                      tcttctTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAaca                                                                               <  2:528064/150‑1 (MQ=255)
                                                                       cttcatTAAACACCCCGTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAacac                                                                              <  2:2422523/150‑1 (MQ=255)
                                                                         tcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACtt                                                                            <  2:1409484/150‑1 (MQ=255)
                                                                           tttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCt                                                                          <  1:1597006/150‑1 (MQ=255)
                                                                           tttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCt                                                                          <  1:2458008/150‑1 (MQ=255)
                                                                            ttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTg                                                                         >  1:878565/1‑150 (MQ=255)
                                                                             tAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGc                                                                        >  2:1431727/1‑150 (MQ=255)
                                                                               aacaacCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGTTCGCGACGAGATAACGCGTTAACACTTCTGCaa                                                                      >  2:1856311/1‑150 (MQ=255)
                                                                                  aacCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAt                                                                   <  2:3012815/150‑1 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                                  >  2:890039/1‑150 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                                  <  2:3357398/150‑1 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                                  <  1:2183735/150‑1 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                                  <  1:4397899/150‑1 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                                  <  2:304418/150‑1 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtc                                                                  <  1:2183563/150‑2 (MQ=255)
                                                                                     ccTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCa                                                                <  2:3558217/150‑1 (MQ=255)
                                                                                           gCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataa                                                          >  1:393100/1‑150 (MQ=255)
                                                                                            cccGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaa                                                         <  1:3232583/150‑1 (MQ=255)
                                                                                            cccGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaa                                                         <  1:2513457/150‑1 (MQ=255)
                                                                                            cccGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaa                                                         <  1:3232520/150‑1 (MQ=255)
                                                                                            cccGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaa                                                         >  2:2256794/1‑150 (MQ=255)
                                                                                              cGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaatc                                                       >  1:394155/1‑149 (MQ=255)
                                                                                              cGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaata                                                       <  1:610995/150‑1 (MQ=255)
                                                                                              cGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaata                                                       >  1:394021/1‑150 (MQ=255)
                                                                                               gTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataaataa                                                      >  2:3777919/1‑150 (MQ=255)
                                                                                                tGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAAc                                                     <  1:242364/150‑1 (MQ=255)
                                                                                                   ggTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtg                                                  >  1:1838369/1‑150 (MQ=255)
                                                                                                    gTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtgt                                                 >  2:3945687/1‑150 (MQ=255)
                                                                                                        ccAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACt                                             <  1:10038/150‑1 (MQ=255)
                                                                                                            aTCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGAt                                         <  1:2849549/150‑1 (MQ=255)
                                                                                                              cGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGt                                       >  2:2773625/1‑150 (MQ=255)
                                                                                                                ttGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtct                                     <  1:1378527/150‑1 (MQ=255)
                                                                                                                 tGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctc                                    >  2:4159624/1‑150 (MQ=255)
                                                                                                                 tGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctc                                    >  1:4008104/1‑150 (MQ=255)
                                                                                                                 tGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAAGATAAATAACGTGTAACTTGATGtctc                                    >  2:4045301/1‑150 (MQ=255)
                                                                                                                  ggTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCTGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctct                                   >  1:3157257/1‑150 (MQ=255)
                                                                                                                  ggTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctct                                   >  1:3036076/1‑150 (MQ=255)
                                                                                                                   gTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctctc                                  >  1:1786512/1‑150 (MQ=255)
                                                                                                                   gTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctctc                                  <  1:739275/150‑1 (MQ=255)
                                                                                                                     tAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCaa                                >  2:2978209/1‑150 (MQ=255)
                                                                                                                     tAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCAAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCaa                                >  2:4356264/1‑150 (MQ=255)
                                                                                                                           ctgATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTa                          <  1:3895532/150‑1 (MQ=255)
                                                                                                                             gATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACa                        >  2:492939/1‑150 (MQ=255)
                                                                                                                             gATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGAGAAATAACGTGTAACTTGATGTCTCTCAACACCTACa                        >  1:4460149/1‑150 (MQ=255)
                                                                                                                                aaaaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAt                     >  1:3364676/1‑150 (MQ=255)
                                                                                                                                aaaaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAt                     <  1:2579825/150‑1 (MQ=255)
                                                                                                                                  aaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtag                   <  2:1281546/150‑1 (MQ=255)
                                                                                                                                  aaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtag                   <  1:3827850/150‑1 (MQ=255)
                                                                                                                                  aaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCAGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtag                   <  2:1281782/150‑1 (MQ=255)
                                                                                                                                    ttCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtagta                 >  1:1989200/1‑150 (MQ=255)
                                                                                                                                      ctcCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGa               <  1:492939/150‑1 (MQ=255)
                                                                                                                                         cACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGaaaa            >  1:2724066/1‑150 (MQ=255)
                                                                                                                                            gCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACt         <  1:1022260/150‑1 (MQ=255)
                                                                                                                                              aTTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGt       >  1:1776341/1‑150 (MQ=255)
                                                                                                                                                tGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTAc     <  1:3089936/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGT  >  CP009273/3785768‑3786064

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: