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,8270CA96.0% 842.4 / 21.1 285pseudogene (145/951 nt)yddKpseudogene, leucine‑rich protein; putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base A (142/130);  total (146/139)
Fisher's exact test for biased strand distribution p-value = 8.00e-02
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 285 total aligned reads displayed.

CCTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCAT  >  CP009273/1539682‑1539971
                                                                                                                                                 |                                                                                                                                                
ccTGATAAATCACTGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGGCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAt                                                                                                                                                <  2:1505480/148‑1 (MQ=255)
ccTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAt                                                                                                                                                >  2:2269/1‑148 (MQ=255)
 cTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAtt                                                                                                                                               <  1:3444866/148‑1 (MQ=255)
    aTAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggt                                                                                                                                            <  1:138929/148‑1 (MQ=255)
    aTAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggt                                                                                                                                            >  2:3381599/1‑148 (MQ=255)
      aaaTCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggttg                                                                                                                                          >  2:2171537/1‑148 (MQ=255)
       aaTCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAAttggttgg                                                                                                                                         >  2:1606492/1‑148 (MQ=255)
          cAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAg                                                                                                                                      >  2:2820456/1‑148 (MQ=255)
            aTGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTa                                                                                                                                    <  2:3074062/148‑1 (MQ=255)
             tGGTCACTTCATTAGTGTCATCAATTATATTACTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGtta                                                                                                                                   >  1:2570301/1‑146 (MQ=255)
                 cACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTg                                                                                                                               <  2:219432/148‑1 (MQ=255)
                 cACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTGGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTg                                                                                                                               <  1:3543332/148‑1 (MQ=255)
                   cTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGtt                                                                                                                             >  1:3518476/1‑148 (MQ=255)
                   cTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGtt                                                                                                                             >  1:1178893/1‑148 (MQ=255)
                   cTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGtt                                                                                                                             <  1:266507/148‑1 (MQ=255)
                    ttCATTAGTGTCATCAATTATAGTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTa                                                                                                                            <  1:658280/148‑1 (MQ=255)
                     tCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTAt                                                                                                                           <  1:2870899/148‑1 (MQ=255)
                      cATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTAtt                                                                                                                          >  1:2292385/1‑148 (MQ=255)
                          aGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCtt                                                                                                                      >  1:1207819/1‑148 (MQ=255)
                          aGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCtt                                                                                                                      >  1:1205225/1‑148 (MQ=255)
                           gtgtCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTa                                                                                                                     <  1:3111859/148‑1 (MQ=255)
                            tgtCATCAATTATATTTCTTACTGCAGTCCAGTAATGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTaa                                                                                                                    >  2:1963248/1‑148 (MQ=255)
                              tcatcaATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGt                                                                                                                  >  1:837812/1‑148 (MQ=255)
                              tcatcaATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGt                                                                                                                  >  2:2016271/1‑148 (MQ=255)
                                  caATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTaa                                                                                                              <  2:1062845/148‑1 (MQ=255)
                                    aTTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAat                                                                                                            >  2:2172580/1‑148 (MQ=255)
                                     ttATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAata                                                                                                           >  1:2095105/1‑148 (MQ=255)
                                      tataTTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAatat                                                                                                          <  2:1213470/148‑1 (MQ=255)
                                       atatTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATAtt                                                                                                         <  1:768158/148‑1 (MQ=255)
                                        tatTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTc                                                                                                        >  2:1023786/1‑148 (MQ=255)
                                         atTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTcc                                                                                                       >  2:923642/1‑148 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                     <  1:422038/148‑1 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                     >  1:3300313/1‑148 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                     <  1:1049897/148‑1 (MQ=255)
                                           ttCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAg                                                                                                     <  1:1793770/148‑1 (MQ=255)
                                              ttACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGttt                                                                                                  >  2:1566878/1‑148 (MQ=255)
                                               tACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGtttt                                                                                                 >  2:2124928/1‑148 (MQ=255)
                                               tACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGtttt                                                                                                 <  2:388636/148‑1 (MQ=255)
                                                   gCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAg                                                                                             <  1:878772/148‑1 (MQ=255)
                                                     aGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAgtg                                                                                           <  1:2287805/148‑1 (MQ=255)
                                                      gTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAgtgt                                                                                          >  1:1480185/1‑148 (MQ=255)
                                                       tCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGtt                                                                                         >  2:2907737/1‑148 (MQ=255)
                                                         cAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGtttt                                                                                       <  2:496822/148‑1 (MQ=255)
                                                         cAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGtttt                                                                                       <  2:1452316/148‑1 (MQ=255)
                                                          aGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGttttt                                                                                      <  1:3189050/148‑1 (MQ=255)
                                                          aGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGttttt                                                                                      <  1:1508663/148‑1 (MQ=255)
                                                           gTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTg                                                                                     >  1:762550/1‑148 (MQ=255)
                                                           gTACTGAATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTg                                                                                     <  2:3486674/148‑1 (MQ=255)
                                                            tACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTgg                                                                                    >  2:1179743/1‑148 (MQ=255)
                                                            tACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTgg                                                                                    >  2:1269575/1‑148 (MQ=255)
                                                               tGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGt                                                                                 <  2:1519110/148‑1 (MQ=255)
                                                                   tttGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTaa                                                                             <  2:2140289/148‑1 (MQ=255)
                                                                    ttGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTaaa                                                                            <  2:2646918/148‑1 (MQ=255)
                                                                      gtgCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAAtg                                                                          <  1:2062445/148‑1 (MQ=255)
                                                                          aaGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCa                                                                      <  2:577806/148‑1 (MQ=255)
                                                                          aaGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCa                                                                      <  2:3459459/148‑1 (MQ=255)
                                                                             tGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAata                                                                   >  1:938748/1‑148 (MQ=255)
                                                                             tGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAata                                                                   >  1:2720690/1‑148 (MQ=255)
                                                                              ggATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAatat                                                                  <  2:1701824/148‑1 (MQ=255)
                                                                              ggATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAatat                                                                  <  2:2564341/148‑1 (MQ=255)
                                                                              ggATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAatat                                                                  <  1:983091/148‑1 (MQ=255)
                                                                                   aTTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTg                                                             <  1:2402119/148‑1 (MQ=255)
                                                                                   aTTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTg                                                             <  1:1023521/148‑1 (MQ=255)
                                                                                    ttgttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGt                                                            <  1:3381599/148‑1 (MQ=255)
                                                                                      gttgATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCg                                                          <  1:3144283/148‑1 (MQ=255)
                                                                                                ttAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGa                                                <  1:2921272/148‑1 (MQ=255)
                                                                                                ttAACTAATACATGCCAAAGATGTTGTGATGATATTAGTTGATCAATGTAATTGGTAGGCAGTAAATTGTTATTGCTAAAGTTTAAATAATCAAGTTTTCTAGTGTTTTTGGCGTTTAAATGTTCAATATGGTTGTCGTTGAAATTGa                                                >  1:2460333/1‑148 (MQ=255)
                                                                                                    cTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGt                                            <  1:1323828/148‑1 (MQ=255)
                                                                                                    cTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGt                                            >  1:2892977/1‑148 (MQ=255)
                                                                                                      aaTACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGttt                                          <  2:435627/148‑1 (MQ=255)
                                                                                                         aCATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCa                                       <  2:1648614/148‑1 (MQ=255)
                                                                                                          cATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCAt                                      <  2:1394213/148‑1 (MQ=255)
                                                                                                          cATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCAt                                      <  2:240512/148‑1 (MQ=255)
                                                                                                           aTGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATc                                     <  1:2820462/148‑1 (MQ=255)
                                                                                                           aTGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATc                                     <  1:1630965/148‑1 (MQ=255)
                                                                                                              ccAAAGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTg                                  >  1:869078/1‑148 (MQ=255)
                                                                                                                aaaGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGAGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTggg                                <  1:559680/148‑1 (MQ=255)
                                                                                                                aaaGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGAGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTggg                                <  1:1004333/148‑1 (MQ=255)
                                                                                                                aaaGATGTTTTGATGATATGAGTTGATCAATGTCATTGGTTGTCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTggg                                <  1:113459/148‑1 (MQ=255)
                                                                                                                 aaGATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGga                               >  1:2362760/1‑148 (MQ=255)
                                                                                                                   gATGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgatg                             >  2:1684103/1‑148 (MQ=255)
                                                                                                                    aTGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgatga                            <  1:1037882/148‑1 (MQ=255)
                                                                                                                    aTGTTTTGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGgatga                            >  2:268648/1‑148 (MQ=255)
                                                                                                                       ttttGATGATATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTTTTTTAATCCTGGGATGATtg                         >  1:1599576/1‑148 (MQ=255)
                                                                                                                          tgatgaTATTAGTTGATCAATTTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTATAAATATTCCAGATTATTAGTGTTTTTGTCGTTTAAATGTTCAATATGGTTCTCGTTTAAAATGATTTTTTTCATCATTGGATTATTatgt                      >  2:3087501/1‑148 (MQ=255)
                                                                                                                          tgatgaTATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTgtgt                      <  2:1474394/148‑1 (MQ=255)
                                                                                                                           gatgatATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTa                     <  2:1205225/148‑1 (MQ=255)
                                                                                                                            atgatATTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTaa                    >  2:672120/1‑148 (MQ=255)
                                                                                                                               atatTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAatat                 <  2:3219775/148‑1 (MQ=255)
                                                                                                                                 atTAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATaa               >  1:2181683/1‑148 (MQ=255)
                                                                                                                                   tAGTTGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAgg             <  1:3286568/148‑1 (MQ=255)
                                                                                                                                      ttGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCa          >  1:1753474/1‑148 (MQ=255)
                                                                                                                                       tGATCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAg         <  2:3236241/148‑1 (MQ=255)
                                                                                                                                         aTCAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTg       <  2:3300196/148‑1 (MQ=255)
                                                                                                                                           cAATGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGAt     <  2:83318/148‑1 (MQ=255)
                                                                                                                                            aaTGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATc    <  2:1311037/148‑1 (MQ=255)
                                                                                                                                            aaTGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATc    >  2:229081/1‑148 (MQ=255)
                                                                                                                                             aTGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCa   <  2:535276/148‑1 (MQ=255)
                                                                                                                                              tGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCat  <  1:1963248/148‑1 (MQ=255)
                                                                                                                                              tGTAATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCat  >  2:1807493/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCTGATAAATCAATGGTCACTTCATTAGTGTCATCAATTATATTTCTTACTGCAGTCCAGTACTGCATTTGTGCAAGTGGATCATTGTTGATGCCGTTAACTAATACATGCCAAAGATGTTTTGATGATATTAGTTGATCAATGTCATTGGTTGGCAGTAAATTGTTATTGCTTAAGTTTAAATATTCCAGTTTTGTAGTGTTTTTGGCGTTTAAATGTGCAATATGGTTGTCGTTTAAAGTGATTGTTTTCATCCTGGGATGATTGTGTAATATAAGGTCAGTGATCAT  >  CP009273/1539682‑1539971

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: