Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,539,827 C→A 100% pseudogene (145/951 nt) yddK ← pseudogene, leucine‑rich protein; putative glycoportein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,539,8270CA98.3% 982.0 / 5.1 308pseudogene (145/951 nt)yddKpseudogene, leucine‑rich protein; putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (134/167);  total (138/170)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 308 total aligned reads displayed.

CCTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCA  >  CP009273/1539682‑1539970
                                                                                                                                                 |                                                                                                                                               
 cTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAtt                                                                                                                                              <  2:920734/148‑1 (MQ=255)
   gATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttgg                                                                                                                                            <  2:1370759/148‑1 (MQ=255)
    aTAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggt                                                                                                                                           <  1:2855572/148‑1 (MQ=255)
    aTAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggt                                                                                                                                           <  1:2018048/148‑1 (MQ=255)
        aTCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGc                                                                                                                                       >  2:61219/1‑148 (MQ=255)
         tCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCa                                                                                                                                      <  1:1164579/148‑1 (MQ=255)
         tCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAGTGGTTGGCa                                                                                                                                      <  1:2366739/148‑1 (MQ=255)
          cAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAg                                                                                                                                     <  2:2014087/148‑1 (MQ=255)
             tGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTaa                                                                                                                                  >  2:1614842/1‑148 (MQ=255)
              ggTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTaaa                                                                                                                                 <  1:347924/148‑1 (MQ=255)
               gTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAAt                                                                                                                                >  1:1027675/1‑148 (MQ=255)
               gTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAAt                                                                                                                                <  2:3543190/148‑1 (MQ=255)
               gTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAAt                                                                                                                                <  2:101059/148‑1 (MQ=255)
                   cTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGtt                                                                                                                            <  2:78466/148‑1 (MQ=255)
                    ttCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTa                                                                                                                           <  1:62346/148‑1 (MQ=255)
                     tCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTAt                                                                                                                          >  1:3473537/1‑148 (MQ=255)
                         tAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCt                                                                                                                      <  1:1831794/148‑1 (MQ=255)
                         tAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCt                                                                                                                      <  1:1595382/148‑1 (MQ=255)
                           gtgtCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTa                                                                                                                    <  2:402179/148‑1 (MQ=255)
                             gtCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAg                                                                                                                  <  1:1454472/148‑1 (MQ=255)
                               catcaATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGtt                                                                                                                >  2:1821551/1‑148 (MQ=255)
                               catcaATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGtt                                                                                                                >  2:2736571/1‑148 (MQ=255)
                                 tcaATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTa                                                                                                              >  1:1929117/1‑148 (MQ=255)
                                  caATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTaa                                                                                                             >  1:2542381/1‑148 (MQ=255)
                                  caATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTaa                                                                                                             >  2:1059442/1‑148 (MQ=255)
                                   aaTTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTaaa                                                                                                            >  2:3424688/1‑148 (MQ=255)
                                      tataTTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAatat                                                                                                         >  2:557542/1‑148 (MQ=255)
                                        tatTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTc                                                                                                       >  1:2241156/1‑148 (MQ=255)
                                         atTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTcc                                                                                                      >  2:1758263/1‑148 (MQ=255)
                                         atTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTcc                                                                                                      <  1:1608645/148‑1 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                    <  2:2128951/148‑1 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                    <  2:336375/148‑1 (MQ=255)
                                               tACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGtttt                                                                                                >  2:2625065/1‑148 (MQ=255)
                                                aCTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTg                                                                                               >  2:1414035/1‑148 (MQ=255)
                                                aCTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTg                                                                                               <  2:2598212/148‑1 (MQ=255)
                                                  tGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTa                                                                                             <  2:2099950/148‑1 (MQ=255)
                                                  tGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTa                                                                                             >  2:2882057/1‑148 (MQ=255)
                                                      gTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAgtgt                                                                                         <  1:658540/148‑1 (MQ=255)
                                                          aGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGttttt                                                                                     >  1:1638670/1‑148 (MQ=255)
                                                            tACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTgg                                                                                   >  1:2467874/1‑148 (MQ=255)
                                                             aCTGCATTTGTGCAAGTGGATCATTGTTTATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTAAAATATTCCAGTTTTGTATTGtttttttc                                                                                  >  2:3034587/1‑145 (MQ=255)
                                                             aCTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGc                                                                                  >  2:681497/1‑148 (MQ=255)
                                                             aCTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGc                                                                                  >  2:2466543/1‑148 (MQ=255)
                                                             aCTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGc                                                                                  >  1:3083143/1‑148 (MQ=255)
                                                                gCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGGTGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGtt                                                                               >  2:2585665/1‑148 (MQ=255)
                                                                gCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGtt                                                                               <  2:2416120/148‑1 (MQ=255)
                                                                   tttGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTaa                                                                            <  2:2664723/148‑1 (MQ=255)
                                                                     tgtgCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAAt                                                                          <  2:3324269/148‑1 (MQ=255)
                                                                        gCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTTCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAAtgtg                                                                       <  2:1828816/148‑1 (MQ=255)
                                                                        gCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAAtgtg                                                                       >  1:2695611/1‑148 (MQ=255)
                                                                            gTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAat                                                                   >  1:3399311/1‑148 (MQ=255)
                                                                            gTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAat                                                                   <  1:2669678/148‑1 (MQ=255)
                                                                             tGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAata                                                                  >  2:1142078/1‑148 (MQ=255)
                                                                             tGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAata                                                                  <  2:317400/148‑1 (MQ=255)
                                                                              ggATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGTCGTTTAAATGTGCAatat                                                                 >  2:1692553/1‑148 (MQ=255)
                                                                                 tCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGt                                                              >  1:535308/1‑148 (MQ=255)
                                                                                    ttgttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGt                                                           >  1:1087620/1‑148 (MQ=255)
                                                                                      gttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCg                                                         >  1:949667/1‑148 (MQ=255)
                                                                                      gttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCg                                                         <  2:1173599/148‑1 (MQ=255)
                                                                                       ttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGt                                                        <  2:70108/148‑1 (MQ=255)
                                                                                           tGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTaa                                                    <  1:1017184/148‑1 (MQ=255)
                                                                                           tGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTaa                                                    <  1:823233/148‑1 (MQ=255)
                                                                                           tGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTaa                                                    >  2:1428580/1‑148 (MQ=255)
                                                                                            gCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTaaa                                                   <  2:2002807/148‑1 (MQ=255)
                                                                                            gCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTaaa                                                   >  2:3115335/1‑148 (MQ=255)
                                                                                            gCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTGTGGCGTTTAAATGTGCAATATGGTTGGCGGGTaaa                                                   >  2:3285293/1‑148 (MQ=255)
                                                                                                 tAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGAt                                              >  2:1232863/1‑148 (MQ=255)
                                                                                                  aaCTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGAtt                                             <  2:1562694/148‑1 (MQ=255)
                                                                                                   aCTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTg                                            >  2:2533260/1‑148 (MQ=255)
                                                                                                    cTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGt                                           <  2:3553802/148‑1 (MQ=255)
                                                                                                     tAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGtt                                          <  1:1758263/148‑1 (MQ=255)
                                                                                                     tAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGtt                                          <  1:157219/148‑1 (MQ=255)
                                                                                                      aaTACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGttt                                         <  1:362004/148‑1 (MQ=255)
                                                                                                        tACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTc                                       <  2:2851745/148‑1 (MQ=255)
                                                                                                        tACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTc                                       <  2:2913058/148‑1 (MQ=255)
                                                                                                         aCATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCa                                      <  2:652249/148‑1 (MQ=255)
                                                                                                         aCATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCa                                      <  2:1929117/148‑1 (MQ=255)
                                                                                                          cATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCAt                                     <  2:2690841/148‑1 (MQ=255)
                                                                                                            tGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTTTTTTCATcc                                   >  1:606882/1‑148 (MQ=255)
                                                                                                                  aGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgat                             <  1:2093681/148‑1 (MQ=255)
                                                                                                                   gATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgatg                            >  1:1353738/1‑148 (MQ=255)
                                                                                                                    aTGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgatga                           >  2:1530498/1‑148 (MQ=255)
                                                                                                                        tttGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATtgt                       >  1:1499691/1‑148 (MQ=255)
                                                                                                                           gatgatATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTAATCCTGGGATGATTGTGTa                    >  1:2779654/1‑148 (MQ=255)
                                                                                                                             tgatATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAat                  <  1:3292624/148‑1 (MQ=255)
                                                                                                                              gatATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAata                 >  2:748591/1‑148 (MQ=255)
                                                                                                                               atatTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAatat                <  2:2305806/148‑1 (MQ=255)
                                                                                                                                tatTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAAtata               >  2:726146/1‑148 (MQ=255)
                                                                                                                                tatTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAAtata               <  1:3136181/148‑1 (MQ=255)
                                                                                                                                 agtAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTGGTCGTTTAAAGTTATTGTTTTCATAATGGGATGATTGTGTTATATaa              >  1:2710679/3‑148 (MQ=255)
                                                                                                                                   tAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGATTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAgg            >  2:501193/1‑148 (MQ=255)
                                                                                                                                    aGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGt           <  2:1033575/148‑1 (MQ=255)
                                                                                                                                    aGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGt           <  1:2332149/148‑1 (MQ=255)
                                                                                                                                    aGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGt           <  1:1821551/148‑1 (MQ=255)
                                                                                                                                     gTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTc          >  2:2522219/1‑148 (MQ=255)
                                                                                                                                      ttGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTTTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCa         <  1:1851259/148‑1 (MQ=255)
                                                                                                                                      ttGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCa         >  2:823436/1‑148 (MQ=255)
                                                                                                                                           cAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGAt    <  1:2211703/148‑1 (MQ=255)
                                                                                                                                             aTGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCa         >  2:1461699/1‑141 (MQ=255)
                                                                                                                                             aTGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCa  <  1:1541010/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCA  >  CP009273/1539682‑1539970

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: