Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,181,300 C→T 5.7% pseudogene (521/2445 nt) yqiG → pseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,181,3000CT 5.7% 1205.3 / 38.4 332pseudogene (521/2445 nt)yqiGpseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (161/152);  new base T (9/10);  total (170/162)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 337 total aligned reads displayed.

GCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGC  >  CP009273/3181152‑3181447
                                                                                                                                                    |                                                                                                                                                   
gCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTc                                                                                                                                                     >  1:4134802/1‑149 (MQ=255)
  aCCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCt                                                                                                                                                   <  1:1359935/149‑1 (MQ=255)
  aCCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCt                                                                                                                                                   <  1:4912407/149‑1 (MQ=255)
  aCCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCt                                                                                                                                                   <  2:3786576/149‑1 (MQ=255)
   ccTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTa                                                                                                                                                  >  1:4244143/1‑149 (MQ=255)
   ccTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTa                                                                                                                                                  <  2:4953896/149‑1 (MQ=255)
      cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATgg                                                                                                                                               <  2:931755/149‑1 (MQ=255)
      cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATgg                                                                                                                                               <  2:5165688/149‑1 (MQ=255)
       aGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATggg                                                                                                                                              >  1:4362645/1‑149 (MQ=255)
       aGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATggg                                                                                                                                              <  2:3539677/149‑1 (MQ=255)
          cATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAg                                                                                                                                           <  1:2572704/149‑1 (MQ=255)
          cATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAg                                                                                                                                           <  2:4333539/149‑1 (MQ=255)
           aTGGTTAAAATACACAGATGCAAACTGGTGTCCTCCTTATGTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTGTTTTTTCACCAAATCATGGCGATAATTCGCAAAATATAAGTTTCTTTGTGCAgg                                                                                                                                          <  2:89260/149‑1 (MQ=255)
               ttAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTgg                                                                                                                                      <  1:1010781/149‑1 (MQ=255)
                   aaTACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGtt                                                                                                                                  >  2:4163919/1‑149 (MQ=255)
                      aCAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAAt                                                                                                                               <  2:545456/149‑1 (MQ=255)
                       cAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATc                                                                                                                              >  2:1269895/1‑149 (MQ=255)
                        aaaGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCt                                                                                                                             >  2:1014840/1‑149 (MQ=255)
                        aaaGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCt                                                                                                                             >  1:3478251/1‑149 (MQ=255)
                         aaGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCtt                                                                                                                            <  2:1774049/149‑1 (MQ=255)
                            aTGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTggg                                                                                                                         <  2:2486842/149‑1 (MQ=255)
                            aTGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTggg                                                                                                                         >  1:1902812/1‑149 (MQ=255)
                               cAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGcc                                                                                                                      <  1:4538607/149‑1 (MQ=255)
                                   cTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGgc                                                                                                                  <  2:5085074/149‑1 (MQ=255)
                                     ggACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGgcgc                                                                                                                >  2:681491/1‑149 (MQ=255)
                                      gACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTTCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGcc                                                                                                               >  2:4400764/1‑149 (MQ=255)
                                      gACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGcc                                                                                                               <  1:3868027/149‑1 (MQ=255)
                                         gCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGc                                                                                                            <  2:1123899/149‑1 (MQ=255)
                                           ctcctGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGt                                                                                                          >  1:1552498/1‑149 (MQ=255)
                                               tGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTg                                                                                                      >  1:106148/1‑149 (MQ=255)
                                                 aaCTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAt                                                                                                    <  2:645649/149‑1 (MQ=255)
                                                 aaCTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAt                                                                                                    <  1:3547607/149‑1 (MQ=255)
                                                 aaCTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAt                                                                                                    <  1:5074662/149‑1 (MQ=255)
                                                  aCTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAtt                                                                                                   >  2:4631033/1‑149 (MQ=255)
                                                    tttGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                 >  1:1606615/1‑149 (MQ=255)
                                                       ggAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAg                                                                                              <  1:1148646/149‑1 (MQ=255)
                                                          aTCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTAc                                                                                           >  2:1363341/1‑149 (MQ=255)
                                                             aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGAt                                                                                        <  2:5841/149‑1 (MQ=255)
                                                              tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                       <  2:242011/149‑1 (MQ=255)
                                                                 tATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGt                                                                                    >  1:3981313/1‑149 (MQ=255)
                                                                  aTTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTc                                                                                   <  1:339537/149‑1 (MQ=255)
                                                                   ttGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCa                                                                                  <  1:4599267/149‑1 (MQ=255)
                                                                     gCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                <  2:2796474/149‑1 (MQ=255)
                                                                         gggCATTTCTTGACTACTATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAAc                                                                            <  2:4228855/149‑1 (MQ=255)
                                                                             aTTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATg                                                                        >  1:3002023/1‑149 (MQ=255)
                                                                             aTTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATg                                                                        <  1:3778669/149‑1 (MQ=255)
                                                                              tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATgg                                                                       <  1:1326127/149‑1 (MQ=255)
                                                                                   tgccaaCAATTTATATGCCTCTCATTCTGCACCACATCAGGGCGATACTTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCCTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCaaaa                                                                  <  2:33111/144‑1 (MQ=255)
                                                                                    gACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAg                                                                 <  2:3981313/149‑1 (MQ=255)
                                                                                     aCTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGc                                                                <  2:8259/149‑1 (MQ=255)
                                                                                     aCTACAATTTATATGCCTCTCATTAAGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGc                                                                <  2:8608/149‑1 (MQ=255)
                                                                                      cTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGcc                                                               >  2:2614729/1‑149 (MQ=255)
                                                                                      cTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGcc                                                               <  1:2289521/149‑1 (MQ=255)
                                                                                        aCAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAg                                                             <  2:91793/149‑1 (MQ=255)
                                                                                          aaTTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTGACAATGGCAAAAGCCAGGc                                                           <  1:1457015/149‑1 (MQ=255)
                                                                                          aaTTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGc                                                           >  2:2051662/1‑149 (MQ=255)
                                                                                                 atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAc                                                    <  2:3478251/149‑1 (MQ=255)
                                                                                                 atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAc                                                    >  2:5024689/1‑149 (MQ=255)
                                                                                                 atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAc                                                    >  2:2387849/1‑149 (MQ=255)
                                                                                                 atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAc                                                    <  1:1414226/149‑1 (MQ=255)
                                                                                                      tctcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGa                                               <  1:825433/149‑1 (MQ=255)
                                                                                                       ctcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAt                                              >  2:2662803/1‑149 (MQ=255)
                                                                                                         cATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAttt                                            <  2:4907538/149‑1 (MQ=255)
                                                                                                          aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTTCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAtttt                                           >  2:1636073/1‑149 (MQ=255)
                                                                                                          aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAtttt                                           <  1:1818401/149‑1 (MQ=255)
                                                                                                          aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAtttt                                           >  2:756679/1‑149 (MQ=255)
                                                                                                           ttATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                          >  1:5019400/1‑149 (MQ=255)
                                                                                                           ttATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                          >  1:642938/1‑149 (MQ=255)
                                                                                                           ttATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAGTTTc                                          >  2:4727595/1‑149 (MQ=255)
                                                                                                                caccacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGt                                     >  1:409123/1‑149 (MQ=255)
                                                                                                                 accacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGta                                    <  2:2262421/149‑1 (MQ=255)
                                                                                                                 accacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGta                                    >  2:3266648/1‑149 (MQ=255)
                                                                                                                   cacaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                  <  2:106148/149‑1 (MQ=255)
                                                                                                                   cacaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                  <  2:3590301/149‑1 (MQ=255)
                                                                                                                    acaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                 >  2:4593218/1‑149 (MQ=255)
                                                                                                                    acaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                 >  2:1002630/1‑149 (MQ=255)
                                                                                                                    acaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                 >  1:2663404/1‑149 (MQ=255)
                                                                                                                    acaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGAACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                 >  1:937432/1‑149 (MQ=255)
                                                                                                                     caTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattta                                >  1:978916/1‑149 (MQ=255)
                                                                                                                       tCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatttatt                              <  1:21228/149‑1 (MQ=255)
                                                                                                                       tCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatttatt                              >  1:1761069/1‑149 (MQ=255)
                                                                                                                         aGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTg                            <  2:1946569/149‑1 (MQ=255)
                                                                                                                            gCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGttt                         <  2:3027744/149‑1 (MQ=255)
                                                                                                                             cGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTc                        <  1:3166362/149‑1 (MQ=255)
                                                                                                                              gATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCg                       >  2:1820381/1‑149 (MQ=255)
                                                                                                                                tAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGcc                     >  2:1688708/1‑149 (MQ=255)
                                                                                                                                 aaTTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGccc                    <  1:1964049/149‑1 (MQ=255)
                                                                                                                                  aTTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCa                   >  1:5186313/1‑149 (MQ=255)
                                                                                                                                    tCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAAt                 <  2:1528602/149‑1 (MQ=255)
                                                                                                                                        aaaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCa             >  1:4253291/1‑149 (MQ=255)
                                                                                                                                         aaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAg            <  1:1269895/149‑1 (MQ=255)
                                                                                                                                            tataAGTTTCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaa         >  2:725633/1‑149 (MQ=255)
                                                                                                                                            tataAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaa         >  1:1110315/1‑149 (MQ=255)
                                                                                                                                              taAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatg       >  1:3065701/1‑149 (MQ=255)
                                                                                                                                                 gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaat    >  1:3517385/1‑149 (MQ=255)
                                                                                                                                                 gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaat    <  1:2839518/149‑1 (MQ=255)
                                                                                                                                                 gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaat    <  1:3260186/149‑1 (MQ=255)
                                                                                                                                                 gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaat    <  2:4031386/149‑1 (MQ=255)
                                                                                                                                                  ttCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg   >  2:304500/1‑149 (MQ=255)
                                                                                                                                                  ttCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg   >  2:4124475/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGC  >  CP009273/3181152‑3181447

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: