Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 578,193 A→T 9.9% pseudogene (68/653 nt) ybcY ← pseudogene, DLP12 prophage; methyltransferase homology;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273578,1930AT 9.9% 641.9 / 43.0 216pseudogene (68/653 nt)ybcYpseudogene, DLP12 prophage; methyltransferase homology;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (98/96);  new base T (12/10);  total (110/106)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

GATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCGCAG  >  CP009273/578044‑578342
                                                                                                                                                     |                                                                                                                                                     
gATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCaaa                                                                                                                                                       <  1:1897345/150‑1 (MQ=255)
   atCAGACTACTCTCAGGTACTTGAGTAATGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAcc                                                                                                                                                    <  1:697296/150‑1 (MQ=255)
    tCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAccc                                                                                                                                                   >  1:3856831/1‑150 (MQ=255)
    tCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAccc                                                                                                                                                   >  1:1280341/1‑150 (MQ=255)
    tCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAccc                                                                                                                                                   >  2:1229599/1‑150 (MQ=255)
    tCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAccc                                                                                                                                                   >  2:3968662/1‑150 (MQ=255)
      aGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACccac                                                                                                                                                 >  2:744909/1‑150 (MQ=255)
      aGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACccac                                                                                                                                                 >  1:3990774/1‑150 (MQ=255)
       gACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACccacc                                                                                                                                                <  2:2685269/150‑1 (MQ=255)
        actactCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACATAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACccacca                                                                                                                                               >  2:3135127/1‑150 (MQ=255)
           actCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTc                                                                                                                                            <  2:3933902/150‑1 (MQ=255)
             tctcAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCAt                                                                                                                                          >  2:2221547/1‑150 (MQ=255)
                cAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATaaa                                                                                                                                       <  2:1115482/150‑1 (MQ=255)
                 aGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAg                                                                                                                                      >  2:573457/1‑150 (MQ=255)
                  ggTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGc                                                                                                                                     <  2:4886550/150‑1 (MQ=255)
                   gTACGTGAGTAAGGTAAAACCCAGGTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCt                                                                                                                                    <  2:1658443/150‑1 (MQ=255)
                    tACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCtt                                                                                                                                   <  2:3041175/150‑1 (MQ=255)
                         gAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAggg                                                                                                                              >  1:2172900/1‑150 (MQ=255)
                         gAGTAAGCTAAAACCCAGTTCCACCACCAATATCCAGATGGTTGTTACCTAACTGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCACAAAGCTTTAggg                                                                                                                              <  1:3817092/150‑1 (MQ=255)
                          aGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGt                                                                                                                             >  1:4180623/1‑150 (MQ=255)
                           gTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTa                                                                                                                            >  2:3046951/1‑150 (MQ=255)
                              aGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTGTAGGGTAAgt                                                                                                                         >  2:3902792/1‑150 (MQ=255)
                                 tAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAgtggt                                                                                                                      <  2:3230296/150‑1 (MQ=255)
                                  aaaaCCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAgtggtg                                                                                                                     >  2:567974/1‑150 (MQ=255)
                                   aaaCCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGgtgt                                                                                                                    >  2:4643391/1‑150 (MQ=255)
                                    aaCCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGgtgtg                                                                                                                   >  2:3244795/1‑150 (MQ=255)
                                      cccAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTa                                                                                                                 <  2:4852367/150‑1 (MQ=255)
                                       ccAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTaa                                                                                                                <  2:1492976/150‑1 (MQ=255)
                                       ccAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTaa                                                                                                                <  1:2594276/150‑1 (MQ=255)
                                          gTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAAtt                                                                                                             >  1:3121374/1‑150 (MQ=255)
                                          gTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGAGTTCCTTTGTAGGACATCCCCAAGCAAGCCGATTAGATACTCCCAAAACCCACCAGACAAAAAGCTTTAGTGTAAGTGATGAGTAAAtt                                                                                                             >  1:538706/1‑150 (MQ=255)
                                           ttCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTc                                                                                                            <  1:2944494/150‑1 (MQ=255)
                                             ccAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTa                                                                                                          <  2:1556032/150‑1 (MQ=255)
                                             ccAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTa                                                                                                          <  2:1054281/150‑1 (MQ=255)
                                              cAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAg                                                                                                         >  1:2080388/1‑150 (MQ=255)
                                                   ccAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCca                                                                                                    >  1:3652124/1‑150 (MQ=255)
                                                      atatCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCcatca                                                                                                 <  2:322461/150‑1 (MQ=255)
                                                        atCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCatcatc                                                                                               <  1:1086606/150‑1 (MQ=255)
                                                          ccAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCtg                                                                                             <  1:858703/150‑1 (MQ=255)
                                                          ccAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCtg                                                                                             <  2:646789/150‑1 (MQ=255)
                                                          ccAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCtg                                                                                             >  1:2142852/1‑150 (MQ=255)
                                                           cAGATTGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCtgt                                                                                            >  1:1421208/1‑150 (MQ=255)
                                                           cAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTCAGTGGTGTGTAAATTCTAGCCCCATCATCtgt                                                                                            >  2:1617755/1‑150 (MQ=255)
                                                           cAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCtgt                                                                                            >  2:2541890/1‑150 (MQ=255)
                                                              aTGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGtt                                                                                         >  2:4206966/1‑150 (MQ=255)
                                                              aTGGTTGTTACCTAAAGGGTCCAGAAAGTGTGGAAGAAGGTGGTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGtt                                                                                         <  2:2922772/150‑1 (MQ=255)
                                                               tGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGttt                                                                                        <  1:2369966/150‑1 (MQ=255)
                                                               tGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGttt                                                                                        >  1:5015826/1‑150 (MQ=255)
                                                                   tgttACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGttttttt                                                                                    >  2:2743752/1‑150 (MQ=255)
                                                                      tACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTttatt                                                                                 >  2:5158957/1‑150 (MQ=255)
                                                                       aCCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTttatta                                                                                <  1:1229599/150‑1 (MQ=255)
                                                                       aCCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTttatta                                                                                <  1:3968662/150‑1 (MQ=255)
                                                                        ccTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTaa                                                                               <  1:798971/150‑1 (MQ=255)
                                                                        ccTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTaa                                                                               >  2:2079690/1‑150 (MQ=255)
                                                                         cTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAAt                                                                              <  2:3164248/150‑1 (MQ=255)
                                                                          tAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAAtt                                                                             >  2:4740939/1‑150 (MQ=255)
                                                                          tAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAAtt                                                                             >  2:3790722/1‑150 (MQ=255)
                                                                               gTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCAcc                                                                        <  2:3313302/150‑1 (MQ=255)
                                                                                  ccAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATg                                                                     >  2:3580924/1‑150 (MQ=255)
                                                                                    aGAAAGTGTGGAAGAAGGTGTTCCTTTGTATGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAATCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGtt                                                                   <  1:2538602/150‑1 (MQ=255)
                                                                                    aGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGtt                                                                   >  1:4460014/1‑150 (MQ=255)
                                                                                     gAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTta                                                                  <  1:1732417/150‑1 (MQ=255)
                                                                                     gAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTta                                                                  >  1:714011/1‑150 (MQ=255)
                                                                                      aaaGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTtat                                                                 >  1:2805418/1‑150 (MQ=255)
                                                                                           gtgGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGttt                                                            <  1:332659/150‑1 (MQ=255)
                                                                                            tgGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGtttt                                                           >  2:4503610/1‑150 (MQ=255)
                                                                                            tgGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGtttt                                                           <  2:3574595/150‑1 (MQ=255)
                                                                                              gaagaaGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTAt                                                         >  2:4115972/1‑150 (MQ=255)
                                                                                                 gaagGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCGCCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTtg                                                      <  2:4923481/150‑1 (MQ=255)
                                                                                                  aagGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTtgt                                                     <  2:17262/150‑1 (MQ=255)
                                                                                                   agGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTtgtg                                                    >  2:1323301/1‑150 (MQ=255)
                                                                                                      tgtTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAAt                                                 >  2:2800024/1‑150 (MQ=255)
                                                                                                           cTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAAt                                            <  1:4440921/150‑1 (MQ=255)
                                                                                                             ttGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCa                                          <  1:573400/150‑1 (MQ=255)
                                                                                                                    aCATCCCCATGCAAGCCGATTTTATACTCCCAAAACCCACCAGTCATAAAGCTTTAGTGTAATTGGTGTTTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATgg                                   <  1:4798525/150‑1 (MQ=255)
                                                                                                                    aCATCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATgg                                   <  2:4963886/150‑1 (MQ=255)
                                                                                                                    aCATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATgg                                   >  1:3009107/1‑150 (MQ=255)
                                                                                                                     cATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGc                                  <  1:3046951/150‑1 (MQ=255)
                                                                                                                     cATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTGAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGc                                  <  2:2998901/150‑1 (MQ=255)
                                                                                                                       tCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCaa                                >  1:974465/1‑150 (MQ=255)
                                                                                                                         cccATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATg                              >  1:3570379/1‑150 (MQ=255)
                                                                                                                           cATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGaa                            <  1:2896116/150‑1 (MQ=255)
                                                                                                                           cATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGaa                            >  1:1542745/1‑150 (MQ=255)
                                                                                                                            aTGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAAt                           >  2:5057992/1‑150 (MQ=255)
                                                                                                                            aTGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAAt                           >  2:1022157/1‑150 (MQ=255)
                                                                                                                            aTGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAAt                           >  1:2914257/1‑150 (MQ=255)
                                                                                                                              gCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTa                         <  2:191524/150‑1 (MQ=255)
                                                                                                                                     gATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAgggg                  >  1:1075775/1‑150 (MQ=255)
                                                                                                                                         tGATACTCCCACAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTaa              <  1:649446/150‑1 (MQ=255)
                                                                                                                                         tGATACTCCCAATACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTaa              <  1:3019942/150‑1 (MQ=255)
                                                                                                                                              cTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAAtgct         <  2:4322202/150‑1 (MQ=255)
                                                                                                                                               tCCCAATACCCACCAGTCATAAAGCTTTAGGGGAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAAtgctg        <  2:1795824/150‑1 (MQ=255)
                                                                                                                                               tCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAAtgctg        <  1:1577447/150‑1 (MQ=255)
                                                                                                                                                 ccAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGcc      <  2:2049162/150‑1 (MQ=255)
                                                                                                                                                   aaaaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgc    >  2:4541646/1‑150 (MQ=255)
                                                                                                                                                   aaaaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgc    >  2:1971127/1‑150 (MQ=255)
                                                                                                                                                    aaaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgca   <  2:4656808/150‑1 (MQ=255)
                                                                                                                                                    aaaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgca   <  2:138336/150‑1 (MQ=255)
                                                                                                                                                    aaaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgca   >  1:4545401/1‑150 (MQ=255)
                                                                                                                                                     aaCCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCgcag  >  2:668445/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAATCAATTATGGCAATGAATTACAAGGGGTTAAATGCTGCCGCAG  >  CP009273/578044‑578342

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: