Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,708,861 C→A 6.0% E298D (GAG→GAT)  bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,708,8610CA 6.0% 903.6 / 21.9 240E298D (GAG→GAT) bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base C (124/101);  new base A (10/5);  total (134/106)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

CCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAG  >  CP009273/3708713‑3709010
                                                                                                                                                    |                                                                                                                                                     
ccAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAAtcg                                                                                                                                                      >  2:2315507/1‑149 (MQ=255)
  aaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATctct                                                                                                                                                    >  1:3882379/1‑150 (MQ=255)
    cATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTctgc                                                                                                                                                  <  2:1022954/150‑1 (MQ=255)
        gctgcCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGcc                                                                                                                                              <  2:4280178/150‑1 (MQ=255)
        gctgcCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGcc                                                                                                                                              >  2:3952809/1‑150 (MQ=255)
         ctgcCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGccc                                                                                                                                             <  1:3336945/150‑1 (MQ=255)
         ctgcCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGccc                                                                                                                                             <  1:2874080/150‑1 (MQ=255)
                   gACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGt                                                                                                                                   >  2:475858/1‑150 (MQ=255)
                    aCGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGtt                                                                                                                                  <  2:822169/150‑1 (MQ=255)
                    aCGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGtt                                                                                                                                  <  2:3782255/150‑1 (MQ=255)
                     cGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCATCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCTCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGttt                                                                                                                                 <  1:2542830/150‑1 (MQ=255)
                     cGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGttt                                                                                                                                 <  2:2668505/150‑1 (MQ=255)
                      gATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGtttt                                                                                                                                <  1:2224385/150‑1 (MQ=255)
                       atcatcCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTc                                                                                                                               <  1:1035762/150‑1 (MQ=255)
                         catcCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATATCTGCTGCCCATTCGGCGGTTTTCGc                                                                                                                             >  1:1578287/1‑150 (MQ=255)
                         catcCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATATCTGCTGCCCATTCGGCGGTTTTCGc                                                                                                                             >  1:3837295/1‑150 (MQ=255)
                          atcCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCt                                                                                                                            <  2:1727661/150‑1 (MQ=255)
                                tgtTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCg                                                                                                                      >  1:3520672/1‑150 (MQ=255)
                                  tttttGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTc                                                                                                                    >  1:1212657/1‑150 (MQ=255)
                                   ttttGCTCACCAAACTGTTGGCACTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCa                                                                                                                   <  2:1114608/150‑1 (MQ=255)
                                     ttGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCAct                                                                                                                 <  1:1101238/150‑1 (MQ=255)
                                        cTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACtctc                                                                                                              <  1:4047150/150‑1 (MQ=255)
                                        cTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACtctc                                                                                                              <  2:1822588/150‑1 (MQ=255)
                                        cTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCGGCCCATTCGGCGGTTTTCGCTATTCCGTCACtctc                                                                                                              <  2:60399/150‑1 (MQ=255)
                                          cACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGc                                                                                                            >  2:2148889/1‑150 (MQ=255)
                                           aCCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGcc                                                                                                           >  1:1021022/1‑150 (MQ=255)
                                           aCCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGcc                                                                                                           <  1:906920/150‑1 (MQ=255)
                                           aCCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGcc                                                                                                           >  2:362585/1‑150 (MQ=255)
                                             cAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCa                                                                                                         >  2:4072551/1‑150 (MQ=255)
                                             cAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATATCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCa                                                                                                         <  2:42535/150‑1 (MQ=255)
                                              aaaCTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAg                                                                                                        >  1:1375386/1‑150 (MQ=255)
                                              aaaCTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATATCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAg                                                                                                        >  2:4460410/1‑150 (MQ=255)
                                                  tGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCaaa                                                                                                    >  2:2967431/1‑150 (MQ=255)
                                                  tGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCaaa                                                                                                    >  1:1283207/1‑150 (MQ=255)
                                                     tGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTAGTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATaa                                                                                                 >  2:2473905/1‑150 (MQ=255)
                                                        cgcTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAagag                                                                                              <  1:1782236/150‑1 (MQ=255)
                                                         gcTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCCCTATTCCGTCACTCTAGCCCAGCAAATAAGAgg                                                                                             >  1:3672782/1‑150 (MQ=255)
                                                             cATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGaa                                                                                         >  2:4187661/1‑150 (MQ=255)
                                                             cATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGaa                                                                                         >  1:2142208/1‑150 (MQ=255)
                                                              aTTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAg                                                                                        >  1:743075/1‑150 (MQ=255)
                                                                tCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGAc                                                                                      <  2:1147243/150‑1 (MQ=255)
                                                                   ccAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGc                                                                                   <  2:2505412/150‑1 (MQ=255)
                                                                       ccTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATaa                                                                               >  1:415501/1‑150 (MQ=255)
                                                                         tGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAAcc                                                                             >  2:4169294/1‑150 (MQ=255)
                                                                         tGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAAcc                                                                             <  1:1317822/150‑1 (MQ=255)
                                                                          gCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGATGACGGCTTAACCt                                                                            <  1:301492/150‑1 (MQ=255)
                                                                          gCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCt                                                                            >  2:230035/1‑150 (MQ=255)
                                                                             aTCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCtgtg                                                                         >  2:2689503/1‑150 (MQ=255)
                                                                                aGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTgtggtg                                                                      >  2:2210380/1‑150 (MQ=255)
                                                                                aGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTgtggtg                                                                      >  2:3838952/1‑150 (MQ=255)
                                                                                aGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTgtggtg                                                                      >  1:982345/1‑150 (MQ=255)
                                                                                aGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTgtggtg                                                                      <  2:3363834/150‑1 (MQ=255)
                                                                                 gCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGc                                                                     >  2:1906108/1‑150 (MQ=255)
                                                                                  cATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCa                                                                    >  2:1833126/1‑150 (MQ=255)
                                                                                   aTGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCaa                                                                   >  1:950362/1‑150 (MQ=255)
                                                                                    tggtggTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAAc                                                                  <  1:1778477/150‑1 (MQ=255)
                                                                                    tggtggTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAAc                                                                  >  1:509749/1‑150 (MQ=255)
                                                                                     ggtggtATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAAcg                                                                 <  2:3745994/150‑1 (MQ=255)
                                                                                     ggtggtATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAAcg                                                                 >  1:659140/1‑150 (MQ=255)
                                                                                        ggtATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACgcgc                                                              >  1:334421/1‑150 (MQ=255)
                                                                                         gtATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGcc                                                             >  1:1420279/1‑150 (MQ=255)
                                                                                         gtATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGcc                                                             <  1:1866306/150‑1 (MQ=255)
                                                                                               tggtggAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAAc                                                       <  1:2605017/150‑1 (MQ=255)
                                                                                                    gAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCtgc                                                  >  2:3225818/1‑148 (MQ=255)
                                                                                                    gAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTtc                                                  >  2:1593197/1‑150 (MQ=255)
                                                                                                    gAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTtc                                                  >  2:2575173/1‑150 (MQ=255)
                                                                                                      aaaTAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTtcgt                                                >  2:280110/1‑150 (MQ=255)
                                                                                                       aaTAGCCGCCAGCTCGCGGGTTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTtcgtc                                               >  2:2247141/1‑150 (MQ=255)
                                                                                                       aaTAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTtcgtc                                               >  1:2439738/1‑150 (MQ=255)
                                                                                                             cgccAGCTCGCGGATTTTCGCTGCGCCAACACCAAAAATATCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACACTCGCCCAGCAAATAAGAGGCGAATACGGCATAACCTGTGGTGCAACTCGCCAGAAACGCTTCGTCGTGCCa                                         >  1:33238/1‑150 (MQ=255)
                                                                                                              gccAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCaa                                        >  2:1590455/1‑150 (MQ=255)
                                                                                                               ccAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAAc                                       <  2:210319/150‑1 (MQ=255)
                                                                                                               ccAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAAc                                       <  2:1923216/150‑1 (MQ=255)
                                                                                                                 aGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCa                                     >  1:657917/1‑150 (MQ=255)
                                                                                                                        ggATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTcca                              <  1:4484687/150‑1 (MQ=255)
                                                                                                                          aTTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTccacc                            >  1:4347869/1‑150 (MQ=255)
                                                                                                                            tttCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAg                          >  1:395936/1‑150 (MQ=255)
                                                                                                                            tttCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAg                          >  1:395931/1‑150 (MQ=255)
                                                                                                                             ttCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGc                         <  2:3101086/150‑1 (MQ=255)
                                                                                                                             ttCGCTGCGCCAACACCACAAATATCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGc                         >  2:3869645/1‑150 (MQ=255)
                                                                                                                              tCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCg                        >  1:3554205/1‑150 (MQ=255)
                                                                                                                                  tGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATg                    >  1:3190527/1‑150 (MQ=255)
                                                                                                                                  tGCGCCAACACCACAAATCTCTGCTGCCCAGTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATg                    <  2:3233003/150‑1 (MQ=255)
                                                                                                                                   gcgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATgg                   >  1:3647519/1‑150 (MQ=255)
                                                                                                                                    cgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATggg                  >  1:3739965/1‑150 (MQ=255)
                                                                                                                                     gcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGc                 >  2:3552740/1‑150 (MQ=255)
                                                                                                                                      ccAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCg                <  2:3126002/150‑1 (MQ=255)
                                                                                                                                       cAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGa               <  1:2157028/150‑1 (MQ=255)
                                                                                                                                         acacCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATc             <  2:3703934/150‑1 (MQ=255)
                                                                                                                                          caccacAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcc            <  2:2602898/150‑1 (MQ=255)
                                                                                                                                          caccacAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcc            <  2:3029545/150‑1 (MQ=255)
                                                                                                                                            ccacAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcccc          >  1:1625275/1‑150 (MQ=255)
                                                                                                                                            ccacAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcccc          >  1:1625231/1‑150 (MQ=255)
                                                                                                                                                  aTCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATc    >  2:1351890/1‑150 (MQ=255)
                                                                                                                                                   tctcTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCa   >  1:2278986/1‑150 (MQ=255)
                                                                                                                                                   tctcTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCa   <  2:768101/150‑1 (MQ=255)
                                                                                                                                                   tctcTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCa   <  2:768550/150‑1 (MQ=255)
                                                                                                                                                   tctcTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCa   <  2:3482696/150‑1 (MQ=255)
                                                                                                                                                    ctctGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAg  >  2:2595334/1‑150 (MQ=255)
                                                                                                                                                    ctctGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAg  >  1:1657225/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAG  >  CP009273/3708713‑3709010

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: ‑