Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,275 G→T 22.2% H164N (CAT→AAT)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,851,2750GT22.2% 815.8 / 177.3 313H164N (CAT→AAT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (105/138);  new base T (35/35);  total (140/173)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

TCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCG  >  CP009273/3851127‑3851421
                                                                                                                                                    |                                                                                                                                                  
tagtagtaAAAATCAGCGTATTTCTCTAAATACTCCACCTGACAGGTGAACGTGTGATGCGGCTCATCATACGAATCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgc                                                                                                                                                   <  1:1907400/148‑1 (MQ=255)
tCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTAAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGAAgc                                                                                                                                                   >  2:586155/1‑150 (MQ=255)
 cGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgcg                                                                                                                                                  >  1:1197490/1‑150 (MQ=255)
   tagtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGcc                                                                                                                                                <  2:3730445/150‑1 (MQ=255)
    agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGccc                                                                                                                                               >  2:463447/1‑150 (MQ=255)
     gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCa                                                                                                                                              >  1:4031177/1‑150 (MQ=255)
      taAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAg                                                                                                                                             <  1:561407/150‑1 (MQ=255)
              aGCGTATTTCTCTTAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGt                                                                                                                                     <  1:184962/150‑1 (MQ=255)
              aGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGt                                                                                                                                     <  1:116990/150‑1 (MQ=255)
                cGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGttt                                                                                                                                   <  2:4483513/150‑1 (MQ=255)
                 gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                  <  2:4910323/150‑1 (MQ=255)
                    tttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtc                                                                                                                               >  2:4853842/1‑150 (MQ=255)
                     ttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtcg                                                                                                                              >  1:38164/1‑150 (MQ=255)
                     ttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtcg                                                                                                                              >  2:244704/1‑150 (MQ=255)
                     ttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtcg                                                                                                                              >  1:1118096/1‑150 (MQ=255)
                      tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                             >  1:3204848/1‑150 (MQ=255)
                      tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                             >  2:5023158/1‑150 (MQ=255)
                         ctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGatat                                                                                                                          <  1:1027162/150‑1 (MQ=255)
                          tAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                         >  2:1327523/1‑150 (MQ=255)
                            aaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGt                                                                                                                       >  2:2757398/1‑150 (MQ=255)
                             aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                      <  2:3263196/150‑1 (MQ=255)
                             aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                      <  2:3771998/150‑1 (MQ=255)
                             aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                      >  1:1406954/1‑150 (MQ=255)
                               aCTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                         <  1:1631342/145‑1 (MQ=255)
                                 tCCACCGGACAGGTGACCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGAGACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                  <  1:2027555/150‑1 (MQ=255)
                                 tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                  <  1:3109726/150‑1 (MQ=255)
                                  ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCt                                                                                                                 >  2:2023370/1‑150 (MQ=255)
                                    aCCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGt                                                                                                               >  1:4308763/1‑150 (MQ=255)
                                      cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                             >  1:690864/1‑150 (MQ=255)
                                      cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                             <  2:89898/150‑1 (MQ=255)
                                       ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                            <  1:2178967/150‑1 (MQ=255)
                                       ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                            <  2:816194/150‑1 (MQ=255)
                                        gACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAAt                                                                                                           >  1:3268869/1‑150 (MQ=255)
                                          cAGGTGAACGGGTGATGCGGCTCATCAGACGAACCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGGTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGGTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                         <  2:1432463/150‑1 (MQ=255)
                                               gAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGAc                                                                                                    >  2:1672191/1‑150 (MQ=255)
                                                  cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                 <  1:850356/150‑1 (MQ=255)
                                                  cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                 >  1:2387524/1‑150 (MQ=255)
                                                    ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                               <  2:4156186/150‑1 (MQ=255)
                                                      tGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGttt                                                                                             >  2:3039237/1‑150 (MQ=255)
                                                        aTGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAg                                                                                           <  2:922756/150‑1 (MQ=255)
                                                        aTGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGATCGCCTGTAAATCTTCGACGCTGTTTAg                                                                                           >  1:3767067/1‑150 (MQ=255)
                                                            ggCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCa                                                                                       >  2:4302905/1‑150 (MQ=255)
                                                             gCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCAt                                                                                      <  1:4710494/150‑1 (MQ=255)
                                                               tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                    >  1:34667/1‑150 (MQ=255)
                                                               tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                    <  2:2133511/150‑1 (MQ=255)
                                                                  tcatACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTgcgc                                                                                 >  2:1272316/1‑150 (MQ=255)
                                                                     tACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcac                                                                              <  2:1136574/150‑1 (MQ=255)
                                                                      aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                             <  2:4793822/150‑1 (MQ=255)
                                                                      aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                             <  2:2328208/150‑1 (MQ=255)
                                                                       cGAAACCACCATCAGGAAAGGCTCGTCTGCGTGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                            <  1:2890379/150‑1 (MQ=255)
                                                                         aaaCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGc                                                                          >  1:1132495/1‑150 (MQ=255)
                                                                          aaCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCt                                                                         <  1:369110/150‑1 (MQ=255)
                                                                            ccaccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTaa                                                                       <  2:574823/150‑1 (MQ=255)
                                                                            ccaccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTaa                                                                       <  2:3446007/150‑1 (MQ=255)
                                                                              accaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtc                                                                     <  2:1795859/150‑1 (MQ=255)
                                                                               ccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtct                                                                    <  1:3551951/150‑1 (MQ=255)
                                                                               ccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtct                                                                    >  2:4322523/1‑150 (MQ=255)
                                                                                   cAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttt                                                                >  2:1178501/1‑150 (MQ=255)
                                                                                     ggAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttttt                                                              >  2:210167/1‑150 (MQ=255)
                                                                                       aaaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTcc                                                            <  2:1363862/150‑1 (MQ=255)
                                                                                         aGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGt                                                          >  1:4255892/1‑150 (MQ=255)
                                                                                          ggCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTc                                                         >  1:248431/1‑150 (MQ=255)
                                                                                             tCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGt                                                      <  2:3099210/150‑1 (MQ=255)
                                                                                               gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                    >  1:336617/1‑150 (MQ=255)
                                                                                               gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                    <  2:4670386/150‑1 (MQ=255)
                                                                                               gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                    >  1:3036837/1‑150 (MQ=255)
                                                                                               gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                    >  1:2965096/1‑150 (MQ=255)
                                                                                                  tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                 >  1:4960159/1‑150 (MQ=255)
                                                                                                     gcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAg                                              <  1:2269619/150‑1 (MQ=255)
                                                                                                      cgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGa                                             <  2:4033779/150‑1 (MQ=255)
                                                                                                         gcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAg                                          >  1:769215/1‑150 (MQ=255)
                                                                                                          cAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGt                                         >  1:507377/1‑150 (MQ=255)
                                                                                                           aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                        >  1:2919860/1‑150 (MQ=255)
                                                                                                              cTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGc                                     >  1:2679941/1‑150 (MQ=255)
                                                                                                                   gCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCAt                                >  1:457475/1‑150 (MQ=255)
                                                                                                                   gCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCAt                                <  1:1043852/150‑1 (MQ=255)
                                                                                                                   gCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCAt                                >  2:1155607/1‑150 (MQ=255)
                                                                                                                    cAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCaac                               >  2:323649/1‑148 (MQ=255)
                                                                                                                      gAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGa                             <  2:5145388/150‑1 (MQ=255)
                                                                                                                        aaaTCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAc                           <  2:364115/150‑1 (MQ=255)
                                                                                                                            ccACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGt                       <  1:1557126/150‑1 (MQ=255)
                                                                                                                             cACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTa                      <  2:1074130/150‑1 (MQ=255)
                                                                                                                              aCCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTAGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTaa                     <  2:4954417/150‑1 (MQ=255)
                                                                                                                              aCCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTaa                     >  1:2698321/1‑150 (MQ=255)
                                                                                                                                 gCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                  <  1:3518166/150‑1 (MQ=255)
                                                                                                                                 gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                  >  1:4817199/1‑150 (MQ=255)
                                                                                                                                   ccGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGc                <  2:2583854/150‑1 (MQ=255)
                                                                                                                                    cGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCg               <  1:419145/150‑1 (MQ=255)
                                                                                                                                      aTTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc             <  2:44535/150‑1 (MQ=255)
                                                                                                                                        tGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTccc           <  1:1156697/150‑1 (MQ=255)
                                                                                                                                           tGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt        >  2:3200609/1‑150 (MQ=255)
                                                                                                                                           tGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt        >  2:2914525/1‑150 (MQ=255)
                                                                                                                                            gATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTc       >  1:2202346/1‑150 (MQ=255)
                                                                                                                                            gATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTc       <  2:2489710/150‑1 (MQ=255)
                                                                                                                                             aTACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTcc      <  2:2689510/150‑1 (MQ=255)
                                                                                                                                             aTACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCTGCGTCCCACTcc      <  2:1909763/150‑1 (MQ=255)
                                                                                                                                             aTACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTcc      >  2:219540/1‑150 (MQ=255)
                                                                                                                                              tACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcg     <  2:187837/150‑1 (MQ=255)
                                                                                                                                               aCGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcgg    <  1:3916372/150‑1 (MQ=255)
                                                                                                                                                 gATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggcg  <  1:4884121/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCG  >  CP009273/3851127‑3851421

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: