Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,520,789 G→T 17.1% R321L (CGA→CTA)  yrfF → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,520,7890GT17.1% 738.3 / 110.3 251R321L (CGA→CTA) yrfFinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (113/95);  new base T (14/29);  total (127/124)
Fisher's exact test for biased strand distribution p-value = 1.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 255 total aligned reads displayed.

GAAAACGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGA  >  CP009273/3520641‑3520938
                                                                                                                                                    |                                                                                                                                                     
gAAAACGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGa                                                                                                                                                      >  2:4281003/1‑150 (MQ=255)
 aaaaCGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGat                                                                                                                                                     >  1:2553331/1‑150 (MQ=255)
   aaCGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGatat                                                                                                                                                   >  1:3195822/1‑150 (MQ=255)
   aaCGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGatat                                                                                                                                                   >  2:1293001/1‑150 (MQ=255)
      gATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCtt                                                                                                                                                <  1:2742211/150‑1 (MQ=255)
       aTCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCttt                                                                                                                                               <  1:1833274/150‑1 (MQ=255)
       aTCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCttt                                                                                                                                               >  1:2861473/1‑150 (MQ=255)
        tCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTc                                                                                                                                              >  1:1746770/1‑150 (MQ=255)
         cAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCg                                                                                                                                             <  1:4536122/150‑1 (MQ=255)
         cAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTca                                                                                                                                             <  2:689435/150‑2 (MQ=255)
         cAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCg                                                                                                                                             >  2:2151470/1‑150 (MQ=255)
         cAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCg                                                                                                                                             >  1:1739505/1‑150 (MQ=255)
           ggAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCt                                                                                                                                           <  2:2705006/150‑1 (MQ=255)
                  aTCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATCTCTTTCGCTGCatgat                                                                                                                                    >  2:2695792/1‑150 (MQ=255)
                  aTCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCatgat                                                                                                                                    <  1:675735/150‑1 (MQ=255)
                   tCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCatgatg                                                                                                                                   >  1:3342981/1‑150 (MQ=255)
                   tCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCatgatg                                                                                                                                   <  1:2260185/150‑1 (MQ=255)
                      acaaTATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAg                                                                                                                                >  1:4828772/1‑150 (MQ=255)
                      acaaTATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAg                                                                                                                                <  1:2571442/150‑1 (MQ=255)
                             ttCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTaaaaaa                                                                                                                         <  2:3542943/150‑1 (MQ=255)
                             ttCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTaaaaaa                                                                                                                         >  2:4937723/1‑150 (MQ=255)
                               cACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACt                                                                                                                       <  2:2506784/150‑1 (MQ=255)
                               cACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACt                                                                                                                       <  1:4569732/150‑1 (MQ=255)
                               cACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACt                                                                                                                       >  1:3442559/1‑150 (MQ=255)
                                aCTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAATACtt                                                                                                                      >  2:3751959/1‑150 (MQ=255)
                                aCTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACtt                                                                                                                      >  1:2389503/1‑150 (MQ=255)
                                  tCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTc                                                                                                                    <  1:5037940/150‑1 (MQ=255)
                                   cGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTcc                                                                                                                   <  2:4697980/150‑1 (MQ=255)
                                   cGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTcc                                                                                                                   <  2:2112442/150‑1 (MQ=255)
                                    ggTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCg                                                                                                                  >  1:1255955/1‑150 (MQ=255)
                                      tattatCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGtt                                                                                                                <  2:3011515/150‑1 (MQ=255)
                                         tatCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATCTCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACa                                                                                                             >  2:3301642/1‑150 (MQ=255)
                                         tatCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCTACCGCCAGGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGGTACa                                                                                                             >  2:4833665/1‑150 (MQ=255)
                                         tatCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACa                                                                                                             >  2:2760248/1‑150 (MQ=255)
                                         gatCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACa                                                                                                             <  2:4499580/149‑1 (MQ=255)
                                          atCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAg                                                                                                            <  2:3651734/150‑1 (MQ=255)
                                          atCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAg                                                                                                            >  1:1531221/1‑150 (MQ=255)
                                          atCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAg                                                                                                            >  1:4011107/1‑150 (MQ=255)
                                                  ggTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCt                                                                                                    <  2:4113116/150‑1 (MQ=255)
                                                   gTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCTACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTg                                                                                                   <  2:4334852/150‑1 (MQ=255)
                                                    tCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCTTGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTgc                                                                                                  <  1:3914511/150‑1 (MQ=255)
                                                     cTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTgcg                                                                                                 <  2:312066/150‑1 (MQ=255)
                                                      tATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTgcgc                                                                                                >  2:2546094/1‑150 (MQ=255)
                                                      tATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTgcgc                                                                                                >  2:1196800/1‑150 (MQ=255)
                                                       aTCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCa                                                                                               >  2:349923/1‑150 (MQ=255)
                                                             cacaCTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGa                                                                                         >  1:1390683/1‑150 (MQ=255)
                                                               cacTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGAtt                                                                                       <  2:4178659/150‑1 (MQ=255)
                                                               cacTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGAtt                                                                                       >  1:374738/1‑150 (MQ=255)
                                                               cacTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGAtt                                                                                       <  2:4434630/150‑1 (MQ=255)
                                                                 cTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTAt                                                                                     >  1:288138/1‑150 (MQ=255)
                                                                 cTGGCAGCCATACATGGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTAt                                                                                     <  1:5126298/150‑1 (MQ=255)
                                                                   ggCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCg                                                                                   <  2:4471028/150‑1 (MQ=255)
                                                                        ccATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCg                                                                              >  2:359177/1‑150 (MQ=255)
                                                                         cATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCgg                                                                             <  2:1712734/150‑1 (MQ=255)
                                                                         cATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCgg                                                                             <  2:3440923/150‑1 (MQ=255)
                                                                           tACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGc                                                                           >  1:4306573/1‑150 (MQ=255)
                                                                            aCATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCt                                                                          >  1:3675800/1‑150 (MQ=255)
                                                                             cATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTc                                                                         <  1:4428651/150‑1 (MQ=255)
                                                                              aTTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCg                                                                        <  2:51750/150‑1 (MQ=255)
                                                                                tGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgct                                                                      <  2:3792175/150‑1 (MQ=255)
                                                                                 gCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgctg                                                                     <  2:4764248/150‑1 (MQ=255)
                                                                                 gCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgctg                                                                     <  1:564838/150‑1 (MQ=255)
                                                                                  cTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgctgc                                                                    >  1:1455755/1‑150 (MQ=255)
                                                                                  cTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgctgc                                                                    >  1:1455815/1‑150 (MQ=255)
                                                                                   tCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCgctgct                                                                   >  1:4840949/1‑150 (MQ=255)
                                                                                    cAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGctgctg                                                                  <  1:4214690/150‑1 (MQ=255)
                                                                                    cAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGctgctg                                                                  >  1:2391920/1‑150 (MQ=255)
                                                                                     aGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTgg                                                                 >  1:3111675/1‑150 (MQ=255)
                                                                                       gATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTg                                                               >  2:3934785/1‑150 (MQ=255)
                                                                                       gATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTg                                                               >  1:2733645/1‑150 (MQ=255)
                                                                                        aTCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGc                                                              >  2:38832/1‑150 (MQ=255)
                                                                                          ctcGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTg                                                            <  2:3819553/150‑1 (MQ=255)
                                                                                            cGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGtt                                                          <  1:208655/150‑1 (MQ=255)
                                                                                            cGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGtt                                                          >  1:2915862/1‑150 (MQ=255)
                                                                                              gTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTa                                                        >  2:1053600/1‑150 (MQ=255)
                                                                                              gTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTa                                                        >  2:4581435/1‑150 (MQ=255)
                                                                                               tCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTAt                                                       <  2:250864/150‑1 (MQ=255)
                                                                                                   caaACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTg                                                   <  1:5185640/150‑1 (MQ=255)
                                                                                                   caaACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTg                                                   >  1:3655987/1‑150 (MQ=255)
                                                                                                      aCCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTa                                                >  1:2979068/1‑150 (MQ=255)
                                                                                                       ccGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTAt                                               <  2:2421832/150‑1 (MQ=255)
                                                                                                         gatatcgatatcTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTc                                             <  1:3024805/150‑1 (MQ=255)
                                                                                                          atatcgatatcTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCt                                            <  1:1001246/150‑1 (MQ=255)
                                                                                                           tatcgatatcTATCTCGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTg                                           <  2:3036318/150‑1 (MQ=255)
                                                                                                            atcgatatcTATCTCTACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTgg                                          <  2:18851/150‑1 (MQ=255)
                                                                                                               gatatcTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATc                                       >  1:3503228/1‑150 (MQ=255)
                                                                                                                    ctatctCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCt                                  <  2:2553384/150‑1 (MQ=255)
                                                                                                                      atctCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTgg                                <  2:1461361/150‑1 (MQ=255)
                                                                                                                      atctCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTgg                                >  1:498694/1‑150 (MQ=255)
                                                                                                                       tctcGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGa                               <  1:232278/150‑1 (MQ=255)
                                                                                                                         tcGACCGCCATGTAGTGCGTCAGGGACTATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGata                             >  2:3996061/1‑150 (MQ=255)
                                                                                                                             ccGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGcc                         >  1:2662569/1‑150 (MQ=255)
                                                                                                                                  aTGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGa                    >  1:1022088/1‑150 (MQ=255)
                                                                                                                                     tAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAAt                 >  2:576663/1‑150 (MQ=255)
                                                                                                                                            tCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCActct          >  1:864690/1‑150 (MQ=255)
                                                                                                                                              aGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCActctct        >  2:359672/1‑150 (MQ=255)
                                                                                                                                                 gACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCAt     >  1:2063607/1‑150 (MQ=255)
                                                                                                                                                    gATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGa  >  2:3052601/1‑150 (MQ=255)
                                                                                                                                                    gATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGa  >  2:5095250/1‑150 (MQ=255)
                                                                                                                                                    gATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGa  <  1:163890/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAAACGATCAGGAACAGATCAACAATATTTCACTCGGTATTATCGACCTGGTCTATCCCGCACACTGGCAGCCATACATTGCTCAGGATCTCGGTCAACAAACCGATATCGATATCTATCTCGACCGCCATGTAGTGCGTCAGGGACGATATCTTTCGCTGCATGATGAAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGA  >  CP009273/3520641‑3520938

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: